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In eternity without beginning, the splendor-ray of Thy beauty boasted 
Revealed became love; and, upon of the world, fire dashed. 

Prom that torch, reason wanted to kindle its lamp 
Jealousy's lightning flashed; and in chaos, the world dashed. 

The Persian Poet, Hafez (1325/26-1389/90) 
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Abstract 



Weak-scale supersymmetry is one of the most favoured theories beyond 
the Standard Model of particle physics that elegantly solves various theo- 
retical and observational problems in both particle physics and cosmology. 
In this thesis, I describe the theoretical foundations of supersymmetry, is- 
sues that it can address and concrete supersymmetric models that are 
widely used in phenomenological studies. I discuss how the predictions of 
supersymmetric models may be compared with observational data from 
both colliders and cosmology. I show why constraints on supersymmetric 
parameters by direct and indirect searches of particle dark matter are of 
particular interest in this respect. Gamma-ray observations of astrophysi- 
cal sources, in particular dwarf spheroidal galaxies, by the Fermi satellite, 
and recording nuclear recoil events and energies by future ton-scale direct 
detection experiments are shown to provide powerful tools in searches for 
supersymmetric dark matter and estimating supersymmetric parameters. 
I discuss some major statistical issues in supersymmetric global fits to 
experimental data. In particular, I further demonstrate that existing ad- 
vanced scanning techniques may fail in correctly mapping the statistical 
properties of the parameter spaces even for the simplest supersymmetric 
models. Complementary scanning methods based on Genetic Algorithms 
are proposed. 

Key words: supersymmetry, cosmology of theories beyond the Stan- 
dard Model, dark matter, gamma rays, dwarf galaxies, direct detection, 
statistical techniques, scanning algorithms, genetic algorithms, statistical 
coverage 
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Svensk sammanfattning 



Supersymmctri ar cn av dc most valstuderade teorierna for fysik bor- 
tom standardmodellen for partikclfysik. Den loser pa ett elegant satt 
flera teoretiska och observationella problem inom bade partikelfysik och 
kosmologi. I denna avhandling kommer jag att beskriva de teoretiska 
fundamenten for supersymmetri, de problem den kan losa och konkreta 
supersymmetriska modeller som anvands i fenomenologiska studier. Jag 
kommer att diskutera hur forutsagelser fran supersymmetriska modeller 
kan jamforas med observationella data fran bade partikelkoUiderarc och 
kosmologi. Jag visar ocksa varfor resultat fran dirckta och indirekta 
sokanden efter mork materia ar sarskilt intressanta. Observationer av 
gammastralning fran astrofysikaliska kallor, i synnerhet dvarggalaxer med 
Fermi-satcUitcn, samt kollisioner mod atomkarnor i kommande storskaliga 
direktdetektionsexperiment ar kraftfuUa verktyg i letandct cfter super- 
symmetrisk mork materia och for att bestamma de supersymmetriska 
parametrarna. Jag diskuterar nagra statistiska fragestallningar nar man 
gor globala anpassningar till experimentella data och visar att nuvarande 
avancerade tekniker for att skanna parameterrymden ibland misslyckas 
med att korrekt kartlagga de statistiska egenskaperna, aven for de enklaste 
supersymmetriska modellerna. Alternativa skanningsmetoder baserade 
pa genetiska algoritmer foreslas. 
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Preface 



This thesis deals with the phenomenology of weak-scale supersymmetry 
and strategies for comparing predictions of supersymmetric models with 
different types of observational data, in particular the ones related to the 
identification of dark matter particles. Currently, various experiments, 
either terrestrial, such as colliders and instruments for direct detection 
of dark matter, or celestial, such as cosmological space telescopes and 
dark matter indirect detection experiments, are providing an incredibly 
large amount of precise data that can be used as valuable sources of 
information about the fundamental laws and building blocks of Nature. 
Analysing these data in statistically consistent and numerically feasible 
ways is now one of the crucial tasks of cosmologists and particle physics 
phenomenologists. There are several issues and subtleties that should be 
addressed in this respect, and dealing with those form the bulk of the 
present work. 

The papers included in this thesis can be divided into two general 
categories: Some (Paper I and Paper III) mostly aim to illustrate 
how real data can be used in constraining supersymmetric and/or other 
fundamental theories, and others (Paper II and Paper IV) are more 
about whether existing statistical and numerical tools and algorithms 
are powerful enough for correctly comparing theoretical predictions with 
observations. 

Thesis plan 

This thesis is organised as follows. It is divided into three major parts: 
Part I is an introduction to the theoretical and statistical backgrounds rel- 
evant to my work. Part II summarises the main results we have obtained 
in our investigations and Part III presents the included papers. Part I is 
itself divided into 7 chapters: Chapter 1 is a short and non-technical in- 
troduction to the field and the main motivations for investigating models 
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of physics beyond the Standard Model of particle physics in particular 
supersymmetry, Chapters 2 and 3 discuss the motivations for considering 
supersymmetry as a possible underlying theory of Nature in more detail 
and in demand for explaining both the dark matter problem in cosmol- 
ogy and theoretical issues with the Standard Model, Chapter 4 introduces 
supersymmetry and its theoretical foundations in a top-down approach 
and in a rather technical language, Chapter 5 details the most interesting 
supersymmetric models that are being used in current phenomenological 
studies. Chapter 6 provides a review of different observational sources 
of information that can constrain supersymmetric models and parame- 
ters, and Chapter 7 describes statistical frameworks and techniques for 
analysing supersymmetry. 

Almost all the included papers are written in rather comprehensive, 
self-contained and self-explanatory manners. Therefore, in order to avoid 
any unnecessary repetitions, I have written the introductory chapters such 
that they provide in a rather consistent and coherent way a more general 
and detailed description of the field to which the papers contribute. This 
also provides some additional background material that may not have 
been discussed in detail in the papers. The reader is therefore strongly 
recommended to consult the papers for more advanced and technical dis- 
cussions. 

Contribution to papers 

Paper I focuses on potential experimental constraints one may place 

upon supersymmetric models from indirect searches of dark matter (this 
has been done for the particular case of the Constrained Minimal Su- 
persymmetric Standard Model (CMSSM) as the model, and gamma-ray 
observations of the dwarf galaxy Segue 1 as the data). We have assumed 
that the lightest neutralino is the dark matter particle that annihilates 
into gamma rays observable by our detectors. The instrument for obser- 
vations is the Large Area Telescope (LAT) aboard the Fermi satellite. 
Conventional state-of-the-art Bayesian techniques are employed for the 
exploration of the CMSSM parameter space and the model is constrained 
using the LAT data alone and also together with other experimental data 
in a global fit setup. In preparing and writing the paper, I was mostly in- 
volved in general discussions and edition of the manuscript. I also helped 
Pat Scott in setting up SuperBayeS for the numerical calculations. 

Paper II deals with the issue of efficiently scanning highly complex 
and poorly-understood parameter spaces of supersymmetric models. It 
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attempts to introduce a new scanning algorithm based on Genetic Algo- 
rithms (GAs) that is optimised for frequentist profile likelihood analyses 
of such models. In addition to comparing its performance with that of 
the conventional (Bayesian) methods and illustrating how our results can 
affect the entire statistical inference, some physical consequences of the 
results (in terms of the implications for the Large Hadron Collider (LHC) 
and dark matter searches) are also presented and discussed. The analy- 
ses are done for a global fit of the CMSSM to the existing cosmological 
and collider data. I have been the main author of the paper. The use of 
Genetic Algorithms for exploration of supersymmctric parameter spaces 
was to a great extent my own initiative. I modified SupeBayeS and added 
GA routines to it. I did the numerical calculations, analysed the results 
and produced the tables and figures. I wrote most of the text. 

Paper III aims to predict how far one can go in constraining super- 
symmetric models with future dark matter direct detection experiments. 
The methodology and the main strategy of the paper are very similar 
to the analysis of Paper I: The studied supcrsymmetric model is the 
CMSSM and nested sampling is used as the scanning technique. Both 
profile likelihoods and marginal posteriors are presented. I have been the 
main author of the paper, performed the numerical scans, analysed the 
results and produced tables and plots. Christopher Savage also signifi- 
cantly contributed to the work by providing the background material for 
direct detection theory and experiments, as well as preparing the likeli- 
hood functions for the experiments that I used in the analysis. 

Paper IV studies a rather technical issue in the statistical investi- 
gations of supcrsymmetric models, namely the coverage problem. The 
analysis of this paper was computationally very demanding and required 
a substantial amount of computational power; this made the project a 
rather lengthy and challanging one. I have been the main author for this 
paper as well. I wrote most of the text and produced the results and 
all plots and tables. The numerical likelihood function for the analysis 
was provided by Christopher Savage, but I performed all the scans and 
interpreted the results. 

Yashar Akrami 
Stockholm, April 2011 
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Chapter 1 



Why dark matter and why go 
beyond the Standard Model? 



The visible Universe that we know and love is made up of planets, stars, 
galaxies and clusters of galaxies. We know that these objects exist mostly 
because they emit light or other types of electromagnetic radiation which 
we detect cither by eye or by various telescopes. In addition, the celes- 
tial objects substantiate their existence through their gravitational effects 
which impact the motions of other objects in their vicinity. For most 
nearby astrophysical objects the two sources of information fairly agree 
and are therefore used as complementary ways in studying interesting 
properties of their sources. A problem emerges however when wc look at 
scales of the order of galaxies or larger, where the gravitational effects 
imply the presence of massive bodies that are not detected electromag- 
netically. These objects that exhibit all the gravitational properties of 
normal matter but do not emit electromagnetic radiation (and are there- 
fore invisible) are referred to as 'dark matter' (DM). 

Almost every attempt at explaining the nature of DM with the known 
types of matter has so far failed. This is mainly because DM seems to be 
required in order to consistently explain very different astrophysical phe- 
nomena that have been observed by completely different methods. This 
inevitably leads us to the assumption that DM is composed of new types 
of matter that are beyond our current understanding of the elementary 
particles and their interactions. 

Our present knowledge of the fundamental building blocks of the Uni- 
verse is summarised in the so-called Standard Model (SM) of particle 
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physics (for an introduction, see e.g. ref. [1]). The SM provides a mathe- 
matically consistent (though rather sophisticated) framework for describ- 
ing different phenomena in a relatively large range of energy scales. At 
low energies the model describes the everyday life processes in terms of 
normal atoms, molecules and chemical interactions between them, and at 
high energies it has been capable of explaining various processes observed 
in nuclear reactors, particle colliders and high-energy astrophysical pro- 
cesses with remarkably high precision. The SM is a quantum-mechanical 
description of particles (or fields) and is based on a particular theoretical 
framework called quantum field theory. 

The SM is now extensively tested at colliders and is in excellent agree- 
ment with the current data. However, as we stated earlier, the SM does 
not contain any type of matter with properties similar to the ones we 
need for DM. This simply implies that if DM exists, the SM has to be 
appropriately modified or extended so as to include DM particles with 
required properties. The need for DM is therefore one of the strongest 
motivations for going 'beyond' the SM. 

Apart from the lack of any DM candidates in the SM, there are ad- 
ditional reasons in support of the existence of new physics beyond this 
framework. These reasons are mainly motivated by some theoretically ir- 
ritating characteristics of the model that cannot be explained otherwise. 
Perhaps the most notorious one is that the SM does not contain grav- 
ity. Currently four different type of force have been known in Nature: the 
gravitational force between massive objects, the electromagnetic force be- 
tween charged particles, the strong force that put together neutrons and 
protons inside atomic nuclei, and the weak force which is responsible for 
radioactive processes. While three of these forces, i.e. electromagnetic, 
strong and weak are well described quantum mechanically by the SM, the 
gravitational interactions do not fit consistently into the model. The rea- 
son is that when one attempts to quantise gravity with the known math- 
ematical methods of quantum field theory, the resulting theory contains 
some infinities that cannot be removed in an acceptable manner. This 
is done for the other interactions through the so-called 'renormalisation' 
procedure, a method that breaks down for gravitational interactions. We 
are therefore forced to treat gravity as a classical field which is best de- 
scribed by Einstein's theory of general relativity. This distinction between 
gravity and the other forces does not lead to serious problems provided 
that we do not want to describe gravitational processes at high energies 
where the quantum effects become important. There are however inter- 
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esting high-energy cases where one needs to have a quantum-mechanical 
description of gravity so as to be able to study the physical systems. Two 
important examples are (1) extreme objects such as black holes and (2) 
the physics of the very early Universe. It is therefore commonly accepted 
that the SM must be modified at least at those high energies where gravity 
needs to be quantised. 

In addition, the SM possesses a very special mathematical structure 
that is based on particular types of fields and symmetries. This struc- 
ture, although being crucial for the model to successfully describe different 
phenomena in particle physics, does not find any explanation within the 
theoretical principles of the model. The model also contains some free 
parameters, such as masses and couplings whose values have been de- 
termined experimentally. Some of these parameters take on values that 
require extensive fine-tuning. All these aesthetically vexatious issues and 
a few more give us strong hints that the SM is not the fundamental de- 
scription of Nature and has to be appropriately extended. 

Fortunately, several interesting extensions for the SM exist, the best 
of which are those that address all or most of the aforementioned issues 
simultaneously. One of these proposals is weak-scale supersymmetry. It is 
a very powerful framework in which the SM is conjectured to be modified 
by some new physics that kicks in at energies just above the electroweak 
scale, i.e. the energy scale at which the electromagnetic and weak forces 
are assumed to be unified into one single electroweak force. This new 
physics assumes that all particles of the SM are accompanied by some 
partner particles that are more massive than the original ones. The exis- 
tence of these so-called superpartners provides elegant solutions to many 
of the problems listed above, and paves the way for the resolution of many 
others in some broader theoretical framework. An important example is 
the inclusion of new matter fields with properties similar to what we need 
for a viable DM candidate. 

Supersymmetric models, like any other theories in physics, need to be 
tested experimentally. Indeed, there have been many theoretically fasci- 
nating ideas in the history of physics that were abandoned only because 
they have not been consistent with particular experimental data. For- 
tunately, there are various sources of information from both man-made 
experiments, such as particle colliders, and astrophysical/cosmological 
observations that can be used for testing the supersymmetric models. 
Ideally, all these different types of data should be combined appropriately 
so as to give the most reliable answers to our questions about the validity 
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of particular models and frameworks. This is however not a trivial task, 
because there are usually various sources of complication and uncertainty 
that enter the game and, if not addressed properly, can make any interpre- 
tations completely unreliable. This is exactly where the main objectives 
of the present thesis stand. We would like to examine how a class of in- 
teresting supersymmetric models can be compared with observations in 
the presence of different experimental (and theoretical) uncertainties and 
statistical / numerical complications. 

First, in the following two chapters we give a more thorough (and 
more technical) description of the problems with the SM, including the 
need for DM. In each case, we describe in rather general terms how the 
problem finds appropriate solutions in supersymmetry. The detailed res- 
olutions of some of the problems will be discussed later when supersym- 
metry is defined and concrete supersymmetric models are presented in 
chapters 4 and 5, respectively. In chapter 6 we review important observa- 
tional constraints we have employed in our analyses and describe different 
uncertainties in each case. Chapter 7 will be devoted to a discussion of 
the main statistical and numerical issues that we have dealt with in our 
endeavour. In the last chapter, i.e. chapter 8, we will briefly review our 
major results and present an outlook for future work. 



Chapter 2 



The cosmological dark matter 
problem 



2.1 The standard cosmological model 

The standard model of cosmology (for an introduction, see e.g. refs. [2, 3]) 
is a mathematical framework for studying the largest-scale structures of 
the Universe and their dynamics. In other words, cosmologists attempt to 
answer various fundamental questions about the origin and evolution of 
the cosmos using the fundamental laws of physics. The model is based on 
Einstein's theory of general relativity as the currently best description of 
gravity at the classical level, as well as two important assumptions about 
the distribution of matter and energy in the Universe that are usually 
called together cosmological principles: the homogeneity and isotropy 
on large scales. The cosmological principles immediately imply that the 
correct metric for the Universe has to be of a particular form that is known 
as Priedmann-Lemaitre-Robertson- Walker (FLRW) metric and has the 
following form: 



ds'^ = df - a{tf 



r'^{de'^ + sin^Odct)^ 



(2.1) 



1 — kr"^ 

Here r, 6 and (p denote the spherical coordinates and t is time, a as 
a function of time, is called the scale factor of the Universe and is an 
unknown function that can be determined by solving the Einstein field 
equations 

G^u = SttGT^u- (2.2) 

Here G^j^ is the Einstein tensor which contains all geometric properties 
of spacetime and Tf^i, is the stress-energy-momentum tensor (or simply 
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stress-energy tensor) that includes the information about various sources 
of matter and energy on that spacetime. G is Newton's gravitational 
constant. The time evolution of a therefore depends upon the assumptions 
we make for the matter and energy content of the Universe, k is called 
the curvature parameter and depending on its value, the Universe may 
be closed, open or flat (corresponding to k = +1, k = —1 and A; = 0, 
respectively) . 

The assumption for the stress-energy tensor T^y on the right-hand 
side of Eq. 2.2 is that the matter and energy of the Universe can be well 
described by a perfect fluid that is characterised by two quantities p (its 
energy density) and p (its pressure) . By inserting the stress-energy tensor 
for such a fluid, T^,y = diag{p,p,p,p), into Eq. 2.2 we end up with the 
following simple equations: 



+ 



a 

a J a- 

a _ AttG 
a 3 



SvrG 



ip + 3p). 



(2.3) 



By solving these so-called Friedmann equations, one can obtain the 
dynamics of the Universe in terms of the time evolution of the scale factor 
a. The quantity d/a on the left-hand side of the first equation that gives 
the expansion rate is called Hubble parameter H{t). In order to solve 
Eqs. 2.3, it is essential to also know how p and p are related, i.e. what the 
equation of state (EoS) is for the perfect fluid. For normal non-relativistic 
matter, the energy density p is much larger than the pressure p and one 
can therefore reasonably assume that the EoS is simply pm = 0. For 
relativistic matter (or radiation) on the other hand pr = 3pr, and for the 
vacuum energy pv = —pv (vacuum energy can be effectively written in 
terms of a cosmological constant A in which case p^ = A/SvrG). 

In cosmology it is useful to write the various energy density contri- 
butions to the total density (at present time) in terms of the so-called 
density parameters O^, ^r, and for matter, radiation and vacuum, 
respectively. The same is usually done for the curvature term in the first 
Friedmann equation by defining ilk in an analogous way. These den- 
sity parameters are defined as the ratio of a density p at present time 
(/9°) to a specific quantity called the critical density pc- Pc (defined as 
Pc = 3Hq/8ttG, where Hq is the present value of the Hubble parameter) 
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Figure 2.1. Major change points in the history of the Universe. Credit: 
NASA/WMAP Science Team. 

is the density for which the Universe has an exact flat curvature: 

k A 

ilm = — , = —, = 772' = 7772- (^•^) 

Pc Pc oHq 

The first Friedmann equation in Eqs. 2.3 can be written in the following 
simple form in terms of the density parameters: 

H"^ = Hi {^ril + zf + ^m{l + zf - nk{l + zf + n^], (2.5) 

where z = a^/a — 1 is the redshift with oq being the present value of the 
scale factor usually taken to be 1. There are various ways to measure the 
Hubble parameter H as a function of time from which one can determine 
the values for different density parameters and therefore the energy budget 
of the Universe. 

Thanks to different high-precision cosmological observations, we have 
now been able to not only confirm the relative validity of our standard 
cosmological model, but also determine the values of different parameters 
that enter the mathematical formulation of the model to a high degree 
of accuracy. We now know that (see e.g. Fig. 2.1) the Universe started 
from an extremely hot and dense state about 13.7 billion years ago (a 
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state that we call the Big Bang) and then expanded, cooled down and 
became structured by galaxies, stars and other astrophysical objects. We 
also know that the curvature of the Universe is, to a good approximation, 
flat and also that it has recently entered an accelerated expansion phase. 
Although we still need a quantum theory of gravity to understand what 
exactly happened in the very early moments of the cosmic evolution, we 
have been able to infer some properties of the Universe at those times. For 
example there are various reasons to believe that shortly after its birth 
the Universe has seen a short inflationary phase during which its size has 
grown exponentially: (1) The Universe is (at least approximately) flat. 
(2) The observed cosmic microwave background radiation (i.e. the relic 
radiation from the recombination epoch at which photons that were orig- 
inally in thermal equilibrium with matter could escape the equilibrium 
and freely travel in the Universe) is to a great degree isotropic. (3) The 
Universe is not perfectly homogeneous and structures exist. All these fea- 
tures can be gracefully explained by inflation. The underlying mechanism 
for inflation is yet to be understood, but the evidence for its occurrence 
is so strong that it has now become one of the main paradigms of modern 
cosmology. 



2.2 The need for dark components 

Perhaps the best confirmation of our cosmological picture to date has 
been from observations of the cosmic microwave background (CMB). It is 
extremely difficult (if not impossible) to explain the black-body spectrum 
of the CMB with alternative cosmological models. The measurements 
performed by the NASA satellite Wilkinson Microwave Anisotropy Probe 
(WMAP) have played a central role in this direction [4]. Not only have 
such measurements confirmed the fact that the Big Bang theory is a suc- 
cessful description of the Universe, they have also determined the actual 
values of the density parameters we introduced in the previous section. By 
fitting the model to the so-called angular power spectrum of the CMB for 
the tiny temperature fiuctuations observed on the 7- year WMAP sky map 
(see e.g. Fig. 2.2), it is now known that, for example, = 0.27 ± 0.03 
and l^A = 0.73 ± 0.03. 

The first surprising observation is that the vacuum energy (or the 
cosmological constant) is non-zero and even constitutes about 74% of the 
total energy budget of the Universe. A similar number was for the first 
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Figure 2.2. Temperature fluctuations on the cosmic microwave back- 
ground (CMB) observed by the Wilkinson Microwave Anisotropy Probe 
(WMAP) sateUite (background image), and the angular power spectrum 
of the fluctuations (inset). Credit: NASA/ WMAP Science Team. 

time reported in 1998 by two different measurements of the so-called lu- 
minosity distance (a quantity that is defined in terms of the relationship 
between the absolute magnitude and apparent magnitude of an astronom- 
ical object and can be calculated theoretically for a cosmological model in 
terms of the Hubble parameter for an object with a specific redshift) us- 
ing Type la supernovae (SNe) [5, 6]. The first explanation for this energy 
component that implies a recent transition of the Universe to an acceler- 
ated expansion epoch was that it is just a cosmological constant. From a 
particle physics point of view, however, the vacuum energy density of the 
SM contributes to the cosmological constant and hence affects the expan- 
sion history of the Universe. But the value estimated in this way is much 
larger than the observed one and this poses a serious problem that cannot 
be explained within the SM [7]. It was then proposed that perhaps some 
new physics has made such contributions from the vacuum energy small 
(or zero) and what we observe cosmologically is not the cosmological con- 
stant but rather a new energy source (with an EoS parameter w = p/p 
that is not identically equal to —1) that can be detected only gravitation- 
ally (hence the name dark energy). There are numerous suggestions for 
the nature of the dark energy, most of which come from particle physics 
theories beyond the SM (for a review, see e.g. ref. [8]). 

Although the WMAP results imply that normal matter (with the EoS 
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Figure 2.3. A pie chart of the content of the Universe today. Credit: 
NASA/WMAP Science Team. 

of pm = 0) forms about 26% of the total energy density, the surprise comes 
from the value it has measured for the energy density of baryonic matter 
in the Universe. This is the matter that is composed mainly of baryons 
and includes all types of atoms we know. The baryons' energy density 
can be measured because the CMB angular power spectrum is sensitive 
directly to the amount of baryonic matter: While the location of the first 
peak (see Fig. 2.2) gives us information about the total amount of matter, 
i.e. r^m, the second peak tells us about the total amount of baryonic mat- 
ter 0;,. Estimations then determine Ofe to be 0.045ib0. 003. Comparing the 
values for 0^ and O;, indicates that the usual baryonic matter constitutes 
only about 4% of the energy content of the Universe and about 22% is 
non-baryonic (see Fig. 2.3). All baryons interact with photons and can be 
detected also through non-gravitational effects whereas the non-baryonic 
component has been detected only gravitationally and is therefore named 
dark matter. In order to agree with observations of large-scale struc- 
ture of the Universe, this non-baryonic dark matter must be dominantly 
cold (i.e. almost non-relativistic). This cold dark matter (CDM) together 
with the assumption that dark energy is nothing but the cosmological 
constant A, a hypothesis that is in excellent agreement with all existing 
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Figure 2.4. The concordance cosmological model: 68.3%, 95.4%, and 
99.7% confidence regions in the flm-^^A (left) and flm-w (right) planes 
determined by observations of Type la supernovae (SNe) , baryon acoustic 
oscillations (BAO) and cosmic microwave background (CMB). Adapted 
from ref. [9]. 

observations, contrives the foundations of our current standard model of 
cosmology that is accordingly called ACDM. 

The left panel of Fig. 2.4 shows the currently best constraints on the 
energy densities of matter and dark energy from three important types of 
cosmological observations, i.e. the CMB, Type la SNe and baryon acoustic 
oscillations (BAO) [9]. The latter refers to an overdensity of baryonic 
matter at certain length scales due to acoustic waves that propagated in 
the early Universe. BAO can be predicted from the ACDM model and 
compared with what we have observed from the distribution of galaxies on 
large scales. The right panel of Fig. 2.4 depicts constraints from the same 
set of data but in terms of O-m versus the equation of state parameter w for 
dark energy (if = —1 is for dark energy being the cosmological constant). 
By looking at both plots, it is quite interesting to see that the constraints 
from all these three sources of information are in perfect agreement with 
each other and also consistent with our theoretical model. This model is 
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Figure 2.5. An example of the rotation curves of galaxies (for NGC 
6503) where circular velocities of stars and gas are shown as a function 
of their distance from the galactic centre. Here, the dotted, dashed and 
dash-dotted lines are the contributions of gas, disk and dark matter, re- 
spectively. Adapted from ref. [18]. 

also in harmony with many other observations (such as constraints from 
Big Bang Nucleosynthesis (BBN) on the baryon density [10], gravitational 
lensing [11] and X-ray data from galaxy clusters [12]), and is accordingly 
called the concordance model of cosmology. 

The argument for the existence of dark matter, i.e. the mass density 
that is not luminous and cannot be seen in telescopes, is actually very old. 
Zwicky back in 1933 already reported the "missing mass" in the Coma 
cluster of galaxies by studying the motion of galaxies in the cluster and 
using the virial theorem [13]. A classic strong evidence for dark matter 
existing in the scale of galaxies comes from the study of rotation curves 
in spiral galaxies by Rubin [14, 15, 16, 17]. The observed rotation curves 
are not consistent with the standard theoretical assumptions unless one 
assumes the existence of dark matter halos surrounding all known contents 
of the galaxies, i.e. stars and gas (for an example, see e.g. Fig. 2.5). 

We should note here that some alternative explanations have been put 
forward that claim the anomalous observational data do not necessarily 
lead to the conclusion that dark matter exists. Some of these alternative 
proposals, such as the ones in the context of modified Newtonian dynamics 
(MOND) [19, 20], have been successful in for example explaining the 



2.2. The need for dark components 



15 




Figure 2.6. The colliding Bullet Cluster. This is a composite image that 
combines the optical image of the object with a gravitational lensing map 
(in blue) and X-ray observations (in pink). Optical data: NASA/STScI; 
Magellan/U.Arizona/[21]. Lensing map: NASA/STScI; ESO WFI; 
Magellan/U.Arizona/[21]. X-ray data: NASA/CXC/CfA/[22]. 



rotation curves of spiral galaxies (although in a rather ad hoc way). As we 
saw, the dark matter problem is not limited to astrophysical phenomena 
on particular scales and shows up in different observations from the scale 
of a galaxy to cosmological scales. It is in fact extremely difficult to 
explain all those observations without dark matter. 

Perhaps the best direct evidence for the existence of dark matter is 
the so-called Bullet Cluster[21] (see Fig. 2.6). The Bullet Cluster consists 
of two galaxy clusters that have recently collided. Fig. 2.6 is a composite 
picture that shows (apart from the optical image) two types of observa- 
tions of the cluster: gravitational lensing (in blue) and X-ray observations 
(in pink). Comparing these two cases evidently show that the baryonic 
gas component, which emits X-ray radiation, does not form the total mass 
of the cluster. Most of the mass, mapped by the lensing measurement, 
seem to come from a component that, in contract with the baryons, is 
collisionless: it does not interact with either baryonic gas or itself. These 
properties are all consistent with the assumption of dark matter. 
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2.3 Weakly Interacting Massive Particles 

The astrophysical/cosmological observations we discussed in the previous 
section all imply that dark matter probably exist. The next question we 
need to answer is what is the nature of dark matter, i.e. what are the 
basic constituents of it. We have already inferred some of the properties 
the dark matter components should possess: (1) They must be massive 
otherwise we would not have seen their gravitational effects. (2) They 
must be dark, i.e. they should not emit or absorb electromagnetic radi- 
ation (at least not noticeably), otherwise they would have already been 
detected by our telescopes. (3) They must be non-baryonic (confirmed 
by e.g. the observations of CMB anisotropics and BBN). (4) They must 
be effectively collisionless with respect to both normal matter and them- 
selves, otherwise they would loose energy through electromagnetic (or 
stronger) interactions and form dark matter disks (which contradicts the 
observations of galactic rotation curves). Observations of astrophysical 
systems like the Bullet Cluster could also not be explained in this case. 
(5) Dark matter must be cold(ish) (i.e. almost non-relativistic) , otherwise 
it would have not given rise to proper structure formation as we observe 
on cosmological scales. (6) It must be stable or at least very long-lived 
(compared to the age of the Universe) ; this is required because dark mat- 
ter comprises a significant fraction of the total energy of the Universe at 
the present time (this fraction is given in terms of the dark matter relic 
abundance VL^m)- 

Unfortunately, all attempts at finding a suitable dark matter candi- 
date in the framework of the SM of particle physics have so far failed. 
This is because there are no standard particles that can satisfy all the re- 
quirements we listed above, and this means that cosmology requires new 
particles. This takes us to the realm of particle dark matter, namely that 
dark matter is composed of some new particles that have not been discov- 
ered yet. The need for particle dark matter is one of the main motivations 
for us to go beyond the SM (for detailed introductions to particle dark 
matter, see e.g. refs. [23, 24, 25]). 

Fortunately, several viable dark matter candidates have been proposed 
in the literature (for a review, see e.g. ref. [26]) and most of the interesting 
ones fall into the class of Weakly Interacting Massive Particles (WIMPs) . 
WIMPs are particles that couple to the SM particles only through in- 
teractions that are of the order of the weak nuclear force (or weaker). 
This immediately tells us that WIMPs are electrically neutral, dark, ef- 
fectively collisionless and non-baryonic. They are also massive, usually 
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with masses within a few orders of magnitude of the electroweak scale. 
Having high enough masses also means that they are cold. WIMPs are 
also stable on cosmological timescales and this characteristic comes from 
a (usually imposed) discrete symmetry of the theory that gives WIMPs 
some conserved quantum number. This quantum number then prevents 
WIMPs from decaying into other particles and therefore makes them sta- 
ble. In most scenarios, WIMPs are produced thermally in the early 
Universe [27, 28, 29, 30]. A generic (and highly interesting) feature of 
thermally-produced WIMPs is that they naturally provide the correct 
relic density of dark matter {^dm), i-e. a value that is in excellent agree- 
ment with observations. We explain this intriguing feature in more detail 
below. 

In the early Universe, right after the Big Bang, all the created particles 
(including WIMPs) are in both chemical and thermal equilibrium. Here 
chemical equilibrium refers to the situation where the primordial particles 
are created and destructed with almost equal rates and no net changes in 
their abundances with time. On the other hand, by thermal equilibrium 
(which is also called kinetic equilibrium) we mean that the particles are 
in thermal contact with each other without a net exchange of energy. In 
this latter case the temperatures associated with the particles follow the 
global temperature of the Universe. 

Suppose that the number density associated with our hypothetical 
WIMP particles x is '^x' their relative velocity is v and they annihilate 
into lighter particles with the total annihilation cross-section a. The 
equation governing the evolution of the WIMP density is the Boltzmann 
equation [31] 

^ = _SHn^-{av){nl-nl,g), (2.6) 

where n^^eq is the equilibrium number density of the WIMPs, H is the 
Hubble parameter and the brackets (...) denote thermal average. For 
WIMPs with the mass m^, the equilibrium number density (in the non- 
relativistic limit) at the temperature T reads 

nx,eq = 9xi^f'e—, (2.7) 

where g^. is the number of degrees of freedom associated with the species 
X- 

A direct implication of Eq. 2.6 is that as long as the creation and 
annihilation of the WIMPs is larger than (or comparable with) the ex- 
pansion rate of the Universe (specified by the Hubble parameter), the 



18 



Chapter 2. The cosmological dark matter problem 



0.01 




1 10 100 1000 



x=m/T (time ->) 

Figure 2.7. Chemical freeze-out of WIMPs. Initially when the particles 
are in chemical (and thermal) equilibrium, their actual number density 
follows the equilibrium value Neq- At some later time, the particles fall 
out of chemical equilibrium (or freeze out) and their comoving number 
density becomes fixed. Adapted from ref. [31]. 

particles remain in chemical equilibrium. However, the Universe expands 
and cools, and this means that at some time and temperature, the inter- 
action rate drops below the expansion rate and the equilibrium can no 
longer be maintained. This process during which the WIMPs decouple 
from the other particles is called chemical 'freeze-out'. The number den- 
sity of such thermally-produced WIMPs at the end of chemical freeze-out 
determines the relic density of dark matter today. Obviously, the abun- 
dance of WIMPs at freeze-out (and consequently the dark matter relic 
density) depends on how large the annihilation cross-section is: Larger 
cross-sections cause the WIMPs to remain in chemical equilibrium for a 
longer period and therefore generate a lower relic density (see Fig. 2.7). 

Chemical freeze-out happens at a temperature Tp that for WIMPs 
with weak-scale masses is given approximately as Tp = m;^,/20 [31]. 
After chemical freeze-out, WIMPs still remain in thermal contact with the 
other particles for some time and kinetic freeze-out (or decoupling) hap- 
pens later. The temperature of the WIMPs before this time is the same 
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as the equilibrium temperature, and becomes fixed by kinetic decoupling 
afterwards. This means that the WIMPs will have a temperature lower 
than Tp after kinetic freeze-out and this makes the WIMPs move non- 
relativistically up to the present moment. This characteristic is crucial 
for WIMPs to be 'cold' dark matter. 

In order to obtain the relic density of WIMPs O^, one needs to solve 
Eq. 2.6 numerically. However, to a first-order approximation, it can be 
shown that under very general assumptions does not depend explicitly 
on the WIMP mass and only depends on its annihilation cross-section [23, 
31] in the following way: 

^ (av) 

where h = i?o/100kms~^Mpc~^ w 0.7. For weakly-interacting parti- 
cles with reasonable masses (i.e. with values close to the scale of the 
electroweak symmetry breaking), the quantity (av) can be estimated as 
(av) ~ a'^/m^, where a is the fine structure constant. Assuming a typ- 
ical value of ruy ~ 100 GeV for the WIMP mass, we obtain (av) w 
10~^^cm'^s~^. By inserting this value into Eq. 2.8, we obtain an approx- 
imate value for fi^ with the right order of magnitude. This interesting 
'coincidence', also often referred to as 'the WIMP miracle', means that, 
under the assumption of chemical freeze-out as the actual dark matter 
production mechanism occurred in the early Universe, any particles with 
generic properties of WIMPs can provide a dark matter relic density of 
the correct order. This particular characteristic of WIMPs makes them 
amongst the most interesting and popular dark matter candidates. 

There are a large number of WIMP dark matter candidates on the 
market proposed in different contexts [26], amongst which the lightest 
neutralino in supersymmetry [31, 32, 33], the lightest Kaluza-Klein par- 
ticle in models of Universal Extra Dimension (UED) [34] and the light- 
est inert scalar in the Inert Doublet Model (IDM) [35, 36] are the most 
widely-studied ones. The first one, i.e. the lightest neutralino provides 
arguably the leading dark matter candidate with almost all desired prop- 
erties. A substantial part of this thesis is devoted to the phenomenological 
aspects of the neutralino with particular emphasis on its implications for 
constraining models of weak-scale supersymmetry. 

Before we end this section, let us emphasize that although WIMP dark 
matter proves to be an extremely powerful idea that provides extensive 
scope for phenomenological studies of particle dark matter, there are a 



20 



Chapter 2. The cosmological dark matter problem 



number of other viable dark matter candidates that are either entirely 
non-WIMP or only WIMP-inspired. We do not intend to go through any 
of them here and just provide a list of the most interesting ones and refer 
the reader to the given references for detailed discussions (see also ref. [26] 
for a comprehensive review): axions [37, 38], gravitinos [39], axinos [40], 
sterile neutrinos [41], WIMPzillas [42], Minimal Dark Matter [43, 44], In- 
elastic Dark Matter (iDM) [45, 46], eXciting Dark Matter (XDM) [47], 
WIMPless dark matter [48, 49] and models with Sommerfeld enhance- 
ment [50, 51]. 



Chapter 3 



Theoretical issues with the 
Standard Model 

As stated earlier in chapter 1, the Standard Model of particle physics is 
currently the minimal mathematical description of all known matter parti- 
cles and their interactions that consistently explains various experimental 
observations, and holds over a wide range of energies. This includes phe- 
nomena that we observe in our everyday experiments (i.e. energies of the 
order of a few eV), as well as the ones that can be observed only at high- 
energy colliders and astrophysical processes (i.e. energies of ~ 100 GeV). 
The only key ingredient of this mathematical framework that still needs 
to be confirmed experimentally is the Higgs boson which is thought to be 
responsible for giving masses to the other particles. There are however 
alternative proposals for making the particles massive that although not 
excluded yet, are arguably less motivated (see e.g. refs. [52, 53] for one of 
the most competitive ones). Having said that, it became relatively man- 
ifest soon after its establishment in the 1970s that for purely theoretical 
reasons the SM is incomplete and probably not the end of the story. It 
therefore has to be modified or extended beyond certain energies (which 
are argued to be energies higher than TeV scales). 

As we discussed in the previous chapter, the need for a viable dark 
matter candidate is one pivotal reason for thinking about extensions of the 
SM. We advertised supersymmetry as one of the leading theories beyond 
the SM that provides such candidates. However, the nice thing with 
supersymmetry is that it also helps us circumvent many of the theoretical 
issues with the SM that are not related to the dark matter problem. 

Before we introduce supersymmetry and review supersymmetric mod- 
els, their properties and phenomenological implications in chapters 4, 5 
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and 6, we remind ourselves in this chapter of some of the most notable 
theoretical problems in the SM and corresponding arguments in support 
of the physics beyond the SM, in particular supersymmetry. Clearly with- 
out describing its mathematical foundations and concrete realisations in 
particle physics, we cannot discuss in detail how supersymmetry helps us 
address these problems. We will therefore come back to some of the issues 
raised here in chapter 5 and explain how they can be gracefully resolved 
in some interesting supersymmetric models. 

3.1 The gauge hierarchy problem 

In any quantum field theory, including the SM, all present parameters 
(such as masses and coupling constants) are affected by quantum radiative 
corrections. The amount of the corrections is generically a function of the 
cut-off scale that is used in the process of renormalising the theory or 
removing the divergences arising from various loop integrals. For the case 
of fermions (i.e. particles with half-integer spin) interacting with photons, 
the radiative corrections to the fermion masses nif have a logarithmic 
dependence on the cut-off scale A (here we use A for a Lorentz-invariant 
cut-off): 6mf oc m,j In (A/mj) (see e.g. refs. [54, 55, 56] for a detailed 
discussion). For gauge bosons (i.e. particles with spin 1 in the SM, such 
as photons), by using a gauge- invariant regulator (as is for example used 
in dimensional regularisation) , one can show that the radiative corrections 
to the masses vanish. The reason for the absence of linear, quadratic, or 
higher-order corrections to the masses of fermions and gauge bosons is 
known and attributed to the presence of some particular symmetries of 
the theory: chiral symmetry in the former case and gauge invariance in 
the latter. Such symmetries are said to protect the particle masses from 
large radiative corrections. 

The situation is however different for the scalar fields present in the 
theory, such as the Higgs boson of the SM. Restricting the discussion to 
the SM Higgs mass, the radiative correction to its mass from the self- 
interaction term in the SM Lagrangian reads 



which is quadratically divergent (i.e. when A increases to infinity, the term 
quadratic in A dominates over the others and Srnjj becomes infinitely 
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large). It should be noted that this is not the only quadratically diver- 
gent contribution to the radiative mass corrections for the Higgs boson: 
others come from gauge boson loops and fermion loops. An interesting 
feature of field theory is that the quadratically divergent contributions 
from the fermion loops have opposite signs relative to the contributions 
from the boson loops, an observation that, as we will argue below, plays 
an important role in one of our strongest motivations for extending the 
SM to its supersymmetric version. 

Since the SM is a renormalisable theory, there is in principle no prob- 
lem with the divergent radiative corrections to exist, because they can 
be absorbed into the so-called bare mass parameter. However, in an 'ef- 
fective field theory' interpretation of the SM (for an introduction, see 
e.g. ref. [57]), it is believed that the model is a valid description of parti- 
cle physics up to some particular energy scale which is characterised by 
the cut-off scale A. At energies beyond A, the SM may be modified by 
adding new degrees of freedom (i.e. new fields) that are associated with 
some heavy particles whose effects are neglected at low energies. One 
example of such modifications is the assumption that the gauge group of 
the SM (i.e. SU{?>)c x SU{2)l x U{1)y) is generalised to a larger grand 
unification group such as SU{5) or 50(10). A rather trivial value for A 
beyond which we expect new degrees of freedom to become important is 
the Planck scale Mp ~ 10^^ GeV, but A can certainly be as small as TeV 
scales beyond which the SM has not been tested yet. In this effective field 
theory framework, quadratically divergent corrections pose a theoretical 
problem. 

There are several reasons which indicate that the 'physical' Higgs mass 
(the mass that is measured experimentally) has to be no larger than a few 
hundred GeV. This is the total value after adding the correction given in 
Eq. 3.1 to the bare mass parameter of the theory. If A becomes very large, 
the quadratic term in Eq. 3.1 will dominate over the other terms and this 
effectively means that the physical mass is determined by the bare mass 
and the quadratic term. If one now assumes that the SM is valid below 
the scale of grand unification theories (GUTs) Mgut ^ GeV (see 
section 3.3), the required cancellation of the two large values implies that 
the bare Higgs mass parameter will have to be "fine-tuned" to 1 part 
in 10^^. This becomes even worse if A is as large as the obvious cut-off 
scale of Mp. This fine-tuning problem is often referred to as the 'gauge 
hierarchy problem' of the SM [58, 59, 60, 61]. In other words, the large 
quadratic corrections imply that if A >> 1 TeV, any predictions we make 
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Figure 3.1. An example of how supersymmetry solves the gauge hier- 
archy problem of the standard model. Quadratically divergent quantum 
corrections to the Higgs mass can be cancelled through the presence of 
equal numbers of fermion and boson loops that contribute equally but 
with opposite signs. Here the loop contributions are shown for the top 
quark t and its supersymmetric partner i. 



for physics at TeV energies are highly sensitive to the structure of the 
underlying high-energy theory with the SM being its effective incarnation 
at low energies. 

Although such a fine-tuning of the SM structure is mathematically 
allowed, it has been taken as a strong hint (although not necessarily^) 
that some new degrees of freedom must exist above the electroweak scale 
that 'naturally' cancel the problematic quadratic corrections in Eq. 3.1. 
These new degrees of freedom should then be soon revealed by TeV- 
energy experiments and observations both at colliders and in high-energy 
astrophysical phenomena. 

Weak-scale supersymmetry is arguably the leading proposal that pro- 
vides the required new degrees of freedom and solves the hierarchy prob- 
lem in a simple and elegant way. We mentioned earlier in this section 
that the fermion and boson loops contribute to the dangerous quadratic 
divergences with opposite signs. This immediately suggests that in a the- 
ory with equal numbers of fermionic and bosonic degrees of freedom, the 
quadratic divergences will be cancelled. In order for this idea to work 
at any loop level, the couplings of fermions and bosons are additionally 

^Examples of the alternative approaches include: (1) Simply accepting that Nature 
is actually fine-tuned. (2) Leaving the assumption that elementary scalar fields exist 
in Nature, in models with composite states of bound fermions such as the idea of 
technicolor [52, 53]. (3) Assuming that the Higgs bosons interact strongly (rather than 
perturbatively) with themselves, gauge fields or fermions at the cut-off scale A [62, 63]. 
(4) Making gravitational effects strong at energies close to TeV scales by for example 
assuming the existence of additional compact spatial dimensions [64, 65]. (5) Assuming 
that the quadratic divergences only show up at multi-loop level and not necessarily at 
the lowest order in models such as the Little Higgs [66]. 
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required to be related due to some symmetry. As we will see in the follow- 
ing chapters, both of these requirements are fulfilled in supersymmetry 
as a symmetry that transforms fermions to bosons and vice versa (see 
e.g. Fig. 3.1). 

As we will argue in section 4.5, even if supersymmetry is a correct ex- 
tension of the SM, it has to be broken at least spontaneously (i.e. through 
a mechanism similar to the Higgs mechanism of electroweak symmetry 
breaking). One can show that in a supersymmetric theory where su- 
persymmetry is appropriately broken, the scalar masses all remain sta- 
bilised against radiative corrections and the hierarchy problem is still 
resolved [67]. This observation is so remarkable that it essentially served 
as a watershed in the history of supersymmetry and provided one of the 
strongest motivations for it. 

3.2 Electroweak symmetry breaking 

Electroweak symmetry breaking (EWSB) is an essential ingredient in the 
SM. Through this process all the particles of the model acquire mass, 
a feature that is obviously a crucial requirement for the model to suc- 
cessfully describe the real world. EWSB is realised in the SM through 
the Higgs mechanism: The Higgs boson of the theory is believed to have 
acquired a vacuum expectation value (VEV) which results in the break- 
ing of electroweak gauge symmetry. This is a 'spontaneous' symmetry 
breaking, because the fundamental Lagrangian of the theory (i.e. the SM 
Lagrangian) still remains symmetric while the ground state is no longer 
invariant under the symmetry. In order for the Higgs boson to develop 
an appropriate VEV, a so-called scalar potential of the theory should be 
minimised properly. This requires some particular parameters of the po- 
tential to acquire specific values. Strictly speaking, in order for the EWSB 
mechanism to work, some squared mass parameter for the Higgs boson 
has to be negative and this can be achieved only if some parameters of 
the model possess certain values. Although these values have been set 
experimentally, there is no explanation for such choices and again some 
fine-tuning seems to be necessary. 

As we will discuss in sections 5.1.3 and 5.3.2 for particular supersym- 
metric models, supersymmetry can naturally lead to EWSB and provide 
a deeper understanding of why it happens. This is mainly because in 
supersymmetric models one usually does not have to tune the EWSB pa- 
rameters directly: The conditions of EWSB can be satisfied by setting the 



26 



Chapter 3. Theoretical issues with the Standard Model 



model parameters to some typical values that are motivated for other rea- 
sons. In models for supersymmetry with parameters that are set at some 
high-energy scales (such as the models of sections 5.3.1 and 5.3.2), start- 
ing from a few parameters and evolving them with energy by means of the 
so-called renormalisation group equations (RGEs; see e.g. section 5.1.6) 
can give rise to EWSB at the electroweak scale. This process is often 
referred to as 'radiative electroweak symmetry breaking' (REWSB). 

3.3 Gauge coupling unification 

The SM is constructed based on the gauge group SU{3)c x SU{2)i x 
U{1)y and all particles are different representations of this particular 
symmetry group. But why is this group special? It certainly looks pecu- 
liar and there is no theoretical explanation within the framework of the 
SM for this particular choice. 

The three subgroups of the above gauge group (i.e. 5?7(3)c, SU{2)l 
and C/(1)y), correspond to three forces of Nature, i.e. strong, weak and 
electromagnetic forces, respectively. Each group has a coupling constant 
that determines the strength of its associated force. Experimental mea- 
surements over a wide range of energies tell us that the three forces are 
very different in strength and this is related to the fact that the three 
corresponding coupling constants have very different values. Like any 
other quantity in quantum field theory that in general runs with energy, 
the couplings are also scale-dependent. However, experiments indicate 
that even at energies slightly higher than the electroweak scale where the 
spontaneously broken (sub-)symmetry SU{2)l x U{1)y becomes restored, 
the two associated coupling constants do not unify (see e.g. the left panel 
of Fig. 5.1 in section 5.1.6). 

The peculiar gauge structure of the SM has however important im- 
plications. For example, it prevents the occurrence of some unwanted 
phenomena such as proton decay and large flavour-changing neutral cur- 
rents (FCNCs). Although these characteristics are crucial for the success 
of the model, the way they are achieved in the SM is highly non-trivial 
and seems to be pure luck. In addition, the SM contains many free pa- 
rameters whose values are constrained by experiments. There is however 
no theoretical explanation for such experimentally favoured values. All 
these types of tuning problems, as well as the question about different 
values of gauge coupling constants find reasonable explanations through 
the intriguing idea of 'unification'. 
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In unified theories, the gauge symmetries of the SM are assumed to be 
extended to larger symmetries. For example in the so-called grand unified 
theories (GUTs), that are of particular interest in this respect, the SM 
symmetry group is extended to some simple Lie groups such as SU (5) [68] 
or 50(10) [69, 70, 71]. This extension is largely motivated by the fact 
that the SM field content perfectly fits into multiplets (or representations) 
of these groups, i.e. these larger groups include the SM group as their 
subgroup [72]. This can therefore potentially explain the reason for the 
particular assignment of quantum numbers (such as hypercharges) in the 
SM (which seem to be randomly assigned) . This consequently illuminates 
why dangerous experimental processes are forbidden in the SM. 

One requirement for unification to occur is that all gauge couplings 
of the theory unify to a single quantity. As we mentioned above, this 
is not the case for the SM. We however know that these couplings, as 
well as all other parameters of the model, generally evolve with energy 
through the RGBs. This then gives the hope that although the gauge 
couplings have different values at low energies, they may unify at some 
high energy scale where the underlying larger symmetry group manifests 
itself. If this scenario is true, it provides an appropriate answer to the 
question why different forces of Nature have different strength: this is only 
a natural consequence of running of parameters with energy in quantum 
field theory. In addition, as a bonus, unification usually provides extra 
relations between various parameters of the theory and therefore gives rise 
to a (sometimes dramatical) reduction in the number of free parameters 
of the model. This alleviates the problem with the large number of free 
parameters in the SM. Finally, promoting the peculiar gauge group of the 
SM to a simple group such as SU (5) or 50(10) is on its own an interesting 
feature. 

The problem manifests itself if we now solve the RGEs for the SM 
gauge couplings up to very high energies: the result is that the couplings 
do not unify at any scale (again see e.g. the left panel of Fig. 5.1 in sec- 
tion 5.1.6) and the idea of unification seems to be excluded. However, the 
unification scale (if exists) cannot be chosen arbitrarily and is determined 
by the particle content of the theory and measured values of different 
parameters at some energy scale (e.g. the weak scale). Although for the 
SM, with the known particle content and experimental constraints on its 
free parameters, the gauge couplings do not unify at any scale, a way 
out is to modify the particle content appropriately by adding new degrees 
of freedom to the model. Clearly these new particles should be heavy 
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enough so as to remain hidden at low energies. 

This is exactly where supersymmetry enters the game and turns out to 
be quite helpful. In most interesting versions of weak-scale supersymmet- 
ric models (as we will see in section 5.1.6) the SM field content is modified 
such that the gauge coupling unification can be elegantly achieved. In the 
case of minimal supersymmetric extensions of the SM, the unification is 
obtained typically at a GUT scale of Mqut — 2 x 10^^ GeV with a unified 
gauge coupling of gcvT — 0.7 (see e.g. the right panel of Fig. 5.1 in sec- 
tion 5.1.6) [94, 74, 75, 76, 77]. A detailed and technical discussion of why 
gauge unification is achieved in concrete realisations of supersymmetry 
can be found in section 5.1.6 of this thesis. 

Before we stop our discussion here, let us note that although unifi- 
cation can be obtained this way, it is however highly non-trivial from 
a theoretical point of view. For example, no firm theoretical explana- 
tions exist for why (7gut should remain in perturbative regime, or why 
the GUT scale Mqut resides in a narrow energy range that is required 
for both suppression of proton decay and prevention of possible quantum 
gravitational effects. These characteristics therefore remain as accidental 
properties of the theory. 

3.4 Experimental bounds on the Higgs boson mass 

In the SM, the Higgs mass is set by the quartic Higgs coupling A which is 
fairly unrestricted. Consequently, there are no strong limits on the Higgs 
mass and while the lower limit is set by experiments such as the Large 
Electron- Positron (LEP) collider to be about 114.4 GeV [78], the mass 
can be as large as about 800 GeV. 

On the other hand, as we will see in section 5.1.3, in the most widely- 
studied supersymmetric extensions of the SM, the Higgs mass is not a free 
parameter and is actually a prediction of the theory. In these models, the 
lightest Higgs scalar^ is required to be lighter than about 120 — 135 GeV 
and this much narrower range for the Higgs mass makes the theory more 
falsifiable and therefore phenomenologically more interesting. 

On the other hand, by fitting the SM parameters to the available elec- 
troweak precision data, the favoured value for the Higgs mass (i.e. the 
minimum-x^ point) is fairly low and well below the experimental direct 



^As we will see later, consistency conditions require supersymmetric models to have 
more than one Higgs boson. 
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Figure 3.2. Experimental constraints on tlie Higgs mass. Tiie region 
excluded by direct search limits from the Large Electron-Positron (LEP) 
collider is shown in yellow. The blue band represents the results of fitting 
the Standard Model parameters to the electroweak precision data. The 
red band depicts a similar fit when the SM is minimally extended to its 
supersymmetric version. Adapted from ref. [79] 

limit from the LEP (see Fig. 3.2). Although this discrepancy is not statis- 
tically very significant, it definitely shows some tension. For comparison, 
Fig. 3.2 shows also the result of a typical supersymmetric fit using one 
of the simplest supersymmetric extensions of the SM called the CMSSM 
(see section 5.3.2) [79]. It can be seen from this example that it is possible 
to reconcile theoretical predictions for the Higgs mass with experimental 
data within the supersymmetric extensions of the SM. 

3.5 The need for quantum gravity 

As stated in chapter 1, the SM of particle physics as a framework for 
describing the matter components of the Universe and their interactions, 
has been able to provide such a description in a mathematically consis- 
tent way only for three fundamental forces (out of four). The SM as 
a quantum-field theoretical framework is renormalisable only if gravita- 
tion is not included. Arguably, string theory (for an introcuction, see 
e.g. ref. [80]) has so far been the most favoured candidate for a consistent 
quantum theory of gravitation which is expected to include the SM as its 
effective field theory valid at low energies. It is however highly difficult to 
build a phenomenologically successful string theory that does not require 
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supersymmetry, and this means that supersymmetry is an essential ingre- 
dient of the best quantum description of gravitation so far. This makes 
supersymmetry particularly interesting. 

Even if we do not believe in string theory as a valid description of 
high-energy phenomena, there is yet another intriguing connection be- 
tween supersymmetry and gravity. As will be seen in the following chap- 
ters, the phenomenologically interesting versions of supersymmetry that 
we will consider are all based on a 'global' symmetry. In the language 
of chapter 4 this means that the generators of supersymmetric trans- 
formations are not functions of space and time. It is however entirely 
justified to promote the global symmetry to a local one, in a way anal- 
ogous to the gauge symmetries of the SM. Such a localisation process is 
shown to inevitably lead to the existence of a new spin-2 massless gauge 
field together with its supersymmetric partner, a spin-3/2 particle (see 
e.g. ref. [54]). The former is exactly the particle that is assumed to be 
responsible for gravitational interactions, and is accordingly called the 
graviton. The interesting characteristic of the supersymmetric graviton 
is that its dynamics, which is entirely fixed by local supersymmetry, con- 
tains Einstein's general relativity as our currently best classical theory 
of gravity. Regarding this connection with gravity, local supersymmetry 
is often called 'supergravity' (or SUGRA). Although such a supergravity 
theory is not renormalisable'^, its natural connection to gravity should not 
be ignored. 

3.6 Other issues 

In addition to the issues with the SM we enumerated in the previous 
sections, there are a few other reasons to believe that the SM is not the 
complete theory of Nature. Most of these arguments are again purely 
theoretical (or aesthetic) in nature, but are still highly intriguing so that 
one cannot simply ignore them. We will not attempt to detail these other 
problems here and only list (or briefly introduce) a few interesting ones 
with some references for further reading. 

The first problem comes from the observations of neutrino oscillations. 
These indicate that neutrinos have small but non-zero masses. In the SM, 

^Strictly speaking, this statement may not be correct. There is a particular version 
of supergravity, caUed 'N = S supergravity' (see chapter 4 for the terminology), which 
is conjectured to be renormalisable [81]. As we will point out in the next chapter, these 
versions of supersymmetry are however phenomenologically not very interesting. 
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neutrinos are however massless and this directly imphes that the model 
must be extended so as to accommodate massive neutrinos. In order to 
avoid a 'fine-tuning' problem, this is usually done through the so-called 
'seesaw' mechanism that is generally implemented within the framework 
of grand unified theories discussed in section 3.3 (for a review of neutrino 
masses and mixing, see e.g. ref. [82] and references thein). 

The other problem that is again related to a fine-tuning within the SM, 
is the 'strong CP problem'. This deals with the fact that the quantum 
chromodynamics (QCD) sector of the SM, contrary to the electroweak 
sector, respects the CP-symmetry. This leads to a highly fine-tuned value 
for a parameter called 'vacuum angle' and denoted by 6 [37]. The strong 
CP problem finds natural resolutions in models of physics beyond the 
SM, in particular through the introduction of new particles called axions 
(these are the same particles as the axions we mentioned in section 2.3 in 
our list of viable dark matter candidates) [83] . 

Let us end this chapter by adding to our list of issues two other the- 
oretical speculations on the structure of the SM: (1) Why are there only 
three generations for matter particles, i.e. for leptons and quarks? (2) 
What is the origin of fermion masses? These two may also find appropri- 
ate answers in theories beyond the SM. 
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Theoretical foundations of 
supersymmetry 

In the previous chapters, we attempted to review some answers to the 
question why we are interested in physics beyond the Standard Model 
of particle physics and in particular its supersymmetric extensions. Our 
discussions so far have been based on a very vague understanding of super- 
symmetry. Before we enter the world of concrete supersymmetric models 
in the next chapter and investigate various observational constraints on 
these models in chapter 5, we briefly introduce supersymmetry in this 
chapter and review some of its fundamental properties. In addition, some 
formalisms are discussed and basics of supersymmetric model building are 
presented. This chapter is a rather technical one and the reader who is 
only interested in phenomenological aspects of the field can simply skip 
it and continue directly from chapter 5. 

4.1 Supersymmetry is a symmetry 

All the known elementary particles are either bosons or fermions. Bosons 
are those particles (or fields) that obey Bose-Einstein statistics and this 
means that they can occupy the same quantum state at any given time. 
Fermions, on the other hand, obey Fermi-Dirac statistics and, conse- 
quently, only one fermion can occupy a particular quantum state at a 
time. Although the quantum mechanical distinction between matter and 
force is not a clear cut, fermions are often associated with matter whereas 
bosons are considered as carriers of forces and interactions between the 
fermions. According to the so-called spin-statistics theorem in quantum 
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field theory, bosons have integer spin while fermions possess half-integer 
spin (for an introduction to quantum field theory, see e.g. ref. [84]). 

Elementary fermions that are known to exist in Nature, according 
to the Standard Model of particle physics, are categorised as quarks (6 
particles u (up), d (down), c (charm), s (strange), t (top), h (bottom) 
and 6 corresponding antiparticles iZ, d, c, s, b) or leptons (3 charged 
particles e~ (electron), /^^ (muon), r~ (tau), 3 neutrinos Vf,, I'fj,, Vt and 
6 corresponding antiparticles e"*" (positron), r+, i>e-, i^fj., ^t)- The SM 
also contains 7 elementary bosons in total (if we include the graviton), 
some of which, such as the gauge bosons 7 (photon) , g (gluon) , and Z 
have already been discovered, while the other two, i.e. (Higgs boson) 
and G (graviton) are to be observed experimentally (see e.g. ref. [1] for 
an introduction to the SM). 

As we pointed out in the previous chapters, the mathematical struc- 
ture of the SM that describes its field content and various interactions 
between the fields, is constructed based on some particular symmetries, 
some of which are thought to be fundamental. 

The first symmetry from the latter category is called Pioncare sym- 
metry and is a 'spacetime' symmetry. The Poincare group (for an in- 
troduction to group theory and its applications in particle physics, see 
e.g. ref. [85]) is the full symmetry of special relativity and correspond- 
ingly any relativistic field theory; the SM is no exception. This is a 
10-dimensional noncompact Lie group and is the group of isometrics of 
Minkowski spacetime. The Poincare group includes the Lorentz group 
as a subgroup and is a semi-direct product of translations in spacetime 
and Lorentz transformations (i.e. R^'^ x 0(1, 3), where i?^'^ stands for the 
former and 0(1,3) denotes the latter). Mathematically speaking, all el- 
ementary particles (or fields) are different 'irreducible representations' of 
the group and are specified by two quantities: mass (or four-momentum) 
and spin (an intrinsic quantum number). The Poincare symmetry is con- 
sidered as a fundamental symmetry which every quantum-field theoretical 
framework that describes particles and their interactions should possess 
(including the SM and its potential extensions). 

The other symmetry that is implemented in the SM, and has been 
used as a guiding principle in constructing its theoretical structure, is an 
'internal' symmetry, i.e. a symmetry which is not obviously related to 
space and time. This determines how different components of a theory 
(e.g. different fields in the SM) transform into each other. The SM is 
called a 'gauge theory', and this is because the fundamental Lagrangian 
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of the theory is invariant (or symmetric) under a particular non-abelian 
gauge symmetry: SU (3)c x SU {2)lxU (l)y . This gauge symmetry is an 
example of internal symmetries. 

While spacetime symmetries of a quantum field theory dictate the 
properties of the field components and classify them into various cate- 
gories of scalars, vectors, tensors and spinors, internal symmetries rather 
determine how different terms in the Lagrangian must be written. The 
two symmetries are entirely independent in the SM. 

Supersymmetry or SUSY (for an introduction, see e.g. refs. [54, 55, 56, 
86, 87]) is a symmetry that transforms fermionic degrees of freedom into 
bosonic ones and vice versa. In a supersymmetric theory, every fermion 
has a bosonic 'superpartner' and every boson has a fermionic superpart- 
ner. With this definition, supersymmetry can be considered as a new 
internal symmetry because it gives a new way to transform some compo- 
nents of the theory, say fermions, to some other ones, i.e. bosons. This is 
however not entirely true. Supersymmetry is actually also an extension 
of the Poincare group in the sense that it extends the 'Poincare algebra' 
(and therefore special relativity) through the introduction of four anticom- 
muting 'spinor' generators. In other words, although in supersymmetry, 
fermions are transformed into bosons and vice versa, these transforma- 
tions only modify the particles' spin and this is essentially a spacetime 
property. 



4.2 The supersymmetry algebra 

Supersymmetry, as any other continuous symmetry, is characterised by 
a symmetry algebra, and as we mentioned in the previous section, this 
algebra is obtained by extending the Poincare algebra such that it relates 
two types of fields (bosons and fermions) in a single algebra. The resultant 
algebra is called a Lie 'superalgebra'. 

Let us first look at the Poincare algebra: Since the Poincare group is 
a semi-direct product of the Lorentz group and the group of spacetime 
translations, a general Poincare transformation contains both Lorentz 
transformations and translations. The Lorentz group has 6 generators: 
Jj; i = 1,2,3 (3 rotations) and Ki; i = 1,2,3 (3 boosts). One often 
denotes the generators of translations as P^; ^ = 0,1,2,3. In a more 
covariant looking form, the Lorentz generators are usually written as 
Mny = —Mnu, where Mqi = Ki and Mj,- = eij^Jk- In this notation, 
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the full Poincare algebra can be written as [87] 

[Pt.,Pu\ = 0, 

[M^y,Mpfj\ = igvpM^a - ig^pM^fj - iguaM^p + ig^^Mup, 

[Mp,u,Pp\ = -igp^^Pu + igpuP^- (4.1) 

In supersymmetry, the Poincare algebra is enlarged by generators that 
are 'spinorial', i.e. transform as spinors in contrast to the original 'ten- 
sorial' Poincare generators which transform as tensors. Such generators 
are often denoted by dotted and undotted spinors and Q^. These are 
objects that transform under the group Sl{2,C) as 

Qi ^ Q'i = M/Q^, Qi ^ Q'l = M*^Q'^, (4.2) 

where M G 5/(2, C) and a,/3 = 1,2. Here, the extra index / labels dis- 
tinct SUSY generators in case there are more than one pair. The number 
of such generator pairs is usually denoted by N (i.e. / = 1,...,A^). The 
simplest case with only one pair is accordingly called W = 1 supersym- 
metry'. 

Mathematically, there is no limit on N , but with increasing the 
theories contain particles of increasing spin. Since no consistent quantum 
field theory with spins larger than two exists, this then leads to the condi- 
tion N <8 [87].^ N = 1 supersymmetry (that is also called 'unextended' 
supersymmetry) is of particular interest since it is the only case which 
permits chiral fermions. We know that chiral fermions exist, therefore 
from a phenomenological point of view, any supersymmetric theory of 
particle physics has to be of = 1 type, at least at low energies. We 
therefore restrict our discussions to = 1 supersymmetry. 

The extended Poincare algebra (i.e. the superalgebra) includes the 
following new commutation and anti-commutation relations: 

[Pp,Qi]=0, [Pp,Qi]=0, 

[Mp,, Qi] = i{ap,)£Qj„ [iVV, Q^"] = i(a^,)".Q^^ 

{Qi, Qp = 2a^^.Pp6'J, {Qi, Qj} = e^pZ^^ , {Qi, Q|} = e .^(Z")* . 

Here, Z^-^ = -Z^^ are the so-called 'central charges' of the group, and 
they are the members that commute with all generators of the algebra. 

^This is however the case when gravity is part of the theory, otherwise, spins cannot 
be larger than one and this leads to A'^ < 4. 
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The = 1 (or unextended) SUSY algebra, is the simplest supersymmetry 
algebra which has no central charges. One important property of any 
supersymmetric theory that can be inferred from the above SUSY algebra 
is that the energy Pq is always positive. 

Any irreducible representation of the Poincare algebra is associated 
with a particle. Since the Poincare algebra is a subalgebra of the super- 
algebra, any representation of the latter is also a representation of the 
former. However in general, an irreducible representation of the superal- 
gebra corresponds to a reducible representation of the Poincare algebra, 
and this means that it corresponds to several particles. The particles of 
each SUSY representation are related to each other by the SUSY genera- 
tors Qi and Q^- This means that these particles have spins that differ by 
units of one half, i.e. some are bosons and some are fermions. The spin- 
statistics theorem then implies that the generators and transform 
fermions to bosons and vice versa. The particles that are obtained via 
supersymmetric transformations of other particles, are called supersym- 
metric partners or simply 'superpartners' of the original particles. 

An irreducible representation of supersymmetry that is equivalent to 
a set of supersymmetrically-related particle states is called a 'supermul- 
tiplet'. All particles (or states) belonging to a supermultiplet have equal 
masses, and any supermultiplet contains an equal number of fermionic 
and bosonic degrees of freedom. One can show that for = 1 (i.e. un- 
extended) supersymmetry (with gravity), only three types of massless 
multiplets exist: chiral multiplets (consisting of a Weyl fermion with spin 
1/2 and a complex scalar with spin 0), vector multiplets (consisting of 
a gauge boson with spin 1 and a Weyl fermion) and graviton multiplets 
(consisting of a graviton with spin 2 and a gravitino with spin 3/2). It is 
tempting to also add gravitino multiplets (consisting of a gravitino and 
a gauge boson) to the above set of multiplets, but such a multiplet can 
only happen in an extended (i.e. > 1) SUSY. 

4.3 The Wess-Zumino model 

In this section, we briefly describe how a supersymmetric field theory 
can be constructed at the 'action' level. Although the described model is 
too simple compared to the more sophisticated ones thought to be imple- 
mented in Nature (see chapter 5), it shows the main ingredients of any 
supersymmetric theories including the ones that are of phenomenological 
interest (see also refs. [54, 55, 56, 86]). 
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By looking again at the possible SUSY multiplets we enumerated in 
the previous section, we see that the simplest representation of supersym- 
metric transformations which includes a chiral fermion (with spin 1/2) 
is a chiral multiplet. In addition to a left-handed two-component Weyl 
fermion, that we denote by ip, this multiplet includes a complex scalar 
field, say (p. 

Let us now try to write down a four-dimensional SUSY-invariant ac- 
tion that is composed of M chiral multiplets with M scalar fields and 
M Weyl fermions (^ = 1,.-.,M). We demand the action to be su- 
persymmetric 'off-shell'. This means that the action is invariant under 
supersymmetry even if the classical equations of motion are not satisfied. 
The latter requirement leads to the addition of a set of 'auxiliary' complex 
scalar fields Fi (fields without kinetic terms) to the field content of the 
theory. The resultant Lagrangian density for our SUSY theory reads 

Aree = -8^" cj)* d^cf) + i^^a'' O^i^ + F* F. (4.3) 

This is a theory for massless and free fields (i.e. includes no interaction 
terms) and was first derived by Wess and Zumino [88]. The next step is 
obviously to add non-gauge interaction terms to the Lagrangian such that 
they preserve the supersymmetric property of the action. If we only retain 
the renormalisable interactions (i.e. ones with mass dimension < 4), it can 
be shown that the most general Lagrangian with non-gauge interactions 
for chiral multiplets has to have the following form: 

/:chiral = rfrcc + Ant = - W9^<Ai + 

_1 + - W'W*. (4.4) 

Here Wi and Wij are defined as derivatives of the so-called 'superpotential' 
W that is a function of the scalar fields 

W ^ ^. (4.5, 

= j-j-w. 

OCPiOCpj 

Here, M^^ is a symmetric mass matrix for the fermions, y*-'^ is a Yukawa 
coupling of a scalar and two fermions, and U are some additional pa- 
rameters that influence only the scalar potential of the Lagrangian [56]. 
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Auxiliary fields Fi are eliminated from the expression using their classical 
equations of motion. The term W^W* in Eq. 4.4 is only a function of the 
scalar fields (p and (j)* and is essentially the 'scalar potential' of the theory 
(usually denoted by V{4),(t)*)). The model introduced in Eq. 4.4 is called 
the Wess-Zumino model [89]. 

4.4 Supersymmetric gauge theories 

In the previous section, we showed how a simple supersymmetric theory 
looks like for a chiral supermultiplet. This can be used in constructing a 
model that describes particle physics fermions (i.e. leptons and quarks) 
and scalars (such as the Higgs boson). We however know that in reality, 
at least at low energies, there are other types of fields which should also 
be described in any supersymmetric extension of the SM: gauge fields. 

As pointed out in section 4.2, any SUSY theory for gauge fields should 
include vector (or gauge) multiplets as basic ingredients. These multiplets 
have massless gauge bosons (that we denote by j4^; = 0, 1, 2, 3), as well 
as Weyl fermions A". Here the index a can take on different integer values 
depending on the particular gauge group of the theory (e.g. a = 1, ...,8 
for 5?7(3)c, a = 1, 2, 3 for SU{2)l and a = 1 for C/(l)y). As for the chiral 
multiplet case, one has to also add an auxiliary field to the field content of 
the theory. Such a field, traditionally named D"^, is real and bosonic, and 
is required for the action to be SUSY-invariant off-shell. The Lagrangian 
density for the gauge multiplet is shown to have the following form: 

>Cgaugc = -\f;,F^''"' + iX^'^a^D^X'^ + ^D'^D'', (4.6) 

where F^^, is the Yang-Mills field strength for the gauge fields and 
D^X°^ is the covariant derivative of the A" field [56]. 

As the final step towards constructing a general supersymmetric Largan- 
gian, one needs to consider both contributions from the chiral and gauge 
supermultiplets, as well as any additional interaction terms that are al- 
lowed by gauge invariance and keep the theory supersymmetric. Adding 
the requirement that the interaction terms should be renormalisable (i.e. of 
mass dimension < 4 in four dimensions), our general Lagrangian density 
will have the following form: 

^ — •^chiral ~l~ ■^gauge 

-\/2(7((/.*r»A" - \/2(7A'f"('0"^r» + g{(t)*T''4))D'', (4.7) 
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where i2chirai and 



■gauge 



are defined in Eqs. 4.4 and 4.6, respectively. The 



only difference is that the ordinary derivatives in Eq. 4.4 for the chiral 
supermultiplet Lag rangian /^chiral 8b£Q now replaced by gauge-covariant 
derivatives L*^. T"' are the generators of the gauge group that satisfy 
\T°',T^] = if^-^'^T'^. Here f"'^'^ are the structure constants that define the 
gauge group and g is the 'gauge coupling'. 

The complete scalar potential of the theory in this case is shown to be 
expressible purely in terms of the auxiliary fields Fi and D"^ (which are 
in turn expressible only in terms of the scalar fields (pi): 



The first and second terms in Eq. 4.8 are called F-terms and D-terms, 
respectively. The former are entirely fixed by Yukawa couplings and 
fermion mass terms, while the latter are fixed by the gauge interactions. 
In addition, the scalar potential V{(f),(j)*) can be shown to be bounded 
from below, i.e. it is always greater than or equal to zero. 

Finally, one should notice here that the theory defined in Eq. 4.7 is 
invariant under 'global supersymmetry'. This means that the parameters 
of supersymmetric transformations do not depend on the spacetime co- 
ordinates. We however know that local symmetries also exist and some 
of them have played critical roles in our current description of particle 
physics: The best example is the gauge symmetries of the SM. 'Local 
supersymmetry' also exists and as we pointed out in section 3.5 makes 
an interesting connection between supersymmetry and gravity (in the 
context of 'supergravity'). The theory of supergravity is highly technical 
and we do not detail it in this thesis. We only briefly describe in the 
next chapter (mainly section 5.3) some phenomenologically interesting 
models that have been constructed based on supergravity assumptions. 
We refer the interested reader to the literature for detailed discussions 
(see e.g. refs. [54, 56] and references therein). 

4.5 Spontaneous supersymmetry breaking 

As we will argue in the next chapter, supersymmetry cannot be imple- 
mented in Nature as an exact symmetry and is required to be broken 
appropriately. From a theoretical point of view, arguably the most in- 
teresting way of breaking a symmetry in any quantum field theory is via 
a 'Higgs-like' mechanism, where the symmetry is broken 'spontaneously'. 




(4.8) 



a 
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This idea has seemingly worked very weh in the SM when the electroweak 
gauge symmetry is broken at TeV scales. It is therefore quite interest- 
ing to see how the same idea could work for supersymmetry (see also 
refs. [54, 55, 56, 86]). 

Spontaneous supersymmetry breaking means that while the Lagrangian 
of the theory is SUSY-invariant, the vacuum state |0) is not, i.e. Qa|0) 7^ 
and <5jjj|0) 7^ 0. In an unbroken supersymmetry, the vacuum has zero en- 
ergy (since -f^|0) = 0, where H is the Hamiltonian operator), while in 
a spontaneously-broken supersymmetry the vacuum has positive energy 
(i.e. (0|if|0) > 0). It can be shown from this that if the vacuum expecta- 
tion value (or VEV) of Fi and/or (the auxiliary fields introduced in 
sections 4.3 and 4.4) become non-zero (i.e. {D'^) 7^ and/or (Fi) 7^ 0), 
supersymmetry will be spontaneously broken. 

In SUSY-breaking models in which the vacuum state we live in is 
assumed to be the true ground state of the theory, the structure of the 
models usually imply that the equations Fi = and D'^ = cannot 
be satisfied simultaneously and this breaks SUSY spontaneously. Other 
models exist in which we are not assumed to live in the true ground state 
and instead live in a metastable SUSY-breaking state with sufficiently long 
lifetime (comparable to the current age of the Universe) (see e.g. ref. [90]). 
This metastable state might have been chosen by some finite temperature 
effects in the early Universe. 

Spontaneous SUSY breaking is usually implemented in different mod- 
els either through the Fayet-Iliopoulos (or 'Z)-term') mechanism [91, 92] 
or through the O'Raifeartaigh (or 'F-term') mechanism [93]. 

In the D-term SUSY-breaking mechanism, the gauge symmetry group 
needs to contain a ^7(1) subgroup with a non-zero D-teim VEV. Super- 
symmetry is then broken by introducing the following little extra piece to 
the SUSY Lagrangian: 

-^^Fayet-Iliopoulos = —kD, (4.9) 

which is a term proportional to D (k being a constant). 

In the i^-term mechanism, SUSY breaking occurs due to the existence 
of a non- vanishing F-term VEV that comes from a particular property of 
the superpotential W, namely that there is no simultaneous solutions for 
the equations 

Fi = -W* = 0, (4.10) 
with W* defined in Eq. 4.5 (i.e. W* = 6W*/6^*'). 
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One property of all types of spontaneous global SUSY breaking models 
(with stable or metastable vacuum states), is the existence of a massless 
neutral Weyl fermion as the Nambu-Goldstone mode. This fermion is 
called goldstino, is denoted by G and possesses the same quantum num- 
bers as the broken symmetry generator (which in our case is the fermionic 
charge Qa)- The G can be shown to have the form {{D"') /^/2, {Fi))'^ , 
i.e. its components are proportional to the VEVs of the auxiliary fields 
Fi and D" [56]. 

With this brief introduction to the two aforementioned SUSY-breaking 
mechanisms, we stop our discussion here. We will instead come back to 
the discussion of supersymmetry breaking in section 5.3.1 of the next 
chapter where we discuss various concrete scenarios in connection with 
phenomenologically interesting SUSY models. We will see how some of 
the general strategies described here can be used in constructing real- 
world theories. 



4.6 Superfield formalism 

In order to construct more complex supersymmetric Lagrangians, with 
larger numbers of fields and more complicated interaction terms, one 
needs to develop a rather general procedure that generates SUSY-invariant 
interactions in a systematic way. A compact and convenient way is to 
use 'superspace' and 'superfield' formalism. In = 1 supersymmetry, 
the superspace is the usual four-dimensional spacetime (labelled by the 
four coordinates x^; /i = 0, 1, 2, 3) enlarged by adding four anticommut- 
ing 'Grassmannian coordinates' Oa and 9a- These new coordinates are 
fermionic and transform as a two-component spinor and its conjugate. In 
general, for an extended supersymmetry with N SUSY generator pairs, 
there are 4A extra fermionic coordinates. Superfields are quantum fields 
that differ from the usual ones in that they are defined on the superspace 
rather than the spacetime. Superfields are defined as single objects with 
components being all the different fields (fermionic, bosonic and auxiliary) 
that belong to a supermultiplet. 

The main advantage of using superfield formalism is that the invari- 
ance under SUSY transformations remain manifest during the Lagrangian 
construction; this is because the Lagrangian is defined in terms of inte- 
grals over the superspace. Working with superfield formalism has also the 
advantage that the spacetime nature of supersymmetric transformations 
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is more manifest. Despite all the definite benefits of working with super- 
fields, the formalism is fairly complicated and we do not detail it here. 
We refer the interested reader to e.g. ref. [54] for a detailed introduction. 
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Chapter 5 

Supersymmetry in real life 



In the previous chapter, we described supersymmetry in general and dis- 
cussed various properties of a supersymmetric field theory. This was done 
mainly through the presentation of the simplest possible SUSY models 
with a minimal field content, i.e. the Wess-Zumino model and its gauge 
extension. However, these models are obviously too simple to describe 
the real world. It is the goal of the present chapter to discuss viable 
SUSY models and scenarios that may describe reality. We also argued in 
chapters 1, 2 and 3 why a supersymmetric extension of the SM is helpful, 
although our discussions were limited to rather general arguments. We 
may therefore want to see in a more explicit way how suprsymmetric mod- 
els could address the issues discussed there. We will detail in this chapter 
'some' of those issues in terms of definite SUSY models. Finally, in order 
to examine how observations could enhance our knowledge about super- 
symmetry, its validity and possible implementations in Nature (which has 
been the primary objective of this thesis), we need to have concrete the- 
oretical frameworks to work in. The present chapter also provides these 
frameworks. 



5.1 The Minimal Supersymmetric Standard Model 

There are various strategies in building a SUSY model that has to do 
with reality. In a top-down approach, one looks at some fundamentally 
motivated theories that accommodate supersymmetry, such as string the- 
ory. These theories are usually defined at very high energies that are not 
accessible by experiments. Phenomenological studies can then be carried 
out by extracting an effective supersymmetric field theory valid at low 
energies. 
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In an alternative bottom-up approach, one starts with the SM itself 
and adds all the ingredients that are required for it so as to become 
supersymmetric. It is important that in the latter approach one takes 
into account all phenomenological considerations and constraints in such 
a way that the emergent theory is consistent with observations as well 
as theoretical conditions. It must also give the SM as an effective theory 
valid up to certain energies since we know that the SM is an excellent 
description of particle physics below those energies. 

The simplest phenomenologically-constructed SUSY model (i.e. ob- 
tained through a bottom-up approach) is the so-called 'Minimal Super- 
symmetric Standard Model (MSSM)' [94] (see also refs. [54, 55, 56] for 
comprehensive introductions to the MSSM). It is minimal in the sense 
that it contains the smallest number of new particles (or fields) that can 
be added to the SM in order to make it supersymmetric, and the theory 
still remains consistent with all phenomenological requirements. In this 
section, we describe the MSSM and its properties that are of most interest 
for phenomenological studies of supersymmetry. 

5.1.1 Field content and superpotential 

In every supersymmetric model, including the MSSM, the number of de- 
grees of freedom for particles and their corresponding supersymmetric 
partners (or superpartners) match. This particularly implies that some 
SM particles have more than one superpartner. For example, the ele- 
mentary fermionic spin- 1/2 particles with two degrees of freedom (such 
as leptons and quarks) need two scalar superpartners with one degree of 
freedom each. The superpartners of the SM fermions are called 'sfermions' 
(sleptons for leptons and squarks for quarks) and the superpartners of the 
bosons are called 'bosinos' (gauginos for gauge bosons and Higgsinos for 
Higgs bosons). We also often refer to the superpartners of the SM parti- 
cles simply as 'sparticles'. 

In the MSSM, every known (i.e. SM) particle has a spin 0, 1/2 or 1 and 
must therefore reside, together with its superpartners, in either a chiral 
or gauge supermultiplet (see the previous chapter). We summarise in 
Tab. 5.1 all the particles and spartners in the MSSM. As we see, they are 
divided into two categories of chiral and gauge supermultiplets. Tab. 5.1 
also shows different hypercharges associated with the particles. These 
correspond to the three SM gauge groups. 

One interesting feature of the MSSM is that, contrary to the SM, 
it contains 'two' Higgs doublets (shown as (i/^, H^) and {H^, H^) in 
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Chiral supermultiplets 


Name 


Symbol 


spin 


spin 1/2 


SU{3)c 


5-17(2)^ 


U{1)y 


(s) quarks 
x3 




u 
d 


{UL dL) 


{ul dL) 
4 


3 
3 
3 


2 
1 
1 


1 

6 

2 
3 

1 
3 


(s)leptons 
x3 


L 

e 




{ly bl) 

'^R 


1 
1 


2 
1 


1 
2 

1 


Higgs(inos) 


Hu 
Hi 


{H+ Hi) 
{Hi H-) 


{H+ Hi) 
{Hi H~) 


1 
1 


2 
2 


+ 1 
1 
2 


Gauge supermultiplets 


Name 




spin 1/2 


spin 1 


SU{3)c 


SUi2)L 


C/(l)y 


gluon/gluino 




9 


9 


8 


1 





winos, W bosons 




W° 




1 


3 





bino, B boson 




B° 


B° 


1 


1 






Table 5.1. Chiral and gauge supermultiplets in the Minimal Supersym- 
metric Standard Model. The table is based on similar tables in ref. [56] 



Tab. 5.1) that consequently require two chiral supermultiplets. There 
are two main reasons for this: (1) Only one Higgs chiral supermultiplet 
would introduce a gauge anomaly in the electroweak gauge symmetry that 
would make the theory quantum-mechanically inconsistent. (2) The Higgs 
chiral supermultiplet that has the Yukawa couplings necessary for giving 
masses to the up- type quarks, has a hypercharge that is different from 
the hypercharge of the Higgs chiral supermultiplet that has the Yukawa 
couplings necessary for giving masses to the down-type quarks and the 
charged leptons (see e.g. ref. [56] for more details). 

The existence of two Higgs doublets in the MSSM and the fact that 
every bosonic degree of freedom has a corresponding fermionic degree of 
freedom and vice versa, together imply that the MSSM particle content 
is slightly more than a doubling of the SM particle content. It is also 
important to note that the sparticles presented in Tab. 5.1 are the 'in- 
teraction' (or gauge) eigenstates of the theory and the 'mass' eigenstates 
are in general linear combinations of the gauge eigenstates (we will come 
back to this in section 5.1.4). 

Like any other supersymmetric theory, the SUSY part of the MSSM 
Lagrangian is determined by a superpotential that is defined in terms of 
the chiral supermultiplets (see the previous chapter). The MSSM super- 
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potential is [56] 

Wmssm = uyuQHu - dydQHd - eyeLHa + fJ-HuHd , (5.1) 

where Q, L, u, d, e, Hu, and Hd denote chiral superfields of the theory. 
It is important to notice that there are three generations (i.e. three fam- 
ilies) for quarks/squarks and lepton/sleptons and although not written 
explicitly in Eq. 5.1, the summation over the generations is understood. 
Similarly, all gauge indices and summations are suppressed. The presence 
of three generations implies that the Yukawa couplings yu, yd and ye in 
Eq. 5.1 (which are exactly the same Yukawa couplings as those that enter 
the SM Lagrangian) are 3x3 matrices in the family space. 

The superpotential defined in Eq. 5.1 completely determines the struc- 
ture of the MSSM if SUSY is not broken (see the next section). This 
means that we have now obtained an exactly supersymmetrised version 
of the SM albeit at the cost of introducing one new parameter, i.e. fi. 

5.1.2 SUSY breaking and soft terms 

Simple phenomenological considerations imply that supersymmetry can- 
not be an exact symmetry of Nature (at least at low energies), and if 
implemented in Nature, must be broken spontaneously. In other words 
although the fundamental Lagrangian might be SUSY invariant, the vac- 
uum state that Nature has chosen need not be (see section 4.5 in the 
previous chapter). In a fully supersymmetric theory, masses of particles 
and their corresponding superpartners are equal. This immediately puts 
the theory in trouble if it is to describe reality. For example masses of 
selectrons (i.e. the superpartners of electrons) should be as low as the elec- 
tron mass, i.e. about 0.5 MeV. A particle with such a low mass should be 
easily detected experimentally, as electron is, and should essentially show 
up in our everyday life. This all means that SUSY is a broken symme- 
try. In addition, all sparticle masses should be much higher than the SM 
masses (in order not to have been observed in low-energy experiments). 

As we stated earlier, the MSSM is a phenomenological model in the 
sense that its general structure is not set by any fundamental high-energy 
theory. This clearly means that phenomenological considerations should 
also fix the structure of any SUSY-breaking terms that we may add to 
the MSSM Lagrangian (see e.g. refs. [54, 55, 56, 86, 95, 96]). 

One important guiding principle in determining the SUSY-breaking 
interactions in the MSSM comes from one of the strongest theoretical 
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motivations for extending the SM to its SUSY version, i.e. providing a 
solution to the gauge hierarchy problem. We argued in section 3.1 that the 
quadratic divergences from the radiative corrections to the scalar masses 
can be cancelled out in a supersymmetric theory if fermionic fields and 
their bosonic partners have equal masses. This is clearly not the case in a 
SUSY-broken theory. It can however be shown that if the sparticles have 
masses not much larger than TeV scales, the cancellation of the different 
loop contributions does not require huge fine-tuning and therefore SUSY 
can still provide a solution to the hierarchy problem [67] . 

Supersymmetry is broken in the MSSM by adding the so-called 'soft 
SUSY-breaking terms' to the exact supersymmetric Lagrangian described 
in the previous section. These are the terms that while breaking super- 
symmetry, satisfy four conditions: (1) They do not reintroduce quadratic 
divergences to the Higgs mass (i.e. the gauge hierarchy remains stabilised). 
(2) They preserve the gauge invariance of the SM (and correspondingly 
the SUSY-unbroken MSSM). (3) They do not violate the renormalisabil- 
ity of the theory (which can be achieved by adding only mass terms 
and coupling parameters with positive mass dimensions). (4) They re- 
spect baryon and lepton symmetries of the SM and therefore conserve 
the corresponding quantum numbers B and L. The most general soft 
supersymmetry-breaking Lagrangian then reads 

^^t^ = -\{M^gg + M2WW + MiBB + c.c) 

- {ua-a QHu - dad QHd - eag LHd + c.c.^ 

—Q^m^Q — L^m^L — uTm^li—drv^d — em|W^ 

- mjj^H^Hu - m^^H^Hd - (bH^Ha + c.c.) . (5.2) 



Here, Mi, M2 and M3 are bino, wino and gluino mass terms, respectively. 
Trilinear couplings au, aj and ae are complex 3x3 matrices in the family 
space and are in one-to-one correspondence to the Yukawa couplings yu, 
Yd and ye in the superpotential (see Eq. 5.1). Squark and slepton mass 
terms hiq, m^, m^, 111^ and m| are also 3x3 (Hermitian) matrices 
in the family space with potentially complex entries, rnj^^ and m'jj^ are 
explicit real mass terms in the Higgs sector and 6 is a complex bilinear 
coupling. 
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5.1.3 Electroweak symmetry breaking and Higgs sector 

In order for the Higgs mechanism to work, the scalar potential for the 
Higgs scalar fields needs be minimised and the minimum then breaks 
electroweak symmetry. It can be shown that at the minimum, both 
and can be set to 0, a property that is satisfactory. The reason for 
this satisfaction is that electromagnetism is not spontaneously broken at 
the minimum. Ignoring the terms in the potential that involve or 
HJ , one obtains the following expression for the Higgs scalar potential 
that only contains the neutral Higgs fields H*^ and H^: 

- (6i/X + c.c.) + ^(<7^ + 5'^)(|F°p-|i/°P)^. (5.3) 

Here g and g' are the SU(2)i and U{1)y gauge coupling constants, re- 
spectively, and h is the parameter defined in Eq. 5.2. 

Now, in order to break electroweak symmetry, Vh is required to be 
minimised (with a stable minimum) and the fields -ff^ and acquire 
real and non-zero vacuum expectation values (VEVs). We denote these 
VEVs by Vu and v^, i.e. 

One can show that in order for the scalar potential Vh to develop a 
well-defined local minimum such that electroweak symmetry is appropri- 
ately broken, the following two conditions must be satisfied [54]: 

> ("ik + l/^l')("^L + l/"l')' 
2\h\ > ml^+ml^ + 2\fi\\ (5.5) 

The existence of two Higgs doublets in the MSSM implies that the 
Higgs sector of the theory consists of eight degrees of freedom (two per 
each field for H^, H^, and H^). As in the SM, when electroweak 
symmetry is broken, three of these degrees of freedom are eaten so as to 
make W and Z bosons massive. This means that five degrees of freedom 
remain intact and form five physical Higgs scalars. They are usually 
shown as 

h,H,A,H^, (5.6) 

where h and H are CP-even and neutral (with h lighter than H), A is 
CP-odd and neutral, and are charged (with charges ±1). 
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The quantities Vu and in Eq. 5.4 are related to the masses of the 
Z-boson and M^-boson {mz and mw) and the gauge couplings g and g' 
in the following way: 

vl + vl = ^i% = '^. (5.7) 

In addition, the ratio of the two above VEVs parametrises how the 
total Higgs vacuum expectation value is divided between the two neutral 
Higgs scalars. This is an interesting quantity and is usually denoted by 
tan /3, i.e. 

tan/3=— . (5.8) 

In most phenomenological studies of the MSSM, all the 'effectively 
three' parameters of the superpotential rnj^^ + rn'jj^ + aiid b are 
usually traded for the three more physically interesting parameters (1) 

= v'^ + v^, (2) tan/3 and (3) the mass of one of the physical Higgs 
bosons (conventionally taken to be mA)- Furthermore, since the value of 
v"^ is fixed by the Z- and/or VF-boson masses (Eq. 5.7), one is left with 
only two parameters that need to be determined experimentally: tan /3 
and ruA- 

One main reason for this new parametrisation is that the EWSB con- 
ditions given in Eq. 5.5 can be written as [54] 

{m]j^+m]j^ + 2\^\'^)sm2p 



'''' - tan2 a - 1 ^' ^^-^^ 

the first of which allows us to trade b for tan /3, and the second can be 
used to fix the magnitude (but not the sign) of /i so as to obtain the 
measured value of mz- In addition, itla can be written in terms of the 
original parameters as [56] m\ = 26/ sin 2/3 = rn^H^ + ^^'h^ + 2|/^|^. 

The masses of the other physical Higgs bosons /i, H and can be 
obtained in terms of the above parameters: 



mi 



m\+m\±J {m\+m\Y —Am\m\{cos 2/3)^ 



m]j± = m\ + m^r. (5.10) 

These expressions are however valid only at tree level and for example 
large radiative corrections from the top squark/quark loops can elevate 
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to higher values (see e.g. ref. [97]). These corrections turn out to be rather 
crucial for the model to survive because for example the experimental 
bounds on the Higgs mass indicate that ruh should be larger than niz, 
a condition that is not consistent with what we obtain from the above 
relations: m/i < mz\ cos2/3| (see also section 6.2.2). 

5.1.4 Sparticles and their masses 

We briefly pointed out in section 5.1.1 that in the MSSM, the sparticles of 
Tab. 5.1 are in general different from the mass eigenstates of the theory. 
This is because after SUSY breaking and electroweak symmetry breaking, 
some particles share quantum numbers and can consequently mix. Except 
for the gluino that does not have quantum numbers similar to other parti- 
cles and therefore does not mix with them, the mixing happens for other 
gauginos, higgsinos, squarks and sleptons (gauginos and higgsinos can 
also mix with each other even though they belong to different groups) . In 
this section we review very quickly different mass mixing matrices within 
the MSSM and their corresponding eigenstates (see e.g. refs. [54, 55, 56] 
for more detailed discussions). 



Gluinos 

The simplest case to study is the gluino. Since it is a colour octet fermion, 
and SU{3)c is unbroken, it does not mix with any other MSSM particle 
and is a mass eigenstate. This all means that the gluino mass comes solely 
from the corresponding soft supersymmetry-breaking term of Eq. 5.2 

-^Msgg + C.C., (5.11) 
and the gluino mass is therefore equal to IM3I. 



Neutralinos 

After electroweak symmetry breaking, the neutral higgsinos {H^ and H^) 
and the neutral gauginos {B and W^) mix. The four mass eigenstates 
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corresponding to these fields are called 'neutralinos'. In the gauge eigen- 
state basis {B,W^ ,H^,H^)., the neutralino mass mixing matrix has the 
following form: 

/Ml -g'vd/V2 g'vu/^\ 

M2 gvdl^ -gVu/\/2 . , 

-g'vd/^ gvd/V2 -^x ■ ^ ■ > 

\ g'vu/V^ -gvu/y/2 -fi / 

Here Mi and M2 are the mass term parameters in Eq. 5.2, fj, is the SUSY 
higgsino mass parameter in the MSSM superpotential (Eq. 5.1) and g and 
g' are again the SU (2) ^ and U (l)y gauge coupling constants. The relation 
between g and g' is g' = gtanOw where 9w is the weak mixing angle. Vu 
and Vd are the neutral Higgs VEVs we introduced in section 5.1.3. 

In order to find the mass eigenstates, one diagonalises the mass matrix 
and the four eigenstates that are obtained from this are the four 
neutralinos Xi (^ = 1)2,3,4): 

X° = NaH^ + iV^2^° + N^sWl + ^¥^4^°. (5.13) 

in the above equation is the unitary matrix that diagonalises Mp. 
The lightest neutralino is of particular interest (see section 5.1.5) and 
is the particle that is usually meant by people when they speak of 'the 
neutralino'. It is often denoted by Xi or simply x- 



M70 = 



Charginos 

Following electroweak symmetry breaking, the charged higgsinos {H^ and 
H^) and winos {W~^ and W~) also mix and the mass eigenstates become 
linear combinations of these particles. We call the resultant mass eigen- 
states 'charginos'. 

The mass mixing matrix in this case and in the gauge eigenstate basis 
{W+, H+,W-, H-) reads 

M,. = i"^' ^M, (5.14) 

and can be diagonalised using two unitary 2x2 matrices U and V. The 
mass eigenstates that are obtained this way are 

x; = UiiW~ + Ui2H-, 

Xt = VaW+ + Vi2H^. (5.15) 
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Squarks and sleptons 

We write the squark squared-mass matrices M? and Mj in a basis where 
the squarks and their corresponding quarks are rotated in the same way. 
The squared-mass matrices for the sleptons and sneutrinos, M? and M?, 
respectively, are obtained in a similar way. The squared-mass matrices 
for the squarks and sleptons then read 



+ mutmu + D^^l uiu^ {Au'' - fi* cot P) 
" ~ V (au - /Ucot/3)mu + m^mj + D'^j^l 



^2 ^ / Vtm^ V + mdmd"^ mdt(ad"f - /i* tan /?) 
^ ^ (ad — /xtan/3)md + md^md + i^^^l 



M? 



hil + meme"'' + Dl^l me"''(ae"l' - fi* tan/3) 
(ae - /itan/3)me + mj nie + D^^j^l 



Ml = ml + 01^1- (5.16) 
Here, 

^LL ~ *^os 2/3(T3j — ej sin^ 6'vi/), 

^RR ~ cos 2/3ej sin^ 6*14^) (5-17) 

where T^f is the third component of the weak isospin, e/ is the electric 
charge and 6w is the weak mixing angle. V in Eqs. 5.16 denotes the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix. The matrices hiq, m^, 
m^, m£ and m| are the soft sfermion squared-mass matrices defined in 
Eq. 5.2, au, ad and ag are trilinear couplings, and fi is the /^-parameter 
in the MSSM superpotential (Eq. 5.1). niu, md and me are diagonal 
matrices defined in terms of the quark and charged-lepton masses as 

mu = diag{mu,mc,mt), 
Old = diag{md,ms,mb), 

me = diag{me,mfj,,mr)- (5.18) 

5.1.5 il-parity 

We discussed earlier in this chapter that the soft SUSY- breaking terms are 
chosen such that they conserve both baryon (B) and lepton (L) numbers. 
One immediate consequence of this property is that a discrete symmetry 
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exists between the SM particles and their superpartners. This symmetry 
has a corresponding multiphcative quantum number, known as i?-parity 
which is conserved and can be written in terms of B, L and the particle's 
spin s in the following form: 

i? = (-l)3(«-^)+2«. (5.19) 

It is not difficult to show that all SM particles have i?-parity +1 while all 
their superpartners have i?-parity —1. 

i?-parity conservation in the MSSM has an extremely interesting phe- 
nomenological consequence: the lightest supersymmetric particle (LSP) 
is stable and does not decay into lighter SM states (clearly it also does not 
decay into any other SUSY-particles since by definition it is the lightest 
member of that group). If the LSP is also weakly-interacting and electri- 
cally neutral, it can be a viable dark matter candidate (see e.g. section 2.3 
and references therein) . One example is the neutralino that we introduced 
in section 5.1.4. It is arguably the most favoured dark matter candidate 
and has received the bulk of attention to date. 

B and L conservation, and consequently i?-parity conservation, are 
not fundamental assumptions in the MSSM and there are models with 
i?-partity violation. However, there are good reasons to think that the as- 
sumption is not too ad hoc. For example some grand unified theories (see 
section 3.3), which give the MSSM as their effective low-energy versions, 
accommodate i?-parity conservation in their structure (see e.g. ref. [98] 
for certain 50(10) theories). 

5.1.6 Renormalisation Group Equations 

In the MSSM, analogous to any other quantum field theory, all the pa- 
rameters of the model are subject to running (i.e. evolving with energy 
scale). These include gauge coupling constants, parameters in the super- 
potential and soft SUSY-breaking terms. The running of the parameters 
can be calculated using the renormalisation group equations (RGEs). 

Gauge coupling constants 

Let us first look at the RGE evolution of the gauge coupling constants by 
introducing the commonly used couplings (i = 1, 2, 3): 



(5.20) 
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Figure 5.1. Running of tiie gauge coupling constants in botli tlie Stan- 
dard Model (SM) and its minimal supersymmetric extension (MSSM). 
The left panel depicts the failure of the SM gauge coupling to unify, while 
the right panel shows that the gauge coupling unification can happen in 
the MSSM. Adapted from ref. [99]. 



Here, gi = ^/5/3g' , 92 = 9 and 93 = 9s, where 9', g and gs are the C/(l)y, 
SU{2)l and SU{?>)c gauge couphng constants, respectively. To one-loop 
order, the RGEs for the couplings ai have the form 



d_ 



27r 



(5.21) 



where t = InQ and Q is the energy scale of the running, hi are some 
quantities that are fixed by the structure of the gauge group and the 
matter multiplets (to which the gauge fields couple) of the theory. It can 
be shown that for the SM 6f^^ = -41/10, 6f^^ = 19/6 and 6f^^ = 7, and 
for the MSSM hf^^^'^ = -33/5, b^^^^^ = -1 and ftf -5^^^ = 3 [55]. 

Eqs. 5.21 can be integrated to give the following equations in terms 
of the inverse of the couplings Oj: 



Q 



(5.22) 



Here Qq is the energy scale at which we start the running, and is usually 
taken to be mz since the couplings are well measured at the energy scale 
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mz- Taking into account the effects of two-loop corrections, and other 
subtleties, one can solve Eqs. 5.22 and the results are shown in Fig. 5.1 
for both the SM and the MSSM. 

It is interesting to see how gauge coupling unification can occur in the 
MSSM while it fails in the SM. This provides another explicit example 
of our general discussions in chapter 3, namely that extending the SM to 
its supersymmetric version provides appropriate solutions to some issues 
with the SM (see section 3.3 for details). The other interesting obser- 
vation is the scale Mqut at which the unification happens: it is about 
2.2 X 10^^ GeV, a value that is predicted as the unification scale in many 
grand unified theories but with entirely different motivations (again see 
section 3.3 and references therein). 

Superpotential parameters 

The parameters in the MSSM superpotential include the Yukawa coupling 
parameters yu. Yd) Ye, and the /i parameter (see Eq. 5.1). In order 
to calculate the RGEs for the Yukawa couplings, one usually makes an 
approximation that only the third- family components are important. This 
is because the third-family particles, i.e. the top quark, bottom quark 
and tau lepton, are the heaviest fermions in the SM. In other words, 
one assumes that the Yukawa couplings have the diagonal forms yu = 
diag{0,0,yt), yd = diag{0,0,yb) and ye = diag{0,0,yr). The one-loop 
RGEs for the parameters yt, yt and y-r, as well as the parameter /i are [56]: 
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(5.26) 



Soft SUSY-breaking parameters 

We now look at the RGEs for the parameters in the soft SUSY-breaking 
sector of the MSSM (i.e. Eq. 5.2). These parameters are divided into four 
categories: the gaugino masses Mi, M2 and M3, the trilinear couplings 
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au, ad and ag, the squark and slepton squared- masses hiq, 111^, m|-, 
and m|, and the Higgs squared-mass parameters fn^^, ^'h^ ^-^id b. 

For the gaugino masses Mi {i = 1,2,3), the one-loop RGEs read (in 
terms of the quantities bi and Oi) 

^Mi = -^aiM, =^ ^(Mi/ai) = 0. (5.27) 
at zvr at 

This shows that the three ratios Mi/ui are scale independent, i.e. do 
not change with renormahsation scale. In addition, we observed that the 
gauge couplings unify at the GUT scale Mqut- This suggests to assume 
that the gaugino masses also unify at Mqut • This is a popular assumption 
and the unified value is usually denoted as mi/2 (see also section 5.3). 

For the trilinear couplings a^, ad and ag, it is common to assume that 
they are proportional to the Yukawa couplings yu, yd) ye; and accordingly 
have the diagonal forms au = diag{0,0,at), ad = diag{0,0,ab) and = 
diag{0,0, ar)- In this case, the one-loop RGEs for the parameters at, a^ 
and ttr have the following form [56]: 
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+6abylyr + yA692M2 + -9lM^) ■ (5-28) 

5 J 

The squark and slepton squared-masses, that are 3x3 matrices, are 
usually assumed to have diagonal forms (in order for the potentially dan- 
gerous flavour-changing and CP-violating effects in the MSSM to be sup- 
pressed). For example for hiq this means that mg = diag{mQ^ , mg^ , rriQ^) 
Analogous forms are assumed for the other squarks and sleptons. 

The one-loop RGE expressions for these squared-masses are relatively 
lengthy and we therefore give here, as an example, the expressions only 
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for the third-family quantities [56]: 
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= 2|yt|2(m2,^+m^3+mi3) + 2|at|2, 
Xfe = 2|?/6|2(m|^ +m^3+m^ J + 2|ab|2, 



(5.30) 



Finally, let us look at the RGEs for the Higgs squared-mass parameters 
, '^Ui ^- They are of the following forms [56] : 

3Xt - ^gl\M2\^ - \g\\Mr\^ + '^-g\S 

3Xb + X, - e^llMsl' - Iffi'jMil^ - \g\S 



K^ylVt + Sy^yfe + yXvr - 3^2 - ^5i) 

5 

M(6aty; + 6afey,* + 2a,y; + e^^Ma + \g\M{)\ , (5.31) 
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with Xt, Xfe, Xt and S defined in Eqs. 5.30. 



5.1.7 Parameter space 

Looking at the full Lagrangian of the MSSM, including soft terms (Eq. 5.2), 
one realises that the model, compared to the SM, possesses a much larger 
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number of free parameters, most of which come from the SUSY-breaking 
sector, i.e. soft terms. Indeed only one of the new parameters, fi, belongs 
to the SUSY invariant sector. 

Let us try to count the parameters in the MSSM: In the fermion sector, 
the theory has 5 Hermitian 3x3 mass-squared matrices plus 9 complex 
3x3 trilinear coupling matrices. These give 5x3x3 + 6x3x3x2 = 153 
real parameters. Due to some field redefinitions, this number is reduced to 
110. The gauge sector of the MSSM has 3 usual gauge couplings g' , g and 
gs plus the QCD vacuum angle 6 (see section 3.6 and references therein). 
If we add to these the 6 gaugino masses, this gives 3 + 1 + 6 = 10. Using 
a field transformation, one of the CP-violating masses can be removed, 
reducing the number by 1 and giving rise to 9 parameters in the gauge 
sector. The Higgs terms of the soft SUSY-breaking sector contribute by 2 
real squared masses mj^^ and rn^^, and 1 complex coefficient b, and this 
gives 1 + 1 + 2 = 4 parameters. The SUSY invariant sector contributes 
with the parameter fj, (being complex). A field definition helps reduce 1 
parameter and we are left with 1. We can now count the total number of 
free parameters in the MSSM: It is 110 + 9 + 4 + 1 = 124. The number 
of free parameters in the SM is 19, meaning that the full MSSM has 105 
more parameters than the SM. 



5.2 ... and beyond 

We introduced in section 5.1 the MSSM as a phenomenological super- 
symmetric extension of the SM that contains the minimum number of 
new fields needed for the supersymmetrisation process. Respecting the 
renormalisability, gauge symmetry and B L conservation in the SM, 
the MSSM allows for the most general terms in the soft SUSY-breaking 
sector of the theory. Although there are models with less number of free 
parameters (coming either from particular SUSY-breaking mechanisms or 
from purely phenomenological assumptions, as we will see in section 5.3), 
most of them have usually the same field content as the MSSM while 
defined over particular parts of the MSSM parameter space. The MSSM 
is therefore considered by almost all particle physicists as the most in- 
teresting framework to study supersymmetric extensions of the SM and 
provides extensive scope for various phenomenological investigations of 
such extensions. On the other hand, as we will see in section 5.3, the 
MSSM parameter space is already too large and cannot be fully analysed. 
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Although these all indicate that the main focus of SUSY phenomenolo- 
gists will, for the next few years, be the MSSM and its subclasses, one 
should however notice that attempts for going even beyond the MSSM 
have already begun. This is mainly because the MSSM itself is not flaw- 
less and in fact possesses several problems. In what follows, we give two 
examples of SUSY models beyond the MSSM (and the problems they aim 
to solve) that have attracted much attention in the past few years. 

5.2.1 BMSSM 

In the MSSM, the quartic Higgs coupling constant is given entirely in 
terms of the electroweak gauge couplings g and g'. This can be seen 
e.g. from Eq. 5.3 where this quartic coupling (that we denote by A) is 
expressed as A = {g'^ + g''^)/8. As a result, the value of A becomes small 
(A 0.07). This smallness poses a problem. 

As we discussed in section 5.1.3, the theoretical tree- level lightest 
Higgs mass nih is too low to be consistent with experimental constraints. 
One can show that the value of this tree-level mass is related to the value 
of the quartic Higgs coupling constant A and the smallness of the latter 
implies the same for the former. This all means that in the MSSM, large 
loop corrections are needed to bring the Higgs mass above the experimen- 
tal bounds. The experimental lower bound on the Higgs mass (i.e. 114 
GeV) excludes only a small part of the SM parameter space. On the 
contrary, the requirement of fairly large loop corrections excludes most 
of the natural values for the MSSM parameters. This is called the 'little 
hierarchy problem' of the MSSM (see e.g. ref. [100]). 

Many possible extensions of the MSSM have been proposed so far 
in attempt to address the little hierarchy problem (for a review, see 
e.g. ref. [101] and references in ref. [102]). One example is called 'be- 
yond the MSSM' (or BMSSM) [103] (see also refs. [104, 105, 106, 107] 
and references therein for more phenomenological studies of the model). 

In this model, with an effective field theory approach (see e.g. ref. [57]), 
the field content of the MSSM remains intact, while the quartic Higgs 
couplings are modified by the effects of some new physics that might 
exist at high energies. In the MSSM, the quartic Higgs terms are exactly 
supersymmetric, because in the MSSM the only SUSY-breaking terms 
that are allowed are those with mass dimension three or lower (soft terms) . 
As we discussed in the previous section, this condition is imposed so as to 
prevent the reintroduction of quadratic divergences. In an effective field 
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theory approach however, the effects of new physics can be parametrised 
by non-renormahsable terms. 

In the simplest version of BMSSM only two leading terms with partic- 
ular properties are added to the MSSM. One is a dimension-five operator 
W5 that is added to the MSSM superpotential Wmssm- has the 
following form: 

W5 = -^{HuHdf. (5.32) 
A* 

Here ei is a new effective real parameter of the model that is free and 
should be determined experimentally. 

The other new operator of the theory is another dimension-five oper- 
ator that is added to the soft SUSY-breaking sector of the MSSM. This 
term is 

e2{HuHdf + h.c. (5.33) 

with 62 being another new effective real parameter. 

After introducing the new parameters e\ and €2, we get the following 
new terms added to the full MSSM Lagrangian: 

5C = -2ei{HM{HlHu + H^Hd) - e2{HuHd? + p [2{HuHd){HuHd) 

+2{HuHa){HuHa) + {HuHay + (^^^d)^] + h.c. (5.34) 

The above correction terms affect the Higgs potential, neutralino masses 
and mixings, and chargino masses and mixings [104]. The first effect is 
helpful in alleviating the little hierarchy problem of the MSSM since the 
Higgs quartic self-coupling now receives new contributions with param- 
eters that are not fixed by the theory. BMSSM also has interesting im- 
plications for cosmology, in particular dark matter phenomenology and 
baryogenesis that are studied for example in refs. [105, 104]. 



5.2.2 NMSSM 

One other issue that cannot be explained in the MSSM is the so-called 
//-problem [108], which has to do with the /i-term in the MSSM super- 
potential (Eq. 5.1). The n parameter is dimensionful (i.e. with positive 
mass dimension) but its value can be arbitrary and is not associated with 
any particular scale of the theory such as the SUSY-breaking scale. Phe- 
nomenologically, fi is required to have a value close to the EWSB scale, 
and this scale is not 'natural'. 
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The //-problem has served as the primary motivation for proposing an 
extension of the MSSM that is called 'Next-to- Minimal' Supersymmetric 
Standard Model (or NMSSM) (for a review, see e.g. refs. [109, 110]). 

In the NMSSM, a new gauge singlet is added to the MSSM and conse- 
quently an effective //-term is generated dynamically (i.e. spontaneously). 
The introduction of the additional gauge-singlet superfield is essentially 
the price to pay in order to solve the //-problem. The NMSSM however 
provides explanation for the other problems of the MSSM, including the 
little hierarchy problem described in the previous section. This is because 
the Higgs-boson sector of the theory is much less restricted compared to 
the MSSM and the predicted lower mass bound on the Higgs mass is 
in general substantially shifted. The NMSSM contains two additional 
Higgs bosons and one additional neutralino (that is called singlino). The 
model generally offers different Higgs-boson phenomenology, compared to 
what we expect from the MSSM, with interesting implications for collider 
searches. 

5.3 Shrinking the parameter space 

We saw in section 5.1.7 that the full MSSM Lagrangian possesses 124 
free parameters that should be determined experimentally. This huge 
parameter space makes the phenomenological studies of the model prac- 
tically difficult. In addition, the structure of the theory is such that most 
combinations of the parameters give experimental predictions that are ex- 
cluded. For example many non-diagonal terms in the MSSM Lagrangian 
generate too large FCNCs, at levels that are experimentally excluded. 
The parameter space is therefore highly porous and has a very non-trivial 
structure. It is therefore rather crucial to work with sub-models of the 
MSSM with substantially less numbers of free parameters. In this section, 
we briefly review some approaches to this problem and different strategies 
in reducing the size of the parameter space. 

5.3.1 Connections with fundamentals 

As we discussed earlier in this chapter, all the new parameters in the 
MSSM, except one, come from the SUSY-breaking sector, i.e. the soft 
SUSY-breaking terms of Eq. 5.2. We argued that the soft terms are 
introduced based on some theoretically and/or phenomenologically moti- 
vated properties that one expects the MSSM to possess at low energies. 
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Figure 5.2. Schematic description of supersymmetry breaking mecha- 
nisms. 



Although such terms break supersymmetry explicitly, they are widely con- 
sidered as useful low-energy terms that parametrise our ignorance of some 
underlying mechanism that breaks SUSY spontaneously (see section 4.5). 

Various SUSY- breaking mechanisms have so far been proposed, each 
of which imposes its own set of relations between different soft terms 
and corresponding parameters. For example in many SUSY-breaking sce- 
narios some generally non-zero or complex parameters of the MSSM are 
predicted to be vanishing (or extremely small) or real, and universality 
conditions are imposed on some otherwise unrelated parameters. Usually 
in these scenarios, the parameter space is significantly contracted and the 
analysis of the model predictions becomes considerably easier. 

In section 4.5, we discussed general strategies one can use for breaking 
supersymmetry spontaneously. They were classified into two categories 
of F or D-term SUSY-breaking mechanisms. 

Unfortunately, the structure of the MSSM does not allow any of the 
two above strategies for breaking supersymmetry to be realised without 
extending the field content of the theory [56]: (1) Giving a VEV to the 
D-term associated with the U{1)y part of the MSSM gauge structure 
has turned out to give rise to an inappropriate mass spectrum. (2) The 
MSSM does not contain any gauge singlet with a corresponding F-term 
that develops a VEV. 

Even if we extend the MSSM in such a way that it includes new 
supermultiplets (including gauge singlets necessary for i^'-term SUSY- 
breaking) , it is highly difficult to achieve phenomenologically viable spon- 
taneous supersymmetry breaking only using renormalisable interactions 
at tree level. For example, due to the existence of particular 'sum rules' 
for the tree-level mass terms, some squarks and sleptons in the MSSM 
have to have experimentally excluded low masses. 

These types of problems can be evaded if one assumes that supersym- 
metry is broken in a different sector which communicates with the MSSM 



5.3. Shrinking the parameter space 



65 



particles only indirectly i.e. via either 'non-renormalisable' interactions or 
couplings at 'loop level'. The new SUSY-breaking sector is accordingly 
referred to as "the hidden sector" as opposed to the MSSM sector being 
the visible sector. The MSSM soft terms are then thought of as being the 
results of the 'mediation' of SUSY-breaking from the hidden sector to the 
visible sector (see Fig. 5.2 for a schematic picture of viable SUSY-breaking 
mechanisms). 

In order to understand a SUSY-breaking mechanism in full details, one 
needs to know exactly which theory governs the laws of physics at high 
energies. Such an understanding is still lacking and the usual approach 
is therefore to propose models of SUSY-breaking based on rather general 
frameworks and assumptions that capture interesting features of more 
fundamental theories. 

Existing SUSY-breaking models can be categorised into two general 
classes: either they are based on local supersymmetry (or supergravity; 
see sections 3.5 and 4.4) in which SUSY breaking is mediated through 
gravitational non-renormalisable interactions, or the soft SUSY-breaking 
terms are generated only at loop level [56]. As we stated earlier, in both 
classes the problematic sum rules are circumvented and supersymmetry 
is broken spontaneously. 

In the following sections we briefly review three main SUSY-breaking 
scenarios that have been proposed in the literature and have as of yet 
received most of the attention. We describe their theoretical founda- 
tions only in a few words and pay most of our attention to the simplest 
(but practically the most interesting) models within each framework and 
the corresponding phenomenological aspects (for detailed introduction to 
SUSY-breaking mechanisms, see e.g. refs. [54, 56]). 

Planck-scale-mediated SUSY breaking: mSUGRA 

In this class of models (which has been historically the most popular 
one) supersymmetry is broken in a hidden sector which interacts with the 
visible sector only through gravitational effects near the Planck scale, 
hence the names "gravity-mediated" SUSY-breaking or "Planck-scale- 
mediated" SUSY-breaking (PMSB) mechanism [111, 112, 113, 114, 115, 
116, 117]. 

The presence of gravitational interactions in these models requires a 
supersymmetric theory that includes gravity. We mentioned earlier in 
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sections 3.5 and 4.4 that such a theory can be provided if the global su- 
persymmetry is upgraded to a local one, i.e. the parameters of supersym- 
metry transformations become space and time dependent. Supersymme- 
try is a spacetime symmetry and its local version automatically involves 
gravitation and is therefore also referred to as supergravity or SUGRA. 
The phenomenologically interesting versions of SUGRA, like any other 
known quantum field theory that contains gravity has turned out to be 
non-renormalisable. 

The non-renormalizable terms in the SUGRA effective Lagrangian 
are fortunately suppressed by powers of 1/Mp and their effects are conse- 
quently negligible for all phenomomenological studies that normally probe 
energies below TeV scales. Such terms can however provide a mechanism 
through which the hidden and visible sectors communicate and SUSY 
breaking is mediated from the former to the latter. 

The supergravity Lagrangian, amongst other terms, usually contains 
non-renormalisable (NR) terms that look like 

Cnr = -^FJ^{X'^,ct>,)-^FF*g{cp'). (5.35) 

Here, (pi and Xa are the scalar and gaugino fields of the MSSM and T and 
Q are two functions of the fields whose exact forms are not important for 
our discussion here (see e.g. ref. [56] for details). F denotes the auxiliary 
field corresponding to a chiral multiplet that is assumed to exist in the 
hidden sector. By choosing the right functional forms for and Q, one 
can show that Eq. 5.35 is equivalent to Eq. 5.2 for the MSSM soft terms 
with a mass scale of 

m-soft ~ (5.36) 

The SUSY-breaking scale in the hidden sector M^jj^ is determined 
by the VEV {F) such that 

Msws^~y(F). (5.37) 

This means that in order to have mgoft ~ 1 TeV (as we require for weak- 
scale supersymmetry), il^ygy should be about 10^^ GeV. 

Functions J- and Q in Eq. 5.35 in general contain many free parameters 
whose values are to be determined by the exact form of the underlying 
supergravity theory. One popular way of simplifying the PMSB scenario 
is to work in the framework of "minimal supergravity" or mSUGRA (also 
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called the 'supergravity-inspired scenario'). In this case all the soft pa- 
rameters are fully determined by only four parameters 

mi/2,mo,vlo,-B, (5.38) 

plus the yu-parameter of the MSSM superpotential. The soft parameters 
of Eq. 5.2 in this case read 

Ms = Ma = Ml = mi/2, 

2 2 2 2 2 2-1 

mq = niu = nid = "^L = ™e = "^o 1> 

2 2 2 

T^H^ = T^Hi = mo, 

au = ^oyu, ad = AoYd, ae = ^oYe, 

b = Bfi. (5.39) 

It is important to notice that the above unification relations are de- 
fined at the natural scale for gravity i.e. Mp. One then needs to use 
renormalisation group equations (as we discussed in section 5.1.6) and 
evolve the soft parameters down to the electroweak scale so as to obtain 
the interesting low-energy quantities such as the MSSM mass spectrum. 
From a practical point of view however, it has become common to assume 
that the relations (5.39) are valid at the grand unification scale Mqut aud 
then start the RGEs from that scale instead of Mp. This is mostly be- 
cause we do not know much about how the RGEs behave at scales between 
Mqtjt and Mp whereas the observed unification pattern for the MSSM 
gauge couplings (see section 5.1.6) indicates that our understanding of 
the RGEs at scales below Mqut may not be far from reality. It has been 
shown that the effects neglected by using Mqut instead of Mp in the 
RGEs, are probably not significant. They may however cause additional 
important effects [118, 119, 120, 121]. 

The much simpler parameter space of mSUGRA (with only 5 new pa- 
rameters compared to the SM) has made it one of the most popular models 
for SUSY breaking and after some small modifications (as we will see in 
section 5.3.2) the model has received the bulk of attention in phenomeno- 
logical studies of supersymmetry. According to some authors [122], the 
parameter space of mSUGRA is even simpler than what we described 
here. This is achieved by imposing the extra relation Aq = B + niQ which 
eliminates one more free parameter from the MSSM parameter space; this 
parameter is usually taken to be i3. 
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Gauge-mediated SUSY breaking: mGMSB 

In gauge-mediated SUSY breaking (GMSB) [123, 124, 125, 126], in con- 
trast to the PMSB scenario, the effects of SUSY breaking are not medi- 
ated from the hidden sector to the observable sector (i.e. the MSSM) using 
non-renormahsable interactions. GMSB models circumvent the problems 
with the sum rules in a different way, and that is to generate the soft 
terms through 'radiative' interactions, i.e. couplings at 'loop level' rather 
than tree level. 

The interactions responsible for generating soft terms in a GMSB 
setup are purely gauge, similar to the gauge interactions in the SM or the 
MSSM. The idea is the following: Assume that there are some new chiral 
supermultiplets that, on the one hand, couple to the source of supersym- 
metry breaking in the hidden sector, and on the other hand, couple to 
the MSSM particles through the SM gauge bosons and their superpart- 
ners gauginos. These new chiral supermultiplets are called "messengers" . 
Now, when SUSY is broken in the hidden sector, the effects are first con- 
veyed to the so-called 'messenger sector' and then communicated to the 
visible sector through gauge interactions and 'only radiatively' (i.e. only 
at loop level). 

One should notice that even in the GMSB scenarios gravitational in- 
teractions exist between the fields in the hidden, messenger and observ- 
able sectors. However, such gravitational effects are much weaker than 
the gauge effects. In other words, the gravitational effects are dominated 
by the gauge effects and consequently do not play any important role in 
breaking supersymmetry. 

In the simplest GMSB model, called 'minimal GMSB' (or mGMSB), 
the messenger fields couple to a gauge-singlet chiral supermultiplet. The 
scalar component of the supermultiplet (denoted by S) and its corre- 
sponding auxiliary field Fs develop VEVs (S) and (Fs), respectively, and 
therefore break supersymmetry. The scale of SUSY breaking in the mes- 
senger sector is ~ \/~{Ps) and this should be distinguished from 
the SUSY-breaking scale in the hidden sector M^^^^^- that is associated 
with its own auxiliary field F. 

The mGMSB contains 6 free parameters [54] 

A,M,n5,tan^,sgn;U,Cgrav (5.40) 

Here, the two parameters A (that sets the mass scale of the MSSM spar- 
ticles) and Cgrav (which is called gravitino mass parameter) are defined 
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where A is the common messenger-sector Yukawa couphng. The parame- 
ters 77-5 and M (M > A) set the number of messenger multiplets and the 
mass-scale associated with the messenger sector, respectively. All the soft 
SUSY-breaking terms of the MSSM are obtained in mGMSB by means of 
the RGEs that are evolved from the scale of M down to the electroweak 
scale, tan /3 is the ratio of up-type to down-type Higgs VEVs at the elec- 
troweak scale (see section 5.1.3). can in principle be treated as a free 
parameter like in mSUGRA, but its magnitude is usually fixed in this 
model by imposing 'radiative electroweak symmetry breaking' (REWSB) 
condition at the weak scale, while its sign is to be determined experi- 
mentally (see section 5.3.2 for more details about REWSB and a similar 
assumption for a different model). An interesting feature of mGMSB 
is that the trilinear coupling parameters of 5.2, i.e. au, and ae only 
arise at two-loop level and are therefore very small; they are frequently 
assumed to be vanishing [56] . 

Extra-dimensional-mediated SUSY-breaking: mAMSB 

In this class of SUSY-breaking models, the central idea is that the observ- 
able and hidden sectors correspond to two different spacetime manifolds 
that are hovering in a bulk and separated physically (for detailed discus- 
sion, see e.g. refs. [54, 56]). In the simplest models of extra-dimensional- 
mediated supersymmetry breaking (XMSB), it is assumed that each of 
the two sectors of the theory is confined to a 4-dimensional brane with 
a 5-dimensional bulk spacetime between them. The hidden and observ- 
able sectors then communicate in a manner that depends on whether the 
MSSM gauge supermultiplets are allowed to propagate in the bulk or not. 
If so, they can mediate supersymmetry breaking. It has been shown that 
in the simplest version of such models soft SUSY-breaking is dominated 
by gaugino masses; this scenario is therefore usually called "gaugino me- 
diation" (see e.g. refs. [127, 128, 129]). 

One other possibility is that the gauge supermultiples are not allowed 
to freely travel in the bulk and, analogously to the chiral multiplets, are 



70 



Chapter 5. Supersymmetry in real life 



confined to the visible brane (i.e. the MSSM sector). This means that 
like in the case of the gravity-mediated SUSY-breaking, some supergrav- 
ity effects should be responsible for the transmission of SUSY-breaking 
from the hidden to the visible sector. In an interesting class of these 
models, which is called "anomaly- mediated supersymmetry breaking" (or 
AMSB) [130, 131], the MSSM soft terms are generated at loop level due to 
an anomalous violation of a particular symmetry called 'local supercon- 
formal invariance'. Again, as in the GMSB scenario, gravity-mediation 
is present here but its effects are dominated by the anomaly-mediation 
interactions. 

The original AMSB scenario has many unique properties that make 
the scenario quite interesting. For example all the soft terms generated 
from this model can be written in terms of only one free parameter which 
is the gravitino mass The model is however not viable: it can 

be shown that the sleptons have negative squared-masses (i.e. the par- 
ticles are tachyonic). There have been various proposals for modifying 
the theory to circumvent the tachyonic mass problem (see e.g. ref. [56] 
and references therein) . One phenomenologically motivated approach has 
been to add a new parameter mo (usually set at the GUT scale) to the 
model which provides large contributions to the slepton squared-masses 
and makes them positive. This so-called 'minimal AMSB' (or mAMSB) 
is then characterised by four free parameters [54] 

"^o,"^3/2,tan/3,sgn/i. (5.42) 
5.3.2 Phenomenological assumptions 

An entirely orthogonal approach to handle the large MSSM parameter 
space, is to simply impose phenomenologically justified assumptions and 
simplifications to the parameters without relying on any particular un- 
derlying supersymmetry-breaking mechanism. We briefly review some of 
these simplified models in this section. 

Low-energy models: MSSM-7,8,10,11,18,19,24 

It is already known that extensive regions of the MSSM parameter space 
are excluded experimentally. This is mainly due to the fact that many of 
the soft terms in the SUSY-breaking sector of the MSSM can introduce 
large flavour-changing neutral currents (FCNCs) or new sources of CP- 
violation that are strongly constrained by experiments. These dangerous 
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terms should therefore be suppressed. Examples of soft parameters that 
typically produce large FCNCs include off-diagonal entries in the trilinear 
coupling matrices a^, ad and ag and sfermion mass matrices m^, m^, 
m^, and m|. It is common in the SUSY-phenomenology community 
to accordingly approximate many of these off-diagonal parameters to zero 
(see e.g. ref. [132]), although this is certainly not theoretically motivated 
at this level. 

Let us emphasise here that many of the SUSY-breaking scenarios we 
discussed in the previous section, naturally lead to the suppression of dan- 
gerous off-diagonal terms; an example is mSUGRA. Inspired by mSUGRA 
assumptions, in constructing some phenomenological sub- models of the 
MSSM, one can assume additional properties for the mass and trilinear 
coupling matrices such as 'reality' and 'universality'. This contraction of 
the parameter space has turned out to be so helpful that in some popular 
cases the number of new free parameters of the MSSM has seen dramatic 
reductions, e.g. from 105 to 7 (see section 5.1.7). 

The phenomenological version of the MSSM with the largest number 
of free parameters that has been analysed phenomenologically to date, 
has 24 non-SM parameters (compared to 105 for the full MSSM) and is 
called MSSM-24 or simply 'phenomenological MSSM' (see e.g. refs. [133, 
134]). In one version of MSSM-24 the non-SM parameters are 3 trilinear 
couplings at, and ar for the third-family, 5 x 3 = 15 diagonal entries for 
sfermion squared- mass matrices itiq , m^, m^, 01^ and m? corresponding 
to the three different families (no universality is imposed), 3 gaugino 
masses Mi, M2 and M3, the ratio of up- type to down- type Higgs VEVs 
tan/3, the mass of the pseudoscalar Higgs tua and the parameter of the 
superpotential. All these parameters are defined at the electroweak scale. 

By imposing more and more simplifications, one can decrease the num- 
ber of free parameters and this has been done in the literature in various 
ways. The ones that have been widely used so far have 19, 18, 11, 10, 8 
or 7 parameters. 

As one example of the low-dimensional models, let us briefly describe 
MSSM-7. In this model only the trilinear parameters at and a;, of MSSM- 
24 are taken to be generally non-zero [a-r is set to zero), all 15 diagonal 
sfermion squared-mass parameters are assumed to be equal (with the 
universal mass parameter mo), and gaugino mass parameters Mi, M2 
and M3 are related in the following way: 

Ms = — sin^ ewM2 = 1— cos^ OwMi. (5.43) 

a b a 
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Figure 5.3. Renormalisation group evolution of scalar and gaugino mass 
parameters mi/2 and mo in tlie constrained MSSM. Tiie quantity + 
rrijj^ runs negative and breaks electroweak symmetry. Squark and slepton 
squared-mass parameters remain positive. Adapted from ref. [56]. 

The above relations hold at the electroweak scale and are inspired by 
mSUGRA model described in section 5.3.1. These relations are the results 
of the RGEs evolved from the GUT scale down to the electroweak scale 
(for mSUGRA). 

So far we have 4 parameters at, ab, rriQ and, say, M2. If we add /i, rriA 
and tan (3 to these, we obtain the full set of free parameters for MSSM-7. 

High-energy models: CMSSM and NUHM 

All phenomenological assumptions we made in the previous subsection 
were imposed on the MSSM parameters at 'low energies', i.e. directly 
at the electroweak scale. They are therefore different from the relations 
imposed by certain SUSY-breaking scenarios of section 5.3.1 in that the 
latter are usually applied at 'high energies', in particular the GUT scale. 
The hybrid approach is then to assume purely phenomenologically in- 
teresting simplifications and relations between the model parameters at 
high energy scales. The so-called "constrained MSSM" (or CMSSM) is 
an example [135]. 

In the CMSSM, inspired by mSUGRA model, various universality 
assumptions are imposed on gaugino and scalar mass parameters, as well 
as trilinear a-term couplings. One central assumption in the CMSSM 
that makes it different from mSUGRA is that the condition of 'radiative 



5.3. Shrinking the parameter space 



73 



electroweak symmetry breaking' (REWSB) is required to be fulfilled when 
RGEs evolve the parameters from the GUT to the electroweak scale. 

In section 5.1.3, we described the Higgs sector of the MSSM, as well 
as the conditions that must be fulfilled in order for the Higgs mechanism 
to break electroweak symmetry spontaneously. In fact, EWSB can be 
achieved in the MSSM if any of the Higgs squared-mass parameters Tn^^^ 
and m^H^ in the soft SUSY-breaking sector (see Eq. 5.2) becomes negative. 
This can be certainly done in the MSSM by hand, namely by setting the 
parameters directly at the electroweak scale (of course with the condition 
that the relations 5.5 must be satisfied). The process so far has been very 
similar to the SM case, with the difference that the latter has only one 
squared-mass parameter to be tuned. 

The nice feature of the MSSM that makes it remarkably interest- 
ing from the EWSB point of view, is the observation that even if one 
gives positive values to the Higgs squared-masses m^j^ and at some 
high energy, the renormalisation group equations can make negative 
at the weak scale. This includes models with universal mass assump- 
tions, such as mSUGRA where all scalar mass parameters are unified at 
the GUT scale (see Fig. 5.3). In this latter case, rn^^^ which has the 
same value as the other scalar squared-mass parameters of Eq. 5.2 (see 
Eqs. 5.39), runs negative and makes the quantity + rn^u^ also nega- 
tive. The latter condition then breaks electroweak symmetry. It can be 
observed from Fig. 5.3 that the squark and slepton squared-masses remain 
positive at all scales. This mechanism through which rn^^j^ turns negative 
via the renormalisation group evolution, is called radiative electroweak 
symmetry breaking [54]. It has been shown that REWSB can occur over 
a wide range of MSSM parameters if the top quark mass mt has a value 
between about 100 and 200 GeV, which is obviously the case {mt — 172 
GeV) [54]. This observation has provided one of the strong motivations 
for supersymmetrising the SM (see also section 3.2). 

The first equation in Eqs. 5.9 of section 5.1.3 implies that if elec- 
troweak symmetry is broken appropriately, the value of the mSUGRA 
parameter B can be fully determined in terms of ^, the quantity tan /3 
and the mass parameters of the theory that their weak-scale values are 
calculated by the RGEs (remember that h = Bfi). In addition, the second 
equation fixes the magnitude of the parameter in terms of tan /3 and 
the masses only leaving its sign to be determined experimentally. 

In the CMSSM, assuming that Eqs. 5.9 hold, the high-scale parameter 
B is eliminated in favour of tan/3. This leads to the following set of free 
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parameters for the model (fom' continuous and one discrete): 

mo,mi/2,-4o,tan/3,sgn/i. (5.44) 

One concrete prediction of the CMSSM, as well as mSUGRA, is that 
(as we have already indicated in Eq. 5.43) the three gaugino mass pa- 
rameters Ml, M2 and M3 are in a particular ratio. This is often shown 
approximately as Mi : M2 : M3 ~ 1 : 2 : 7 (according to Eq. 5.43). 

We mentioned earlier (when discussing the RGEs in section 5.1.6) that 
one reason why the gaugino masses are assumed to unify at the GUT 
scale is the intriguing fact that their ratios to the corresponding gauge 
couplings (which are unified at the GUT scale) do not evolve with energy. 
Such a motivation does not exist for the scalar masses and the unification 
assumption on their GUT values seems to be a rather arbitrary choice. 
This means that one can make other equally reasonable choices. 

One popular example is the so-called "non-universal Higgs model 
(NUHM)" (see e.g. ref. [136]). In this slightly less restrictive framework, 
only squak and slepton masses are given a universal mass mo at the GUT 
scale whereas mn^ and m//^ are treated as being independent parameters. 
This seems to be a reasonable relaxation of the universality condition be- 
cause even in grand unified theories the Higgs scalars do not necessarily 
belong to the same multiplet as the sfermions and there is therefore no 
reason for treating them on the same footing. The relaxation made in 
the NUHM model, introduces two new parameters mu^ and tuHci which 
are commonly traded for the more phenomenologically interesting weak- 
scale quantities and m^. This gives the model a total number of 7 free 
parameters.^ 



^Strictly speaking, the model we described here is the so-called NUHM2 model. The 
NUHMl model is slightly different e.g. in that it has only 1 more free parameter than 
the CMSSM. 



Chapter 6 



Observational constraints on 
supersymmetry 

In the previous chapters, we gave a review of supersymmetry in general, 
as well as various supersymmetric models that have been put forward 
as solutions to different problems in particle physics and cosmology. Su- 
persymmetry, as any other theories, has to be tested observationally. In 
fact, questions such as "whether supersymmetry is a correct description 
of Nature at high energies" and if so "which supersymmetric model pro- 
vides the best such description" all need comparison of SUSY predictions 
with real experimental data. In addition, we noticed that even within 
the framework of each supersymmetric model, one usually encounters a 
large number of free parameters that need to be determined experimen- 
tally. One of the main objectives of the present thesis has been to provide 
powerful tools and techniques for SUSY parameter estimation when its 
predictions face different observational data. We therefore in this chap- 
ter briefly review some of the most important observational strategies in 
testing supersymmetric models and constraining their parameter spaces. 

6.1 Supersymmetric WIMPs 

We discussed in section 2.3 that weakly interacting massive particles 
(WIMPs) provide an elegant solution to the dark matter problem. We 
also mentioned that supersymmetric theories contain viable dark matter 
candidates such as the lightest neutralino (that we will simply call "the 
neutralino" from now on), gravitino and axino, of which the neutralino 
has been the most popular WIMP dark matter candidate to date (see also 
section 5.1.4). 
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A copious number of experiments are now looking for WIMPs and 
various observational constraints on their properties can naturally provide 
potential tests of supersymmetric models that contain WIMPs. In this 
section, we give an overview of the main strategies and techniques used to 
search for WIMPs, as well as the major existing observational constraints 
on their properties. Since the neutralino has been the only dark matter 
candidate we have studied in all analyses done in this thesis, we therefore 
restrict our discussions to its properties whenever we speak of a particular 
type of WIMPs. 

6.1.1 Cosmological relic density 

One of the accurately measured observables that significantly constrain 
the parameter spaces of most SUSY models, is the present amount of dark 
matter in the Universe, i.e. VLdm = — ^6 introduced in section 2. We 
discussed in section 2.3 that if the thermally-produced- WIMP scenario is 
correct, the present dark matter particles should have been created at the 
freeze-out period in the early Universe when the expansion rate of the 
Universe hits the WIMP interaction rate and WIMPs fall out of chemical 
equilibrium (and decouple). We also gave a brief exposition of some of the 
basic principles in calculating the dark matter relic density from WIMPs. 

The value for the relic density we presented in Eq. 2.8 was however 
only approximate and one needs to solve the Boltzmann equation 2.6 so 
as to obtain the exact value for a particular WIMP candidate. The relic 
density of a WIMP in general depends on both its mass and annihilation 
cross-section which in turn contain information from the particle nature of 
the WIMP. These provide a connection between the actually observed relic 
density and the theoretical value predicted for example by a particular set 
of model parameters. Fig. 6.1 shows an example analysis of the CMSSM 
(see section 5.3.2) where all but two parameters mo and mii2 are fixed 
and only the relic density constraint is imposed [140] . It is interesting to 
see that this observable provides such a tight constraint on the CMSSM 
parameter space (cyan and dark blue strips) even if no other constraints 
are used. 

It is important to note here that the processes that give rise to the 
final value of the relic density for WIMPs are in most cases much more 
sophisticated than the simple case we have discussed so far. This means 
that for example the Boltzmann equation 2.6 has to be modified. One 
example is that the relic density of dark matter depends upon the history 
of the expansion rate in the early Universe before and during freeze-out. 
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Figure 6.1. The mi/2-mo planes for the CMSSM with tan/3 = 35 and 
assuming fi < 0. The region with cyan shading shows the allowed values 
of mi/2 and mo when the cosmological constraint 0.1 < Q^h^ < 0.3 has 
been imposed, and the region in cyan displays the allowed values when 
a tighter constraint (0.094 < fl-^h^ < 0.129) has been imposed. Adapted 
from ref. [140]. 

Strong modifications of that epoch can substantiaUy impact the value of 
the rehc density today. One other important example is when at the 
time of freeze-out some other particles exist that have masses close to the 
WIMPs and also share a quantum number with them. These particles 
interact with the WIMPs and in many cases can enhance the annihila- 
tion process and therefore change the relic density dramatically. These 
effects are usually referred to as 'coannihilations' [141]. The effects of 
including coannihilations in the calculations of relic density has been ex- 
tensively studied in the literature (see e.g. refs. [142, 143]). The presence 
of coannihilations make the relic density calculations quite difficult and 
they have to be done numerically. For example, for the case of the su- 
persymmetric neutralino WIMPs the relic density is usually calculated 
using advanced computer packages that take various effects into account, 
including coannihilations. The most popular publicly available numeri- 
cal codes that calculate relic densities (amongst other things) are Dark- 
SUSY [132] (available from ref. [137]) and micrOMEGAs [138] (available 
from ref. [139]). 

Finally let us mention here that our best estimation of the value of 
the dark matter relic density at the moment comes from the observations 
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Figure 6.2. A WIMP scattering off normal nuclei in a dark matter direct 
detection experiment. 



of the cosmic microwave background (CMB) by the WMAP satelHte [4]. 
The value is ^dm = 0.222 it 0.026. We have used a similar value in our 
analyses of this thesis. A tighter constraint is expected to be provided by 
e.g. the Planck satellite [144]. 



6.1.2 Direct detection 

One important property of WIMPs is that they interact with the SM par- 
ticles, although very weakly. This property can be used to directly search 
for WIMPs, i.e. by looking for any interactions between them (that are 
supposed to fill our galaxy if they form dark matter) and some normal 
matter particles on Earth [145, 31]. This simple idea, as one of the most 
promising search strategies, has stimulated many experimental groups to 
build different small and large scale detectors looking for WIMP-SM in- 
teraction signals. Such interactions are sought for in 'direct detection' 
experiments by recoding nuclear recoils when WIMPs scatter off the de- 
tector nuclei (see e.g. Fig. 6.2). This is usually done in the detectors in 
different ways that are usually various combinations of three different de- 
tection techniques: 'ionisation', in which the atoms of the target material 
becomes ionised by the transferred 'recoil energy'; 'scintillation', in which 
particular materials known as scintillators are used to measure the fiuo- 
rescent radiation produced by electrons in the target material when they 
decay after getting excited by the transferred recoil energy; and measure- 
ment of 'phonon' excitations generated by nuclear recoils. 
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WIMP Mass [GeV/c^] 

Figure 6.3. Existing strongest exclusion limits on the spin-independent 
elastic WIMP-nucleon cross-section asi versus the WIMP mass m^, pro- 
vided by XENONIOO [146]. The thick blue curve shows the XENONIOO 
limit at 90% C.L. and the dark and light shaded areas depict the la and 2a 
sensitivity of the experiment. The limit is derived with the Profile Likeli- 
hood method (see e.g. section 7.1.1 of the present thesis) where different 
systematic uncertainties are taken into account. For comparison, lim- 
its from some other experiments or analyses are also given: XENONIOO 
(2010) [147], EDELWEISS [148] and CDMS [149]. The 90% C.L. regions 
favored by the two experiments CoGeNT [150] and DAMA [151] are also 
provided. Expected interesting 68% and 95% C.L. regions of the CMSSM 
are shown in shaded gray [152]. Adapted from ref. [146]. 

Constraints provided by direct detection experiments are usually pre- 
sented in terms of exclusion limits on either 'spin-dependent' or 'spin- 
independent' cross-sections of WIMPs and normal nuclei, a'^^ and cr^^ , 
respectively, versus the WIIVIP mass rriy.. There are currently a large 
number of direct detection experiments looking for WIMP signals and 
many others are under construction or planned for construction in the 
near future. The strongest available limits so far have been provided by 
the XENONIOO experiment [146], and, as can be seen in Fig. 6.3, they 
have just started probing interesting regions of SUSY parameter space. 

No positive signal from WIMPs has been detected so far, except for 
the detection of an annual modulation signal by the DAM A/LIBRA ex- 
periment [153, 154] that has not been confirmed by other experiments yet. 
With the upcoming direct detection experiments however, it is expected 
that a substantial fraction of the parameter space for various dark matter 



80 



Chapter 6. Observational constraints on supersymmetry 



models (including SUSY models such as the CMSSM) will be tested. This 
was the main motivation for us in writing Paper III that examines the 
prospects for constraining SUSY models (in the context of the CMSSM) 
with future ton-scale direct detection experiments. Paper III is one of 
the few works that has compared direct detection data (i.e. the number 
of observed events and corresponding recoil energies) directly with the 
model predictions in a full likelihood setup instead of just using available 
exclusion bounds. 

The theory and phenomenology of dark matter direct detection, has 
been reviewed in great detail in Paper III, we therefore do not discuss 
those here. Various experimental issues that should be considered for 
correctly interpreting experimental results, as well as uncertainties in dif- 
ferent nuisance parameters have all been discussed in Paper III. Direct 
detection has also been our particular case of study in Paper IV where 
certain statistical issues in SUSY parameter estimation (as we will discuss 
in the next chapter) have been discussed. 

6.1.3 Indirect detection 

Let us look at the neutralino, our favourite supersymmetric WIMP (see 
section 5.1.4). Neutralinos are 'Majorana' fermions meaning that they are 
identical with their antiparticles. This means that two neutralinos can 
interact and annihilate into other particles including photons, neutrinos, 
antimatter and other types of cosmic rays, as primary or secondary prod- 
ucts. With 'indirect detection' methods one aims to detect such products 
that we receive from the self-annihilation processes (for an introduction, 
see e.g. ref. [23, 24, 25]). 

Similar to the direct detection case, indirect searches usually provide 
their constraints on the WIMP properties in terms of limits on the anni- 
hilation cross-sections as a function of the WIMP mass. The annihilation 
rate is proportional to the square of the WIMP density p^, leading to 
that the best targets for indirect searches are the ones with the highest 
concentration of WIMPs, such as the Galactic Centre or dwarf spheroidal 
galaxies. 

The Galactic Centre could be one of the best choices, because it is 
nearby and potentially contains a large amount of WIMPs. The problem 
is however that it has a highly complex and poorly-understood structure 
with many different components [155, 156]. The dark matter profile of 
the Galactic Centre is also not entirely known [157, 158]. Dwarf galaxies. 
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Figure 6.4. Exclusion limits on the velocity-averaged annihilation cross- 
section (av) versus the WIMP mass , provided by Fermi observations 
of Milky Way dwarf spheroidal galaxies. The red and blue points show 
mSUGRA models that are consistent with all accelerator constraints. The 
red points are models that also give a cosmological dark matter relic 
density equal to the observed value, while the blue ones are models with a 
lower thermal relic density. In the latter case, the neutralinos are assumed 
to be also produced non-thermally and still constitute all of the dark 
matter. The lines show the Fermi 95% upper limits. Adapted from 
ref. [175]. 



on the other hand, are interesting targets because of their high mass-to- 
Hght ratios, a property that reduces the astrophysical background. The 
problem is however that the flux coming from dwarfs is significantly lower 
than the Galactic Centre [159, 160, 161]. 

The other potentially very interesting targets for indirect detection 
are the so-called "unidentified sources". These include all sources of 
radiation in the sky whose astrophysical properties are not identified 
yet. Some of these objects might be small clumps of dark matter whose 
existence is predicted theoretically in some models of the early Uni- 
verse [162, 163, 164, 159, 165]. The problem with these objects is that 
their predicted properties (such as their number, mass and distance to 
Earth) can substantially vary in different theories. 

Finally, one can go even further and look for WIMP signals coming 
from extragalactic sources such as clusters of galaxies [166] or the extra- 
galactic diffuse background (of for example gamma rays) [167]. 
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An indirect detection experiment usually looks for WIMP self- anni- 
hilation signals in a particular annihilation channel, i.e. with particular 
annihilation products. In some cases, for example when photons (that are 
electrically neutral) are being observed, they directly point toward their 
sources, whereas the trajectories of electrons or antimatter particles (such 
as positrons and antiprotons) are easily affected by magnetic fields and 
consequently they do not give us information about the position of their 
sources. 

Currently, various indirect detection experiments are observing poten- 
tial dark matter self-annihilation products coming from different sources 
in the sky through different annihilation channels. Frontiers are: for 
photons, the Large Area Telescope (LAT), aboard the Fermi gamma- 
ray space telescope [168] and several ground-based large air Cerenkov 
gamma-ray telescopes (ACTs) (such as VERITAS [169], MAGIC [170] and 
H.E.S.S. [171]); for electrons, positrons and antiprotons, the PAMELA 
satellite [172], Fermi, H.E.S.S. and some balloon missions such as ATIC 
[173]; and for neutrinos, IceCube [174]. 

So far however, as in the direct detection case, no major signal excess 
has been observed by any of these experiments, and they have therefore 
been able to only provide some exclusion limits on the annihilation cross- 
section versus the neutralino mass. An example of such limits is depicted 
in Fig. 6.4 where limits from observations of Milky Way dwarf spheroidal 
galaxies by Fermi are shown. This shows that like the previous case of 
direct detection, indirect detection experiments are also approaching the 
interesting regions of the SUSY parameter space. 

In Paper I we used Fermi gamma-ray data obtained from observa- 
tions of a particular (and arguably the most interesting) dwarf galaxy. 
Segue 1, to place constraints on the parameter space of the CMSSM. The 
main difference between our approach and the other popular approaches is 
that, like our previously mentioned work in Paper III for direct searches, 
we compared the model predictions with real data directly and in a full 
likelihood setup (see the next chapter). This was arguably the first pa- 
per that applied such a methodology to indirect detection analyses. In 
addition, we included a full treatment of the instrument response func- 
tion (IRE) and its related uncertainties, as well as detailed background 
models. 
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6.2 Collider constraints 

If weak-scale supersymmetry is a correct extension of the SM at high 
energies, it is expected to show up at current and future TeV colhders 
(or accelerators) such as the Large Hadron Collider (LHC) or the pro- 
posed International Linear Collider (ILC). Supersymmetric particles are 
expected to be discovered at these colliders and any positive or negative 
results can place tight constraints on interesting SUSY models. Currently, 
no sparticles have been found and therefore only lower bounds exist on 
their masses that can be used to exclude parts of the SUSY parame- 
ter spaces. In addition, SUSY contributions can indirectly affect other 
measured particle physics quantities including the so-called electroweak 
precision observables (EWPOs) and observables related to some rare pro- 
cesses such as i?-physics observables (BPOs). We review in this section 
some of the most important collider constraints that are widely used in 
SUSY parameter estimation and model selection, including our analyses 
in this thesis. 



6.2.1 Electroweak precision observables (EWPOs) 

The existence of any new physics, including supersymmetry, at close-to- 
electroweak energy scales can indirectly affect various precision observ- 
ables at those energies, and such effects can be probed if high-precision 
experimental data are available [176]. Some of these effects can be large 
such that the absence of significant deviations from the SM predictions 
place strong constraints on models of physics beyond the SM. As far as 
the supersymmetric extensions of the SM are concerned, electroweak pre- 
cision observables have turned out to serve as powerful tools for testing 
the SM and its SUSY extensions by probing indirect effects of SUSY par- 
ticles on those observables. Current precision experimental data, in all 
but only a few exceptional cases, fit the SM predictions very well and this 
implies that the data only put lower bounds on SUSY masses. The small 
deviations from the SM predictions can however be used to favour some 
regions of the SUSY parameter space which explain such deviations by 
higher order corrections that are caused by SUSY particles. 

In addition to the EWPOs that are affected by contributions from 
SUSY particles, the virtual effects of these particles can be detected also 
by measurements of some 'rare' processes. For example, there are some 
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decays that represent flavour-changing neutral currents (FCNCs) and oc- 
cur in the SM only at loop level [177]. This virtually means that such pro- 
cesses happen only very rarely. Examples of these rare decays are some 
B decays such as B ^ Xsj and Bs — )• H''^ IJ'~ ■ In many SUSY models 
however, loop contributions from SUSY particles can be large (compared 
to the SM) and therefore make these processes happen more frequently, 
i.e. with higher rates. The experimental tight constraints on their decay 
rates in these cases can therefore constrain the SUSY models. Obviously, 
the measurements of the observables associated with rare processes can 
also be considered as EWPOs if those measurements are of high preci- 
sion. We however follow the conventions of some authors [176] and discuss 
them under the name i?-physics observables (BPOs) in section 6.2.4. In 
addition, an observable like the anomalous magnetic moment of the muon 
(see below) corresponds both to an EWPO since it has been measured 
very accurately, and to a rare process. We again follow some conventions 
and consider it as an EWPO. Furthermore, quantities like the mass of the 
Higgs boson (if it exists) will be measured with high accuracy in the near 
future. It can therefore be also considered as an EWPOs, but regarding 
the existing limits on its value and the fact that the Higgs has not been 
measured yet, we discuss its constraints on SUSY models separately in 
section 6.2.2. 

Amongst different EWPOs, perhaps the most important ones that 
have been widely used in constraining SUSY models, are the M^-boson 
mass mw, the effective leptonic weak mixing angle siv^O^^, the total 
Z-boson decay width P^, and the anomalous magnetic moment of the 
muon a^. In the rest of this subsection, we briefly describe each quantity, 
its current experimentally measured value and its relevance for SUSY 
parameter estimation. 

Vt^-boson mass mw 

The theoretical prediction for my/ can be expressed as [176] 

m^(l-4^) = ^(l + Ar), (6.1) 



m 



V2G 



where a and Gp are the fine structure constant (calculated at the rriz 
renormalisation energy scale) and the Fermi weak coupling constant, re- 
spectively, mz is the mass of the Z-boson. The quantity Ar in this 
expression stands for all radiative corrections from new physics effects, 
including SUSY (see e.g. refs. [179, 180] and references therein)). 
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Currently the best measurement of the VF-boson mass has been pro- 
vided by LEP [178] as 

mw = 80.398GeV ± 25MeV. (6.2) 

In the case of supersymmetric extensions, the above experimental 
value is so accurate that even two-loop corrections from the superpar- 
ticle effects can be probed. This means that it is quite important that we 
have an accurate theoretical prediction for mw for a given set of SUSY 
parameters. Currently there are some numerical tools for calculating this 
quantity that include two-loop corrections from the SM as well as the 
MSSM. An example is SUSYPOPE [181, 182] (see also ref. [183]). In order 
to compensate various higher-order approximations made in calculating 
the predicted value of mwi one usually considers a theoretical uncertainty 
of about 10 — 15 MeV in addition to the experimental uncertainty given 
above. These types of uncertainties are often used in fitting SUSY models 
to different EWPOs. 



Effective leptonic weak mixing angle sin^ O^ff 

One important group of EWPOs are the Z-boson observables of which 
we briefly discuss here one of the most important ones for SUSY phe- 
nomenology, namely the 'effective leptonic mixing angle' at the Z-boson 
resonance sin^ 6^g. sin^ is a precision observable with high sensitivity 
for probing physics at electroweak scales. It is defined in terms of the 
ratio of the effective vector and axial vector couplings gv and qa, when 
the Z-boson couples to leptons through the vertex ilj^{gv — gAlb)lZ^, 
where I denote leptons [184]. The expression for s\v? then reads 

sin2 0,ff=l(l-i?e^). (6.3) 
^ 9A 

It can be shown that sin^ 0ggf, if calculated at tree level, is approxi- 
mately equal to sin^ & where 6w is the weak mixing angle with sin^ - 
1 — rriy^/m^ in the on-shell renormalisation scheme. The higher-order 
(loop) corrections (from the SM or beyond) are all accommodated through 
the couplings gy and g^ [185, 176]. 

The experimental value of sin^ can be measured at the electron- 
positron colliders through the measurements of different assymmetries 
around the Z-boson peak when QCD effects have been removed [186]. 
Currently the best estimate of the actual value of sin^ is [178] 

sin^ = 0.2324 ± 0.0012. (6.4) 
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Total Z-boson decay width Tz 

Another important Z-boson observable is the 'total decay width' of the 
Z'-boson when it decays to different fermions, i.e. in processes Z — t- //. 
In this case Tz is defined as = Fj, where Fj are the 'partial decay 
widths'. Here only decays to fermions are considered because other decay 
channels are relatively insignificant. The partial decay width Tf can be 
expressed as [176] 



^ 127r 



Ottz^ -i 

Wv? + \9{\\^ - ^) (1 + 6qED) + AF^^^, 



where gy and qa are the effective coupling constants defined in the pre- 
vious subsection, A*"^ is the colour factor being 1 for leptons and 3 for 
quarks, is the mass of the fermion into which the Z-boson decays. 
Sqed '^s some photonic QED correction, and AFg^-,^ denotes the standard 
gluonic QCD corrections plus possibly supersymmetric QCD corrections 
that involve virtual gluinos and quarks. 

Another way of modifying the standard decay width Tz by SUSY pro- 
cesses is when Z-boson can decay 'invisibly' into some new sparticles such 
as neutralinos (if they are sufficiently light). The partial decay width of 
the Z-boson to the dark matter particles with non-zero branching fraction 
is called 'invisible Z width'. In this case, a large deviation of the total 
width Tz from the standard value is expected. 

The currently best experimental estimation for the total decay width 
Tz is [178] 

Tz = 2.4952GeV ± 2.3MeV. (6.5) 



Muon anomalous magnetic moment 

The 'anomalous magnetic moment of the muon' = {g — 2)^/2 is de- 
fined theoretically through the photon-muon vertex function F^^^p in the 
expression [176] 

u{p)Tf,f:,Ap{p, -v\ q)u{p) = u{p) [7pFy(g2) + {p + p)pFM{q^) + ■■] u{p), 

where = — 2m^i<Af (0). At three level, F^i[q^) is vanishing, and non- 
vanishing values are induced by quantum corrections at loop level. Cur- 
rently the best experimental estimation of the actual value of (to eight 



6.2. Collider constraints 



87 



significant figures) fias been provided by the Muon G-2 collaboration 
through an experiment at Brookhaven National Laboratory [187]: 

al'^P = 1165920.80 ± 0.63 x 10~^ (6.6) 

This measured value shows an about 3a deviation from the SM pre- 
diction af,^ = 1165917.85 ± 0.61 x 10"^ [188]. It is common to speak of 
the observed deviation in terms of the quantity (5a^: 

6a^ = a^^^P - a^^^ = 29.5 ± 8.8 x 10"^°. (6.7) 

Supersymmetric models typically well explain the above deviation. To 
show this, let us just present here a simple example where only one-loop 
corrections are considered and also all SUSY particles are assumed to 
have equal masses of the value MsusY- The SUSY contribution to the 
muon anomalous magnetic moment in this case reads [189] 

^WsuSY 

with tan/3 and signfi defined for example as the CMSSM free parameters 
in section 5.3.2. Assuming > 0, and lOOGeV < MsusY < 600GeV 
(depending on the value for tan/3), it can be seen that SUSY effects can 
easily provide the required contribution to the theoretical value of in 
order to reconcile it with observations. It is interesting to also notice 
that the SUSY contributions in some cases may instead deteriorate the 
situation by generating very large or even negative values for a^^^"^. This 
should however be considered a plus for the muon anomalous magnetic 
moment as an observable because it can place strong constraints on the 
SUSY parameter space. 



6.2.2 Experimental bounds on the Higgs mass 

We saw in section 5.1.3 that the predicted value for the MSSM lightest 
Higgs boson h satisfies the condition m/i < mz ~ 91 GeV at tree level. 
This is much lower than the currently strongest lower bound on the Higgs 
mass 114 GeV provided by the LEP direct search [78]. Loop corrections 
can however increase rrih up to 135 GeV. 

Imposing the lower limit 114 GeV can consequently exclude regions 
of SUSY parameter space that cannot provide the required loop correc- 
tions to avoid the limit. To be more precise, searches for the MSSM 
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Particle mass 


Commonly-used LLs 


Conservative LLs 
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> 100 GeV 
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> 95 GeV 


> 73 GeV 




> 87 GeV 


> 73 GeV 


mi, 


> 94 GeV 


> 43 GeV 




> 95 GeV 


> 65 GeV 




> 95 GeV 


> 59 GeV 


ruq 


> 375 GeV 




m~g 


> 289 GeV 





Table 6.1. Experimental lower bounds on the Higgs and supersymmetric 
particles [192]. 

Higgs bosons at LEP have been performed using two production pro- 
cesses e~^e~ — )■ hZ and e~^e~ — )■ hA where A is the CP-odd Higgs bo- 
son. Both processes are mediated by s-channel Z-boson exchange [190] 
and the cross-sections are proportional to the quantities sin^(/3 — a) and 
cos^(/3 — a), respectively. Here /3 is as usual the ratio of the up- type to 
down-type Higgs VEVs and a is the mixing angle between the two CP- 
even Higgs bosons h and H and the Higgs interaction eigenstates and 
H^. One can therefore be more stringent and impose the experimental 
bounds on SUSY models depending on the values for — a (see also 

Tab. 6.1). 

In addition to the LEP lower bounds on the Higgs mass, the CDF 
and DO collaborations have recently excluded an additional mass range 
for the Higgs, but at higher values, using pp collisions at the Fermilab 
Tevatron at ^/s = 1.96Tey. Their combined results exclude the region 
158GeV < ruh < 173GeV at 95% C.L. [191]. These semi-upper bounds 
however do not really put any new constraints on the MSSM parameters 
becasue such high Higgs masses cannot be acheived within the MSSM. 

6.2.3 Experimental bounds on sparticle masses 

As in the case of the Higgs boson, negative results from the current collider 
searches for SUSY particles places lower limits on their masses. These 
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limits can then be used to exclude the points in the supersymmetric pa- 
rameter space that predict masses violating the lower bounds. 

Currently the best lower limits on the sparticle masses have been 
provided by the Particle Data Group [192]. These limits are however 
not completely model-independent and particular (usually CMSSM-like) 
SUSY models are often assumed in extracting the limits. In addition, 
the constraints on each sparticle mass is often dependent on its difference 
with the mass of the lightest supersymmetryc particle (LSP) which for 
the case of the Particle Data Group limits is assumed to be the neutralino 
Xi- Commonly used limits on the Higgs and sparticle masses (and the 
more conservative ones applied under specific conditions) are summarised 
in Tab. 6.1. 

6.2.4 S-physics observables (BPOs) 

In this section we describe four major observables involving S-mesons [177] 
and their decays into other particles that are widely used in the global fits 
of SUSY models to the experimental data. Those are the branching ratios 
BR{h — )■ 57), Bg — )■ Ijl^ , Bu — >• tu and the mass difference between B^ 
and Bg. There are B-physics observables other than these four, but they 
are of less interest and we do not consider them in this thesis. 



Branching ratio for B Xsj 

Perhaps the most important i?-physics observable is the experimentally 
measured value of the decay rate for the flavour changing process B — )■ 
Xgj. In the SM, the main loop contributions are from the VF-boson and 
top quark, and the SM prediction for the branching ratio of the process, 
that we denote by BR{b — )■ 57), is in excellent agreement with experiment. 
Additionally, in the MSSM, mainly chargino/stop and charged Higgs/stop 
loops contribute to the branching ratio [193] and the contributions from 
neutralino loops are comparatively small [194]. The two types of SUSY 
contributions can be large individually, but they can interfere destruc- 
tively with each other giving rise to a value that is not significantly differ- 
ent from the SM prediction. This means that the observable BR{b — )• 57) 
can place stringent constraints on the SUSY parameter space. The con- 
tributions from SUSY are particularly large when charged Higgs bosons 
are light and the parameters fi or tan /3 have large values. 
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The current world-average experimentally measured value for the branch- 
ing ratio is given by the Heavy Flavor Averaging Group (HFAG) as [195] 

BR{h s-i)exp = (3.55 ± 0.221°:°^ ± 0.03) x lO""^, (6.9) 

whereas the SM contribution to the process at next-to-next-to leading 
order in QCD is theoretically predicted to be [196, 197, 198, 199] 

BR{b s-i)sM = (3.28 ± 0.23) x 10"^, (6.10) 

which shows only a la difference from the experimental value. 



Branching ratio for Bs M M 

Another important S-physics observable in SUSY phenomenology is the 
branching ratio for the flavour changing process — )• /U"*"//". The SM 
prediction for this observable is [200] 

BR{Bs fJ^'^fi')sM = (3.42 ± 0.54) x 10~^ (6.11) 

which is well below the experimental upper bound 

BR{Bs ^+Ai")exp < 5.8 X 10"^, (6.12) 

a value given by CDF II data at 95% C.L. [201]. 

The MSSM however predicts that if tan /3 is large, neutral Higgs 
bosons can contribute to the branching ratio and enhance it by several 
orders of magnitude. This can easily violate the above experimental up- 
per bound and therefore, the observable BR{Bs — )• can impose 
important constraints on the MSSM parameter space. 



Branching ratio for tv 

The last -B-physics process that we discuss here is B^ — )• ti/, for which 
the SM prediction reads [202] 

BR{B^^ ^ Tu)sM = - 4)^11^.1^^. (6.13) 

87r 

Here Gi? is the Fermi weak coupling constant, ms and m,- are the B- 
meson and r masses, respectively, Jb is the so-called S-meson decay 
constant, Vub is the u^-element of the Cabibbo-Kobayashi-Maskawa quark 
mixing matrix and tb is the i3-meson lifetime. 
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The experimental value of the branching ratio for this process has been 
mainly measured by the two experiments BELLE [203] and BABAR [204] 
and an often-used world-average value is [205] 

BR{Bu Tv)exp = (1.31 X 0.49) x lO""^. (6.14) 

This value does not quite agree with the SM prediction and the devi- 
ation can be expressed as 

= 0.93 ± 0.41. (6.15) 

In the MSSM, the main contributions to BR{Bu — ?• Tiy) are from the 
direct-exchange of a virtual Higgs boson that decay into ti^. In scenarios 
with minimal flavour violation (such as the CMSSM and NUHM) one can 
show that [205] 



BR{Bu tv)mssm 



^ ml tan^/3 ^ 
BR{Bu^Tu)sM L ^m2^±^l + eotan/3J ' ^ ' ' 

where eo is the effective coupling that parametrises the correction to the 
down-type Yukawa coupling from the gluino exchange (see below for its 
mathematical expression) , is the 5- meson mass and mij± is the mass 
of the charged higgs boson. 



Bs — Bs mass difference 

Finally, another interesting i?-physics quantity that is used in comparing 
new physics predictions (including SUSY) with experimental data is the 
mass difference Bs — Bg that is often denoted by AM^^. This quantity 
is interesting because the frequency of oscillation between particle and 
antiparticle states of the neutral i?-meson (which is measurable exper- 
imentally) is proportional to AM^^. In addition, this quantity can be 
predicted theoretically both in the SM and the MSSM and therefore be 
used to constrain SUSY parameter spaces. 

AMs^ has been measured experimentally by the CDF collaboration 
to have the value [206] 

(AA/sJexp = 17.77 ± 0.12ps-^ (6.17) 

The deviation from the SM prediction is [133] 



(6.18) 
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In the MSSM, the major additional contribution to AM^^ comes from 
the exchange of neutral Higgs bosons [205]. The full MSSM prediction 
compared to the SM one is expressed as [202] 

{AMbJmssm _ ^ 647r sin^"' mbmbmsmbjeY tan^ /?)^ 

{AMbJsm ~ am^5o(m2/m2,) [l + (eo + ey)tan/3]2[l + eotan/3]2 

Here the masses of the bottom and strange quarks and rus are cal- 
culated in the modified minimal subtraction renormalisation scheme MS 
and at the renormalisation scale rub- Sq is a function with the form 

= 4(1 -x)2 - 203^- ^'-^'^ 
eo and ey have the forms 



eo — - 7, H2 — Y ' — Y ' ey — - ^ H2 



and are the effective couplings parametrising the corrections to the down- 
type Yukawa couplings, is the usual supersymmetric Higgs mass term, 
and yt and at are the Yukawa and soft trilinear couplings that involve 
stopts. H2 is a function with the following form: 

/ X xlogx ylogw 

H2{x,y) = - rf -+ 6.20 

{l-x){x-y) [l-y){y-x) 
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Statistical subtleties 



One principal objective of the present thesis has been to provide some 
additional strong support for the fact that SUSY models (even in cases 
where the number of free parameters is dramatically reduced) exhibit 
rather complex structures in their parameter spaces. In Paper II and 
Paper IV we have essentially tried to further demonstrate that constrain- 
ing models using existing experimental data is under no circumstances an 
easy or straightforward task, and therefore care must be taken when par- 
ticular statistical frameworks are being used in such analyses. Indeed, 
we have shown that current data do not sufficiently constrain the model 
parameters in a way completely independent of the employed statistical 
framework and scanning techniques. We have extensively discussed some 
of the important statistical issues in the papers, concluding for example 
that the current scanning techniques may yet have some distance to go 
in this respect. The statistical frameworks and scanning strategies are 
also discussed in detail in the papers, in particular Paper II and Pa- 
per IV. In this chapter we only go into some general descriptions of the 
frameworks and algorithms and refer the reader to the papers for more 
details. 



7.1 Statistical frameworks 

In order to make any meaningful statistical inference about a theoretical 
setup, such as our favourite SUSY models, when its predictions are com- 
pared with experimental data, one needs to first make a decision about 
the statistical framework and formalism that should be used. This for 
example tells us how to make statistically significant statements about 
the model parameters. 
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There are two commonly-used but fundamentally different approaches 
to statistics that are based on entirely different interpretations of 'prob- 
ability' (for a detailed discussion, see e.g. ref. [207]): the so-called "fre- 
quentist" and "Bayesian" statistics (or inference). The first approach 
deals with relative frequencies, while the second talks about subjective 
probabilities. Usually, the former is employed for assigning statistical 
errors to measurements, whereas the latter is used to also quantify sys- 
tematic uncertainties. The basic difference between the two frameworks 
is that in Bayesian inference one is interested in probabilities of some 
model parameters to have particular values when some data are given, 
whereas in frequentist statistics the only meaningful quantity to work 
with is the probability of observing particular data when a specific set 
of model parameters is assumed. In this section we look into the defi- 
nitions and fundamental ingredients of the two statistics and introduce 
some statistical measures that we used in the papers. 

7.1.1 Frequentist inference: profile likelihood 

A frequentist accepts the most common interpretation of probability that 
defines it as a "limiting relative frequency". Let us assume for example 
that we measure a statistical variable through a process that is repeat- 
able. In addition, suppose that we are interested in a particular outcome 
of the measurement that we denote by O. The frequentist then defines 
the probability corresponding to the outcome O, denoted by -P(O), as the 
fraction of times that O occurs when we repeat the measurement pro- 
cedure an infinite number of times. In reality however, no experiments 
can be repeated 'an infinite number of times', and consequently the prob- 
abilities defined this way always remain hypothetical. One of the main 
tasks of a frequentist statistician is then to provide an estimation of the 
exact probabilities based on a finite set of observed data. Having this 
given, the next step would then be to establish an appropriate method 
to investigate the compatibility of a particular theoretical model, which 
predicts the probabilities, with the experimental data and also place some 
constraints upon the model parameter space. 

Let us see how this is usually performed in a frequentist setup by 
assuming that n measurements of a random variable O are made. We 
additionally suppose that the probability density function (PDF) for the 
variable O (denoted by p{0)) is not known a priori. Our aim is now to 
'infer' properties of p{0) based on the set of the n observed data points 
di, ...,dn (that we denote by D). One may in particular be interested in 
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constructing some functions of D so as to appraise different character- 
istics of p(0). In most cases, a hypothetical form of the function p(0) 
is available in terms of some undetermined parameters — {0\^ ...^O^f^. 
This means that the PDF can now be shown as p(0; 0). The objective 
is then to estimate the values of the parameters 0, and this is done by 
constructing particular functions of the data points D. Such functions 
are called 'estimators' and are often denoted for the parameter set by 
0. One requirement for the estimator is that it converges to the 'true' 
(whose actual values are, and may forever remain, unknown) when 
the number of data points n becomes infinitely large. The entire process 
through which the actual model parameters are being estimated from 
the empirical data points D is called parameter estimation or parameter 
fitting. 

The estimators 0(di, d^) are themselves new random variables. 
The reason is that the estimators are functions of the measured values 
D = {di, dn), and obviously, if we repeat the experiment several times, 
each time we obtain in general different values for the measured quan- 
tities. As a result, the estimators 0(D) also receive different values in 
different repetitions of the measurement. One can now define another 
PDF, say q(Q;Q), that corresponds to the statistical distribution of the 
estimators 0; this PDF is called sampling distribution. Studying vari- 
ous statistical properties of estimators using the corresponding sampling 
distributions is another major task of any statistical analysis. This in par- 
ticular includes estimating different experimental errors and theoretical 
uncertainties associated with the estimators. 

One important point about the estimators is that there are various 
ways to construct them for a particular set of model parameters. There 
are however some properties that are said to be desirable for an estima- 
tor. This for example include small (ideally zero) biases and variances. 
Unfortunately one cannot always achieve all the desired properties for 
an estimator simultaneously (an example is the trade-off between bias 
and variance [207]). There are however some methods that provide es- 
timators with optimal properties (e.g. with reasonably small bias and 
variance). Examples are the methods of 'maximum likelihood' (ML) and 
'least squares' (LS) that are arguably the most popular ones in various 
fields of scientific data analysis. Here we are in particular interested in 
the ML method and continue our discussions with describing the method 
and some of its interesting properties. 
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The central ingredient of the ML method (as it is clear from its name) 
is the likelihood function C that is defined as 



where p{0; G) is the PDF according to which our random variable (pre- 
viously shown as O) is distributed, and di,...,dn are the corresponding 
values we have obtained in n repetitions of the experiment. We addition- 
ally assume that the functional form of p(0; 0) is determined although 
the values of the parameters = (^i, ...,9m) are not known. C, if consid- 
ered as a function of the data points D = di, dn, is precisely the joint 
PDF for D. However, one can consider this quantity instead as a function 
of the parameters 0; the variables D are assumed to be fixed in this case. 
Having the likelihood function defined, the ML estimators = 6i, ...,9^. 
corresponding to the parameters = 9i,...9m are then defined as the 
values of the parameters at which C{Q) is maximised. 

As we mentioned earlier, the ML estimator for a parameter of the 
model is itself a random variable and this brings us to the point where we 
should find an appropriate way of quantifying uncertainties in the fitted 
value of the parameter. The simplest way is to calculate and report the 
variance (or the standard deviation) of the estimator, e.g. by simulating 
several experiments with the same number of samples in each as we had 
in the actual measurement and then calculate the best-fit values (i.e. ML 
estimators) in each of those experiments. Then it would be possible to 
work out an estimator for the variance of the obtained parameter esti- 
mators. This procedure however, is not always adequate enough, and the 
statistical uncertainty of a measurement must be given in some other more 
sophisticated ways. This is commonly performed by introducing the 'con- 
fidence intervals' for the parameters in the following way: Assume that 
the experiment is repeated a large number of times, and an interval [a, h] 
contains the true value of a parameter in a particular fraction of times, 
say, 7. The interval [a, b] is then said to be a confidence interval at a 
confidence level (C.L.) 7 [207]. 

For the case of one parameter, the definition above is fine. Now the 
question is: 'can one generalise the concept of confidence interval to the 
case of m parameters = (^1, 0^) in a straightforward way?' One way 
of doing this might be to introduce an m-dimensional confidence interval 
[a, b] in which a = (ai, am) and b = (61, ...,bm) in such a way that each 
[ai , bi] is defined separately as a confidence interval for the i^^ parameter 
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9i . This recipe is however not what is often used in the statistical inference 
(it is computationally difficult [207]). One instead constructs a so-called 
'confidence region' in the parameter space. This is defined as the region 
which contains with a particular probability a point in the parameter 
space corresponding to the true values of the parameters 0. Clearly, 
the form of this region is in general different from the aforementioned 
n-dimensional confidence interval [a, b]. It turns out that if n — )• oo, this 
region approaches an ellipse for m = 2 parameters and an m-dimensional 
hyperellipsoid in general [207]. 

Now, the question is 'how to find in practice such confidence regions 
for a given model at hand with unknown parameters for which we have 
found ML estimators'. One can show that both of the joint PDF g(0|0) 
for the estimator © = (6*1, Om) and the likelihood function £(0) become 
Gaussian in the large sample limit (i.e. when n — )• oo). That is, contours of 
constant g(0|0) and £(0) are hyperellipsoids (ellipses in two dimensions) 
in ©-space and ©-space respectively, centred about the true parameters 
© and the ML estimators ©, correspondingly [207]. It can be shown that 
for this likelihood function of Gaussian form, the regions in the ©-space 
specified by the contours 

ln£(©) = ln£™,,-^ (7.2) 

2 

of the log-likelihood function ln£(©) (or —\) are nothing but the previ- 
ously defined confidence regions with confidence levels 7. The quantities 
Q^^m are the 'quantiles' of orders 7 of the distribution, i.e. 

f^2{z;m)dz = -f, (7.3) 

where f^2{z;m) indicates the distribution for m degrees of freedom. 
Values of the Q-y^m, for la and 2a C.L.s (i.e. 7 = 68.3% and 95.4%) and 
m = 1, 2, 3, 4, 5 fitted parameters, are given in Tab. 7.1. 

The prescription given by Eq. 7.2 for determining the confidence re- 
gions is completely true only in the case of a Gaussian likelihood. How- 
ever, this method is also employed for non-Gaussian functions (including 
the SUSY likelihood function) as a reasonably appropriate approximation 
to the actual confidence regions. In Paper IV, we discuss one main issue 
with this approximation for the cases where the model parameter space is 
large and highly complex, namely the 'statistical coverage' problem. We 
do not discuss the issue here. 
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Q 



m = 1 m = 2 m = 3 m 



= 4 m = 5 



68.3% (la) 1.00 2.30 3.53 4.72 5.89 
95.4% (2cj) 4.00 6.17 8.02 9.70 11.30 

Table 7.1. Values of the quantile Q-y,m for 68.3% and 95.4% (i.e. la and 
2a) confidence levels and for m = 1, 2, 3, 4, 5 fitted parameters. 

Now, all one needs to do in order to determine the best-fit points and 
the associated uncertainties is to map the likelihood function of the model 
given the experimental data. From the plotting point of view, it is a very 
straightforward procedure if there are only one or two fitted parameters in 
the model with one-dimensional (ID) and two-dimensional (2D) likelihood 
functions, respectively. However, for the higher-dimensional functions, 
one should have a good recipe for summarising the statistical inference 
(i.e. the best-fit points and the errors) for each parameter separately, or 
in a 2D plane for one parameter versus the other. A nice way of doing this 
is to make use of the "profile likelihood" (see e.g. ref. [208] and references 
therein) that is defined e.g. for one parameter 9i as 



where C{Q) is the full likelihood function. This definition can easily be 
generalised if one is interested in two-, three-, or higher dimensional profile 
likelihoods, although it is often sufficient to calculate ID and 2D profile 
likelihoods only. 

This definition simply means that a frequentist eliminates unwanted 
parameters by maximising the likelihood along the hidden dimensions. In 
other words, the profile likelihood is nothing but the likelihood function 
of the reduced set of parameters with the unwanted parameters at their 
conditional ML estimates. Now the interesting point about the profile 
likelihood recipe is that, the approximate confidence regions can be set 
using exactly the same prescription of Eq. 7.2 just as in a standard 
fit, but now with the quantiles Q^y^i and (57,2 for the ID and 2D plots, 
respectively. 

7.1.2 Bayesian inference: marginal posterior 




(7.4) 



A fundamentally different way of interpreting probabilities is that of 
Bayesian or subjective statistics. Contrary to the previous framework 
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of frequentism where possible outcomes of a measurement constitute the 
sample space, here one instead talks about the hypothesis space, which 
consists of hypotheses. These are statements that are either true or 
false [207]. Therefore, one of the crucial features of the Bayesian frame- 
work is that a probability can be assigned to a hypothesis, say H. This 
is not possible under the frequentist framework, where a hypothesis can 
only be rejected or not rejected. This probability P(]HI), in the Bayesian 
context, is interpreted as the degree of belief that the hypothesis HI is 
true. This means that, the Bayesian probability can in particular be as- 
sociated with the values of unknown parameters in a theoretical model 
of interest where these parameters themselves are considered as random 
variables. These subjective probabilities reflect our degree of confidence 
that the parameter values reside in certain intervals. Again, a probability 
for an unknown parameter is not meaningful within the context of the 
frequentist statistics. 

The above description of the Bayesian inference then leads to a more 
practical definition of it, namely as a statistical inference in which evi- 
dence or observations are used to update or to newly infer the probability 
that a hypothesis may be true (for an introduction to general applica- 
tions of Bayesian inference in physics, see e.g. ref. [209], and for reviews 
of its applications in cosmology, see e.g. refs. [210, 211, 212, 213]. The 
name Bayesian comes from the frequent use of "Bayes' theorem" in the 
inference process, which in fact forms the basis of Bayesian statistics. 

Assume again (analogous to the previous subsection) that the model 
at hand is parametrised by m unknown parameters G = {9i, ...,9m) and 
there are some experimentally provided data D = (di, ...,dn) which are 
supposed to be used for constraining the parameter values. Bayes' theo- 
rem then reads 



Here p{Q), the so-called 'prior PDF' (or simply the 'prior') and usually 
shown as 7r(6), represents our degree of belief (or the state of knowledge) 
that the parameters G are the true values before the consideration of the 
data, p{D\Q) is nothing but the previously introduced likelihood function 
C(Q), i.e. the probability under the assumption of the specific values of 
G for the model parameters, to observe the data D. p{Q\D) is called the 
'posterior PDF' (or simply the posterior), refiecting the probability of G 
as being the true values after seeing the data D. Finally, the quantity in 
the denominator, i.e. p{D), is called the 'Bayesian evidence' (or simply the 




p{D\e)pie) 
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evidence), which for the purpose of constraining parameters of a model, 
is in fact nothing but a normahsation constant, i.e. independent of 0, and 
is often simply dropped. 

It is important to reahse that Baycsian statistics does not offer any 
fundamental rule for choosing priors for a model's parameters; this should 
be done based on other theoretical considerations. After choosing a partic- 
ular prior, one can see how the degree of belief (or the state of knowledge) 
changes, i.e. how the prior is updated to the posterior, when the data are 
used. This upgrading information is provided by the model's likelihood. 

In Bayesian statistics, all of our knowledge about the parameters Q 
is contained in the posterior PDF p{Q\D). Very often however the pa- 
rameter space of the model at hand is multidimensional and one needs to 
summarise the important characteristics of the multidimensional posterior 
in a practically appropriate way. One usually starts with introducing an 
estimator, which is often taken to be either the expectation values of the 
parameters E[Q] corresponding to the posterior PDF (i.e. the 'posterior 
mean'), or the values of Q that maximise it (i.e. the 'posterior mode'); the 
latter coincides with the ML estimator discussed previously if the prior 
7r(G) is taken to be a constant (i.e. for a flat or linear prior). In this case 
the posterior p{Q\D) is proportional to the likelihood C{Q). As the next 
step in making statistical conclusions about the model parameters based 
on the posterior PDF, one should find an appropriate way of quantifying 
uncertainties, similar to what we did by introducing the confidence inter- 
vals and regions in a frcqucntist framework. Perhaps the most natural 
way is then to construct a so-called 'credible region' in the parameter 
space that contains a certain fraction of the total probability given by 
integrating the posterior over the whole space. This can be simply done 
for example for a 2-dimensional parameter space by drawing a contour 
such that, say, 68.3% of the total posterior falls inside the contour. This 
contour then describes a credible region at a confidence level of 68.3% (or 
la). 

Again, like what we did in the previous case of frequentist statistics, we 
should have a good recipe for calculating and showing uncertainties about 
the estimated value of one of the fitted parameters, say 6i, or drawing 
the relevant contours in the 2-dimensional plane of one parameter, say 
6i, versus the other, say 9j. This can be performed in a natural way 
if the corresponding ID and 2D posterior probability density functions 
for those parameters are available in some manner. Paying attention to 
the fact that the full posterior p(Q\D) is nothing but the joint PDF of 
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all the parameters (which are treated just as random variables in this 
framework), the probability densities for a fewer number of parameters 
are just 'marginal' densities, defined as the ones obtained by marginalising 
(i.e. integrating over) the unwanted (or hidden) parameters. That is, for 
a ID posterior PDF of the parameter ^j, this procedure gives a "marginal 
posterior" as 



A 2D posterior is defined in an analogous manner. 

Now that we have constructed these reduced PDFs, credible intervals 
and regions can be defined in the same way as if we have had only one or 
two fitted parameters from the beginning. 

By looking at the mentioned characteristics of the two frameworks of 
Bayesian and frequentist statistics, we realise that, from a practical point 
of view, they are rather different in (1) the ways one defines estimators 
and associated uncertainties and (2) the recipes for discarding unwanted 
parameters in the statistical inference (i.e. marginalization of the posterior 
PDF p{Q\D) in the marginal posterior prescription, and maximisation of 
the likelihood function L(0) in the profile likelihood prescription). It 
is therefore quite clear that the results of the two inferences might not 
coincide in general, even if the prior PDF 7r(B) is taken to be constant, 
i.e. when p[Q\D) is proportional to L(G). This is especially true if the 
model likelihood has a complex dependence on the parameters (i.e. not 
just a simple Gaussian form) and experimental data are not sufficiently 
available. Of course, in the case of the large sample limit, they (must) 
give similar results (since both the statistical measures become almost 
Gaussian in this case); this is why both methods are commonly used in 
data analysis. 

Finally, it is worth stating some aspects of perhaps the cornerstone 
of Bayesian statistics, i.e. including prior functions in the final inference. 
One practically interesting consequence of this is that it gives a powerful 
way of estimating how robust a fit is. That is, strong dependence of the 
posterior on different priors actually means that the data are not sufficient 
or accurate enough to constrain the model parameters. In other words, the 
posterior in this case is dominated by the prior rather than the likelihood 
function. It can be shown also that if a fit is robust in this language, 
the Bayesian and frequentist methods should identify similar regions of 
the parameter space at any particular confidence level. This robustness 
issues are for example investigated in [214] for the CMSSM with different 
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choices of priors, where it is shown that the currently available data are 
not yet sufficiently constraining to determine the best-fit values of the 
parameters independent of the priors. Therefore, one should inevitably 
decide which approach to use in the statistical analysis of the model. 

In presenting and interpreting our direct and indirect dark matter 
detection results in Paper I and Paper III, we have employed both 
Bayesian and frequentist approaches and constructed the corresponding 
confidence and credible regions. In Paper II we have been particularly 
interested in frequentist statistics and the ability of existing scanning 
techniques in correctly mapping profile likelihoods in SUSY global fits. 
Finally, in our Paper IV, we have focused on the issue of statistical 
coverage in SUSY parameter estimation which directly affects frequentist 
inference; we have however analysed the coverage for both profile likeli- 
hoods and marginal posteriors. 

7.2 Scanning algorithms 

Even if one decides on a particular statistical framework to work in, and 
also chooses the appropriate statistical measures that match the frame- 
work, there are other related issues that need to be addressed in any 
statistical analysis of a theoretical model; SUSY models are no exception. 

One of these issues is the limitations of the available computational 
methods in correctly sampling a large and complex parameter space. It 
is for example crucial within the frequentist framework to map the profile 
likelihoods appropriately so as to be able to construct correct confidence 
regions and intervals. In addition, having a good estimate of the globally 
maximum value for the likelihood function is of extreme importance in 
this framework, because the results of the statistical inference strongly 
depend on that value. In a large, complex and poorly-known parameter 
space however, usual scanning methods such as grid or purely random 
scans can easily fail or otherwise take an extremely long time to generate 
useful results. This has therefore led many phenomenologists to employ 
more sophisticated scanning algorithms in their analyses. 

Most of the existing advanced scanning techniques that are also widely 
used in SUSY phenomenology, are based on 'Markov Chain Monte Carlo' 
(MCMC) algorithms (see e.g. ref. [215]) or their different variations such 
as MultiNest [216, 217], which is an algorithm based on nested sam- 
pling [218, 219] (We do not describe these methods here and instead 
refer the interested readers to the given references for introductions and 
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detailed descriptions). The structures of these algorithms are such that 
they exceptionally match the requirements for an appropriate scanning 
technique that is optimised for Bayesian statistics. These techniques have 
however been also used in frequentist inference, for example to map profile 
likelihoods and corresponding confidence regions. Even though there are 
currently various ongoing efforts in improving such techniques for these 
purposes, there are still some severe issues that have not been properly 
answered by them yet. 

In Paper II and Paper IV, we have studied some of these issues with 
the statistical analysis of SUSY models (in the context of the CMSSM) 
that we think stem from the imperfection of the utilised statistical scan- 
ning algorithms in particular cases. In Paper II, we in addition propose 
an entirely different type of scanning techniques for SUSY parameter es- 
timation, based on Generic Algorithms [220] (GAs; for a classic introduc- 
tion, see e.g. ref. [221]; for recent introductions, see e.g. refs. [222, 223, 
224]; for a modern treatment, see e.g. ref. [225]), that can help in a com- 
plementary way the state-of-the-art and powerful Bayesian methods such 
as MultiNest. The specific version of GAs that we use in our analysis is 
described in detail in Paper II and we therefore do not detail it here. 
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Summary of results 



We have so far given an introduction to the field of supersymmetry phe- 
nomenology with an emphasis on different strategies in comparing various 
viable supersymmetric models with observational data. This in particu- 
lar includes advanced statistical techniques and innovative numerical al- 
gorithms in analysing the complex parameter spaces of SUSY models. 
In this chapter, we give a summary review of our contributions to the 
field that have been presented in the papers included in the thesis (see 
Part III). We also discuss some future prospects for constraining super- 
symmetric models and parameters. 

8.1 Major achievements 

Paper I [226] and Paper III [227] deal with the observational constraints 
that direct and indirect searches for particle dark matter can place upon 
the parameter spaces of minimal SUSY extensions of the SM. We restrict 
our investigations to the CMSSM. 

In Paper I, we analyse the CMSSM parameter space when Fermi 
gamma-ray observations of dwarf galaxy Segue 1 are considered. We per- 
form this in a global-fit framework where likelihoods from Fermi data 
are combined with those from the cosmological measurements of dark 
matter relic density (section 6.1.1), electroweak precision observables in- 
cluding the anomalous magnetic moment of the muon (section 6.2.1), 
collider bounds on the masses of SUSY and Higgs particles (sections 6.2.2 
and 6.2.3), and i?-physics observables (section 6.2.4). Our objective is 
to know which regions in the CMSSM parameter space are favoured by 
Fermi data alone or in the presence of other existing constraints. Our em- 
ployed Fermi likelihoods are constructed from both the observed energy 
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spectrum of gamma rays from Segue 1 and the spatial distribution of the 
gamma-ray photons at and around Segue 1. We also convolve the theoret- 
ical spectrum with the instrumental energy dispersion and point spread 
function. We use the Fermi data from 9 months observations of Segue 
1 and also extrapolate our analysis to 5 years of observations assuming 
that no excess events will have been observed after this time (similar to 
the actual 9 months data). 

We scan over the CMSSM parameter space using MultiNest as imple- 
mented in SuperBayeS [228, 229, 230, 231, 215, 214, 232] (available from 
ref. [233]). SuperBayeS is a numerical global-fit package for exploring the 
parameter space of the CMSSM and finding both Bayesian credible and 
frequentist confidence regions when the model's predictions are compared 
with different types of experimental data. It reconstructs the marginalised 
posterior PDFs, as well as the profile likelihoods although the employed 
scanning techniques (i.e. MCMCs and MultiNest) are both optimised for 
Bayesian statistics. For physics calculations, it uses different packages: 
SOFTSUSY [234] (available from ref. [235]), DarkSUSY [132] (available 
from ref. [137]), FeynHiggs [236, 237, 238, 239] (available from ref. [240]), 
Bdecay and micrOMEGAs [138] (available from ref. [139]). 

The number of total CMSSM points in our final set of samples is 
~ 5.5 X 10^. Our results show that some points in the CMSSM parameter 
space are ruled out when only Fermi data are used. These are models 
that give very large annihilation cross-sections {{(Jv) > 10~^^cm^s~^ in 
the best case where 5 years data are used and a large boost factor at 
the source is assumed) and low neutralino masses. These models are 
excluded because no gamma-ray excess signal is detected from Segue 1. 
When we combine Fermi observations with other constraints, in particular 
the observed value of dark matter relic density, we see that most of the 
CMSSM points excluded by Fermi are already disfavoured by those other 
constraints (see Fig. 8.1). 

The analysis of Paper I presents the first direct inclusion of con- 
straints from indirect dark matter searches in global fits of supersym- 
metric models to experimental data. It is based on the full likelihood 
construction for the experiment, i.e. the Fermi data points are directly 
used in the statistical analysis instead of the usual use of exclusion limits 
provided by the experiment. This makes the paper a significant contribu- 
tion to the field and gives a concrete example of what should be done for 
all types of experimental data in order to provide statistically consistent 
constraints on SUSY and other beyond-the-SM models. In addition, the 
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Figure 8.1. Favoured values of the neutralino annihilation cross-section 
and mass for the CMSSM, when: (1) only existing constraints from 
the dark matter relic density, electroweak precision observables, collider 
bounds on the Higgs and sparticle masses, and B-physics observables are 
included (upper panels), (2) 9 months of Fermi observations of the dwarf 
galaxy Segue 1 are added (central panels), and (3) no excess events are 
observed from Segue 1 after 5 years of Fermi observations (lower pan- 
els). The left panels show frequentist profile likelihoods (where yellow 
and red indicate 68% and 95% confidence regions, respectively) and the 
right panels give Bayesian marginalised posterior PDFs (where solid blue 
contours depict 68% and 95% credible regions). Best-fit points and poste- 
rior means are indicated by black crosses and dots, respectively. Adapted 
from Paper I. 



110 



Chapter 8. Summary of results 



paper is of particular interest because it was the first work in which Fermi 
gamma-ray data have been used in searches for dark matter signals. 

In Paper III, we perform a similar analysis where data from direct 
dark matter searches are used in constraining SUSY parameters. Here, 
instead of working with real data, we are mainly interested in constraints 
provided by the future generation of direct detection experiments, namely 
the ones with ton-scale target materials. We want to know how well those 
experiments can detect a supersymmetric WIMP with particular masses 
and cross-sections and in each case, what constraints they put upon the 
fundamental properties of the underlying theory. 

We investigate these by choosing a few benchmarks in the CMSSM 
parameter space with masses and cross-sections that cover an interesting 
range of values. We generate some synthetic data (i.e. the number of 
events detected by an experiment and the corresponding recoil energies) 
and use those data to scan the CMSSM parameter space and obtain the 
favoured regions. 

The experiments considered in this work are ton-scale extrapolations 
of three existing experiments: CDMS (Cryogenic Dark Matter Search) 
[241], XENON [242] and COUPP (Chicagoland Observatory for Under- 
ground Particle Physics) [243]. These experiments are expected to coop- 
eratively explore large fractions of WIMP masses and scattering cross- 
sections, both spin-independent (SI) and spin-dependent (SD). The two 
experiments CDMS and XENON have the capability of measuring the re- 
coil energies produced in the detector by WIMPs that interact with target 
nuclei. These experiments are however not very sensitive to SD scatter- 
ing. COUPP on the other hand, although only measures the number 
of events above a threshold nuclear recoil energy and does not measure 
event energies associated with the event, has a significant sensitivity to 
SD scattering. COUPP is therefore expected to break some degeneracy 
in the CMSSM parameter space by providing both SI and SD measure- 
ments. Obviously this ability increases when the number of observed (SD) 
events becomes larger. For each experiment, we assume 1000 kg- years of 
raw exposure and energy resolution, energy range and efficiency similar to 
the present day version of the experiment. We also assume backgrounds 
at target levels. In a similar way as in Paper I, our analysis is based 
on reconstructing full likelihoods for the experiments, namely, we use the 
experimental data points (i.e. the number of events and recoil energies) di- 
rectly in the likelihood constructions rather than only using the exclusion 
limits. 
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For each benchmark point, we scan over the CMSSM parameter space 
by means of MultiNest again as implemented in SuperBayeS. We however 
do not combine the direct detection hkelihoods with other constraints 
since we are only interested in what one can gain from direct detection 
data. Another important constituent of the paper is that we take into 
account existing uncertainties on hadronic matrix elements for neutralino- 
quark couplings (important in calculating theoretical cross-sections), as 
well as on halo model parameters (important in calculating differential 
recoil rates). We marginalise (maximise) over these nuisance parameters 
in order to find marginal posteriors (profile likelihoods). 

We show that when the likelihoods from all three experiments are com- 
bined, the primary characteristics of the neutralino, namely its mass and 
cross-sections, can be determined with high certainty if the neutralino sits 
at low masses and high cross-sections (first and second rows in Fig. 8.2). 
In particular, it can be seen from our results that when COUPP is added, 
it can break degeneracies in the parameter space and therefore substan- 
tially help us pin down the actual WIMP properties. This can be achieved 
in cases where a substantial fraction of the signal event rate comes from 
SD interactions of WIMPs and nuclei. The results also indicate that the 
uncertainties on the halo parameters may have relatively large effects on 
the reconstruction of both Bayesian credible and frequentist confidence 
regions. The uncertainties on the hadronic matrix elements on the other 
hand do not have significant effects on the credible and confidence regions. 

As far as the constraints on the CMSSM parameters are concerned, 
our analysis indicates that using direct detection data alone does not place 
strong constraints upon the values of those parameters. We can see this 
from the third and fourth rows in Fig. 8.2 where la and 2a contours are 
given in the mo — mi/2 and Aq — tan /3 planes for the benchmark WIMP 
with low mass and high cross-sections: direct detection experiments only 
determine the gaugino mass parameter mi/2 with relatively high precision. 
This is mainly due to the fact that this parameter has a strong correlation 
with the neutralino mass that can be measured almost accurately by direct 
detection experiments. 

Perhaps the bottom line of Paper III is yet another indication that 
no single type of experiments is able to present conclusive information 
about the fundamental properties of a dark matter model (e.g. supersym- 
metry in our case) . There is always a high degree of degeneracy between 
different parameters of the model that cannot be broken by only one type 
of experimental data. In other words, in order to test and characterise 
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Figure 8.2. Upper two rows: Reconstruction of the mass and scat- 
tering cross-sections of the neutrahno for a CMSSM benchmark when 
hkehhoods from ton-scale versions of three direct detection experiments 
CDMS, XENON and COUPP are combined and synthetic data are used 
in the scan. Lower two rows: Constraints on the CMSSM parameters from 
the combination of the experiments. Left panels: 2-dimensional profile 
likelihoods (yellow and red indicate 68.3% and 95.4% confidence regions, 
respectively). Right panels: 2-dimensional marginal posteriors (inner and 
outer contours represent 68.3% and 95.4% confidence levels, respectively). 
Black dots and crosses show the posterior means and best-fit points, re- 
spectively, and benchmark values are marked with green stars. Adapted 
from Paper III. 
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e.g. weak-scale SUSY as a valid extension of the SM that also provides 
"the solution" to the dark matter problem, a combination of informa- 
tion from complementary experiments such as direct, indirect and collider 
searches would be required. It is however important to note here that fu- 
ture generations of direct detection experiments are highly promising in 
this endeavour and constitute one of the main strategies for identifying 
the nature of dark matter as well as characterising various models beyond 
the SM that include a dark matter candidate. 

In Paper II [244] and Paper IV [245], we discuss the ability of 
existing scanning techniques in correctly exploring the parameter spaces 
of SUSY models in order to provide favoured regions when the models' 
predictions are compared with observational data. We are in particular 
interested in how well they find the confidence regions in a frequentist 
approach. 

As we have discussed in section 7.2 of this thesis, the most powerful 
algorithms that are currently employed in SUSY parameter estimation 
are MCMCs and nested sampling. The structures of these techniques are 
such that both of them are optimised for Bayesian statistics, i.e. they can 
accurately enough map the posterior PDFs and accordingly the marginal 
posteriors. In the absence of advanced scanning methods that are specifi- 
cally designed for frequentist statistics, MCMCs and nested sampling are 
nowadays also exploited in the latter case, i.e. to map profile likelihoods in 
supersymmetric analyses. It is known that (see section 7.1) profile like- 
lihoods and marginal posteriors provide very similar results given that 
high-statistics data are available (i.e. when the effects of priors are dom- 
inated by likelihoods). This also implicitly means that in these cases, 
the scanning algorithms that are optimised for one statistical framework 
can be appropriately used for the other. However, in many realistic cases 
(at least in SUSY phenomenology), the model at hand exhibits a very 
complex parameter space and experimental data are not strong enough. 
As a result, Bayesian credible and frequentist confidence regions are in 
general very different from each other. As we have remarked explicitly in 
Paper I, we think that in order "to gain as complete a picture as pos- 
sible of the preferred regions in an insufficiently-constrained parameter 
space like the CMSSM", both statistical measures should be considered. 
That is, rather than arguing that one statistical approach is more 'cor- 
rect' than the other, we think that additional insight can be obtained 
if one considers both for analysing any given parameter space. This is 
however a useful strategy only if one can accurately map the confidence 
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and credible regions. This is precisely where our points in Paper II and 
Paper IV stand: The methods that work properly for one approach do 
not necessarily work appropriately for the other. 

In Paper II, we propose a new scanning technique based on Genetic 
Algorithms (GAs) (see section 7.2) that is optimised for frequentist profile 
likelihood analysis of complex supersymmetric models. We use existing 
cosmological and collider constraints on SUSY predictions (the same data 
as in Paper I except for Fermi observations of Segue 1) in a global-fit 
setup and analyse the CMSSM as a testbed. We implement the GA 
code PIKAIA [246] (available from ref. [247]) in SuperBayeS and compare 
our results with those of MultiNest (with the standard configuration as 
is employed e.g. in ref. [214]) when the same set of experimental data is 
used. 

Our results are quite surprising: the best-fit point found by GAs has 
a substantially higher likelihood value compared to the one found by 
MultiNest. In addition, many new CMSSM points with high likelihood 
values show up in our scans. These together dramatically impact the 
inferred confidence regions for the CMSSM parameters, as well as some 
other derived quantities such as the neutralino mass and cross-sections 
(see Fig. 8.3). Our investigation therefore indicates that the conven- 
tional Bayesian scanning techniques, including nested sampling (which 
is arguably the best existing one), may give unsatisfactory results in a 
frequentist context, except they become appropriately modified or recon- 
figured. See e.g. ref. [248] where a more appropriate MultiNest configu- 
ration for profile likelihoods is presented and its results are in excellent 
agreement with the results of Paper II. This new configuration however 
requires a significantly larger computational effort. We make the point 
here that even though our application of GAs for the profile likelihood 
analysis of the CMSSM appears to be quite successful, it is not flawless. 
The algorithm is primarily designed to find the global maximum of a 
complex function and is not optimised to accurately 'map' the function 
around the global maximum. This is why the confidence regions found in 
our results (see Fig. 8.3) look somewhat noisy. Perhaps this displeasing 
property of the algorithm can be remedied by combining it with other 
scanning algorithms that are designed to map a function in the vicinity 
of a given point. 

Paper II also has some physics implications. One of the most in- 
teresting ones is the observation that contrary to the findings from some 
MCMC analyses, our best-fit point lies in the focus point region (i.e. at 
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Figure 8.3. Best-fit points in the CMSSM parameter space (upper pan- 
els), as well as in the a'p^ —m^o (central panels) and {av)—m-^o (lower pan- 
els) planes found by Genetic Algorithms (left panels) and MultiNest (right 
panels). The blue and green regions represent iso- likelihood contours that 
correspond to la and 2a confidence regions in the MultiNest scan. The 
dotted circles and triangles indicate the global best-fit points in the GA 
and MultiNest scans, respectively, and the dotted square shows the best-fit 
co-annihilation point in the GA scan. The GAs find more high-likelihood 
points, while MultiNest better maps the likelihood around the moderately 
good fits that it finds. The big circle in the lower left panel indicates a 
section of the stau co-annihilation region found by GAs and missed in 
previous scans. Adapted from Paper II. 
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high mo) rather than the stau co-annihilation region. We however think 
that this discrepancy is hkely to come from different physics codes used 
in calculating different observables rather than from the differences in the 
employed scanning techniques. In addition, our results uncover the exis- 
tence of a section of the stau co-annihilation region at large mo that seems 
to have been commonly neglected in all previous scans (see Fig. 8.3). The 
paper also presents some other implications for Higgs and sparticle masses 
at the LHC, as well as quantities important in direct and indirect searches 
for dark matter. 

We go on to further investigate the ability of scanning algorithms in 
correctly providing the confidence regions for SUSY models. This time, 
in Paper IV, we look at a central requirement for properly constructed 
confidence regions, namely the 'statistical coverage'. This means that a 
region corresponding to a specific confidence level must include the 'true' 
values of the model parameters at the stated confidence level when the 
experiments are repeated infinitely many times. We study the coverage 
for the CMSSM when its parameter space is explored with Bayesian scan- 
ning techniques, in particular MultiNest. We perform the analysis for a 
simple case where only constraints from a direct detection experiment 
are imposed on the CMSSM parameter space. Two CMSSM benchmarks 
are chosen and their corresponding parameters are used as hypothetically 
true parameters. We then generate 100 sets of synthetic direct detection 
data in each case and scan over the CMSSM parameters for each set. In 
order to obtain the degree of coverage, we construct one-dimensional con- 
fidence intervals and count how many times (out of 100) a true parameter 
falls within the intervals. We assess the coverage when two types of pri- 
ors are imposed on the parameter space: (1) flat priors on all parameters 
of the model and (2) logarithmic priors on the scalar and gaugino mass 
parameters mo and mj^/2- For comparison, we also examine the coverage 
for Bayesian credible intervals although we do not expect proper coverage 
for these cases. 

Our results indicate both over- and under-coverage that for some cases 
strongly vary when benchmarks or priors are changed (see Fig. 8.4 for an 
example of scans with under-coverage). The possible reasons for the ob- 
served poor coverage can be twofold: (1) The profile likelihood is mapped 
correctly by the scanning algorithm, but its validity as a proper approxi- 
mation to the full Neyman construction (for details, see e.g. Paper IV or 
ref. [249, 250]) of the frequentist confidence intervals breaks down for com- 
plex parameter spaces such as the CMSSM. In this case more sophisticated 
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Figure 8.4. Two-dimensional profile likelihoods (left) and marginalised 
posterior PDFs (right) for one typical coverage scan over the CMSSM 
parameter space with significant under-coverage. The inner and outer 
contours in each panel represent 68.3% (Ict) and 95.4% {2a) confidence 
levels, respectively. Plots are shown in mo-mi/2 and a^' — m^o planes. 
Black dots indicate the benchmark values which in this case lie outside 
the la and 2a regions. 



methods such as the 'confidence belt' construction [250] could significantly 
improve the coverage. The problem with these alternative techniques is 
that they are rather difficult to implement numerically. (2) The profile 
likelihood is still a good approximation but the scanning algorithm has 
failed to correctly map it. This can clearly lead to poor coverage because 
the confidence regions and intervals are not properly constructed. Our 
analysis in Paper IV indicates that in our particular case of study, while 
option 1 may have played some role, it is more likely that a substantial 
fraction of the over- and (especially) under-coverage comes from the fact 
that the employed scanning algorithm is not optimised for the profile like- 
lihood analysis (i.e. option 2). For another coverage study of the CMSSM 
using a different type of experimental likelihood, see ref. [251]. 
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8.2 Outlook 

The coming decade in particle physics and cosmology will definitely be 
an exciting period in the history of mankind's great endeavour for un- 
derstanding Nature at its most fundamental level. There are currently a 
large number of theories and models as extensions of the Standard Model 
that equally well describe all known phenomena while give different de- 
scriptions of physics at high energies. These theories are all waiting for 
experimental verification, and it would be of extreme interest to know 
whether any of them could be confirmed or excluded. Fortunately, many 
of these new physics frameworks have several predictions that could be 
tested observationally. With the advent of many new experiments and 
computational techniques in recent years with unprecedented power and 
precision, the prospects for verifying and constraining these theories are 
quite promising. 

Most of the new physics models possess rather complex structures 
with large numbers of free parameter. This means that any phenomeno- 
logical analysis that aims to properly compare new physics theories with 
observations should in principle be composed of four main elements: (1) 
a concrete theoretical model whose predictions for any set of free param- 
eters are calculable, (2) relevant observational data corresponding to the 
theoretical predictions with the associated uncertainties, (3) the statis- 
tical framework and formalism for the analysis, and (4) an appropriate 
scanning algorithm. 

The present thesis centres upon one of the most interesting new physics 
frameworks, namely weak-scale supersymmetry, and presents a number 
of powerful strategies and approaches for comparing its predictions with 
observations. A concrete SUSY model is studied (the CMSSM), various 
experimental data are employed (from direct and indirect searches for 
particle dark matter as well as collider constraints), different statistical 
frameworks are used and advanced scanning techniques are utilised. 

Our results (in Paper I and Paper III) indicate that even though 
the current data do not give conclusive information about the validity of 
even the simplest versions of supersymmetric models and the preferred 
values of their parameters, the situation could dramatically change in the 
near future. We explicitly investigate the prospects for dark matter di- 
rect detection experiments, and our results show great promise for the 
next generations of these experiments in characterising supersymmetric 
models. Future indirect detection experiments are also potentially very 
encouraging especially because there are various targets to look at and 
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different annihilation channels to use. However, perhaps the landmark 
in these directions would be the LHC results from searches for physics 
beyond the SM, including SUSY. Analyses similar to ours could be per- 
formed using data from the LHC (see e.g. refs. [252, 134, 152]). Ideally, all 
different sets of cosmological, astroparticle (direct and indirect searches) 
and collider constraints should be eventually combined in a statistically 
consistent global-fit setup where are sources of uncertainties are taken 
into account. 

As far as the statistical frameworks and scanning algorithms are con- 
cerned, our investigations in Paper II and Paper IV imply that there 
are several subtleties that should be considered when working in any 
framework and with any algorithm. Supersymmetric models in particu- 
lar, exhibit quite non-trivial parameter spaces that are constrained differ- 
ently when different statistical formalisms are employed. This also makes 
it difficult for a scanning method to properly provide favoured values of 
the parameters in all statistical frameworks. Methods optimised for one 
framework usually give unsatisfactory results for the other. Some of the 
existing scanning techniques are highly powerful, efficient and relatively 
fast, but their widely- used versions are usually optimised for Bayesian 
statistics. There have been some successful attempts to reconcile them 
with frequentist framework, but at the price of (significantly) increasing 
the computational effort. The latter is indeed a very important point. It 
should not be forgotten that the full phenomenologically interesting mod- 
els (such as the MSSM in SUSY extensions) possess very large parameter 
spaces. Constraining such models requires a sufficiently fast scanning 
technique that provides the results in reasonable times. We therefore 
think that the efforts in constructing efficient scanning techniques for 
both statistical frameworks should still continue. 

Finally, we should make the point here that what we have done so far 
have all been in the context of SUSY 'parameter estimation', namely that 
a SUSY model is given (the CMSSM in our case) and we attempt to fit it 
to the data so as to find the most favoured values for the free parameters 
of the model. The statistical frameworks we have used can however be 
employed equally well in a slightly different context, namely 'model se- 
lection'. As we broadly discuss in this thesis, there are several models of 
supersymmetry that are constructed based on different assumptions and 
motivations. Perhaps a natural way to go, when more constraining data 
become available, would be to compare different models in a statistically 
consistent way so as to see which models are completely excluded by ex- 
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periments or at least which ones are more favoured. Ideally, this should 
not be restricted to SUSY models and other physics-beyond-the-SM theo- 
ries should also be analysed. We do not go through a discussion of various 
techniques and strategies in this direct and only note that such attempts 
have already been started in the community although most of the current 
results appear to be inconclusive (see e.g. ref. [253]). 
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